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Time 
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of gravity in its deformed position 
Blade cross-sectional mass center-of -gravity 
offset from the elastic axis, shown in Fig. 28. 
Positive for c.g. before E.A. 

Blade cross-sectional aerodynamic center offset 
from elastic axis, shown in Fig. 28. Positive 
for A. C. before E.A. 


Coefficients in perturbation expansion of X 

H 

Coupling term in flap equation 

Value of X at 0 = 6 
H c 

Function used in phase-plane investigation 
of condition for suppressing secular terms 


Coefficients in the perturbation expansion 
° fY G 

Coupling term in lag equation 
Value of Y„ at 8 = 0 



Y(S,n) 


Function used in phase-plane investigation 
of condition for suppressing the secular terms 


Greek Symbols 


W° 3 ' a 3 


Quantities defined in Appendix G 
Angle of attack of the whole rotor 


3 


f*l ' (?2 ' 


Inclination of feathering axis with respect to 
the hub plane measured in a vertical plane, 
(angle of built-in coning, shown in Fig. 27a 
Quantities defined in Appendix G 

Angle of preconing, shown in Fig. 27b 


Y 


Y F 


Locke number (y = 2p ^ ■ *) 

*b 

mth inplane bending mode 

Flight path angle with respect to horizontal 


6 


Quantity used in phase-plane investigation 
of condition for suppressing secular terms 


e 

e 


Perturbation parameter 

Perturbation quantity in uniformly valid expansion 
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Symbolic quantity having the same order of 
magnitude like the displacements v and w 


?/ n 


n il ,r, 12' n 2l' n 22 


r^d) 


n 1 (D 


n SF. ,n SL. 
1 1 


Coordinates attached to the blade cross section 
origin located at elastic axis, ri coincides with 
the blade chord, 5 perpendicular to the blade 
chord, shown in Fig. 28. q also used in the 
phase-plane investigation for suppressing the 
secular terms. 

Coefficients in the perturbation expansion of 

damping 

H, at x =1 
to 

Quantity defined in Appendix F 

kth flapwise normal bending mode 

Structural damping coefficients defined in 
Appendix B ' 
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Pitch angle measured from x,y plane 
Critical value of collective pitch 
Coefficients in perturbation expansion of 0 
Constant part of collective pitch 
Cyclic pitch component, multiplied by cos 
sin ip* respectively 

Constants defined by Eqs. G.80 and G.81 

Quantities used in the perturbation method 
defined in Appendix G 

Value of inflow ratio in hover 

Approximate expression defined in Eq. D.9 
Inflow ratio, induced velocity over disk, 
positive down, nondimensionalized with 
respect to OR 

g(x) first cosine component of X 
Constant part of X^ c 

Steady state part of X constant over disc 

Value of X at critical condition 
o 

Coefficients in perturbation expansion of X 


Advance ratio 

Coefficients used in the perturbation expansion 


of the advanced ratio 



Quantity of order one, or less, used to 

characterize the closeness of 0 ) to r or 1 

c 2 

Coefficients of the characteristic equation 
used in Section 10 

Coordinate used in the phase plane investigation 
of the condition for suppressing the secular term 
Defined in Appendix H 

Quantities used in the derivation of the conditions 
for suppressing the secular term 

Quantity associated with limit cycle amplitude 
used in the investigation of the condition for 
suppressing the secular term 

Value of p at ij; = 0 ; 

o 

Value of p for stable limit cycle 
-Density of air 

Mass of blade per unit area used in Appendix K 


Solidity ratio 


Quantity used in phase plane investigation of 
condition for suppressing the secular term 



x 2 »x 2 

X X 'X 2 


Quantities defined in Eq. 5.95 and page 

Amplitudes of oscillation about equilibrium 
position for flap and lag, respectively 


W ••• *n 


Azimuth angle of blade (ip = ftt) measured 

from straight aft position 

Time scales used in multiple time-scale 

expansion 


u> 

c 


(u> ) 

c 


1,2 


00 . f U T . 

Fi Li 


FIN' LIN 


Ct) ,U) 
F10' L10 


Flutter frequency 

_ Flutter frequencies used in the complete 
linearized flap- lag-pitch problem 
Natural frequency of the ith flap or lag mode 

Nonrotating flap and lag frequencies, 
respectively, nondimensional with respect to ft 
Unperturbed value of first flap and lag 
natural frequency, respectively 


ft 


Speed of rotation 
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Special Symbols 

n 


( ) ' 


* 

( ) 


( )* 

* 

n 


( > R , ( ) r 


( n ) 

% 




-i 


[ j 


Nondimensionalized quantity unless otherwise 

stated, lengths associated with elastic 

bending properties, nondimensionalized with 

respect to all others with respect to R; 

frequencies with respect to mass properties 

with respect to I, 
b 

Differentiation with respect to x q 

Differentiation with respect to tp 

Complex conjugate of the quantity in brackets 

Derivative of the product in the bracket 

Subscripts R and I denote, respectively, the 
real and imaginary part of the appropriate 
quantity 

The symbol beneath a quantity denotes a 

vector or a matrix 

Denotes the inverse of a matrix 

Transpose of a matrix [ ] 
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SECTION 1 
INTRODUCTION 


1.1 General 

A helicopter blade in forward flight is exposed to a severe aeroelastic 
environment. Periodic, unsteady airloads act upon the blade due to a combina- 
tion of forward flight and rotation of the blade. Strong inertia, including Cori- 
olis forces, due to the relatively large speed of rotation further complicate the 
problem. Therefore, the aeroelastic stability and the response of the blade 
are of extreme importance for both flutter calculation and vibration level 
estimation. Vibration level predictions are required for both the linear and 

the nonlinear range of blade motion in order to evaluate the fatigue life of 

* 

the blade and the blade supporting structure. The nonlinearities are those 
arising from the inclusion of moderately large deflections in the inertia and 
aerodynamic loading terms. Therefore, during the last thirty years, a con- 
siderable amount of work associated with helicopter blade dynamic and aero- 
elastic problems has been done. 

A good review and an elementary description of the various dynamic and 
aeroelastic problems associated with VTOL vehicles, in general, and helicopter 
blades, in particular, have been given recently by Loewy in Ref. 32. A con- 
siderable amount of the work done up to 1964, in the general area of heli- 
copter blade dynamic response and flutter, has also been reviewed with a 
considerable amount of detail by Bielawa (Ref. 35) . In this report, no at- 
tempt will be made to repeat the reviews given in Refs. 32 and 35. The only 
references cited will be those pertinent to the problem being treated. 

This report is divided, essentially, into three distinct and almost 
independent parts: 

(a) The first part, composed of Sections 1 through 7, deals with 

the somewhat controversial problem of flap-lag-type instability 
of torsionally-rigid, hingeless helicopter blades. This type 
of instability is analyzed in the nonlinear range of blade 
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motion in both hover and in forward flight. The main purpose 
of this part is the identification of the physical mechanism 
of the instability and the determination of the respective 
roles of the forcing function, parametric excitation, and non- 
linear coupling in affecting the coupled flap-lag response. 

(b) The second part, composed of Section 8, presents a con- 
sistent derivation of the general nonlinear equations 

of motion for a hingeless helicopter blade, having flap-lag 
and torsional degrees of freedom. The torsional degree of 
freedom is represented by elastic root torsion; thus dis- 
tributed torsion and built-in twist are not treated. All 
second-order terms (in terms of the displacements) in the 
flap and. lag equations are retained, while third-order terms 
are neglected. In the feathering equation, some important 
nonlinear effects are included by retaining third-order terms. 

The main purpose of this section is to emphasize the various 
approximations involved in obtaining the elastic, inertia, 
and aerodynamic loads. 

Hopefully, these equations will serve as a starting point for 
future work in this field. 

(c) The third part. Sections 9 through 11, is devoted to the 
investigation of the stability of coupled flap- lag-pitch 
blade motion. This treatment is limited to the case of 
hovering flight. Due to the novel aspects of the various 
effects included in the equations of motion, both divergence 
and flutter boundaries had to be. obtained from the linearized 
equations of motion. 

The main purpose of this part is to illustrate how the 
stability boundaries obtained in the first part of this 
report are affected by the addition of the torsional de- 
gree of freedom. 

Each part of this report will have its own introductory section in which 
the pertinent literature will be surveyed. 
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1.2 Introduction to the Flap-Lag Stability Problem 


1.2.1 Brief Review of Past Work 

In the first part of this report, the flap- lag-type instability of 
torsionally-rigid hingeless blades in the nonlinear range of blade motion 
will be treated. This problem was first treated by Young (Ref. 1) with a 
restrictive analytical approach. Modal equations of motion were obtained, 
but the numerical results were evaluated for a blade represented by a cen- 
trally-hinged, spring-restrained, equivalent model. Young concluded that the 
triggering mechanism of the flap- lag-type instability is the lag degree of 
freedom. 

Hohenemser (Ref. 2) treated the same problem, using a somewhat uncon- 
ventional numerical integration scheme. Due to the various approximations 
made in Ref. 2, the results presented there are of a qualitative nature. 

A good treatment of the linear stability of the blade in" hover has 
been made recently by Ormiston and Hodges (Ref. 3). In this work, both the 
centrally-hinged, spring-restrained and modal-elastic representation of a 
hingeless blade were used. Stability boundaries for the linear case were 
obtained. 

The linearized equations of motion in flap lag, at high values of ad- 
vance ratio, were treated by Hall (Ref. 4) . Using a rigid, offset hinged, 
spring-restrained representation of 'the blade, multivariable Floquet theory 
was applied to investigate the stability of blade motion. Reverse flow ef- 
fects were included. The trim conditions associated with the variation of p 
were vaguely mentioned without specifying what they were or how they were 
satisfied. The primary purpose of the Ref. 4 work was the investigation of 
the blade response in the presence of a simple control system; a stability 
investigation in the aeroelastic sense was not intended. 

The flap-lag instability and response for both articulated and hinge- 
less blades was also treated by Elman (Ref. 5) . Reverse flow and stall ef- 
fects were included in this work. Unfortunately, the description of the theo- 
retical part of this work was very brief. 
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The transient nonlinear flap-lag motion of a fully-articulated rotor 
blade was also treated by Jenkins in Ref. 38. Reverse flow effects were in- 
cluded; two-dimensional stall and compressibility effects could also be ac- 
counted for, but were not used in the numerical calculations. The equations 
of motion were solved using direct numerical integration. The physical ex- 
planation of the mechanism of the instability was not attempted due to the 
purely numerical nature of this work. 

The perturbation method in multiple time scales (Ref. 6) has been 
first applied to the nonlinear flap-lag problem by Tong (Refs. 7 and 8). 

1.2.2 Objectives of the Present Study 

A preliminary study of the nonlinear flap-lag probem has been made in 
Ref. 41 using numerical integration. The main purpose of Ref. 41 is to gain 
some physical insight into the problem before applying the perturbation method. 

In the present report, a consistent system of equations representing 
the flap-lag motion of a hingeless elastic blade with moderate nonlinearity 
has been derived. Modes are assumed and the spatial variables are eliminated, 
using Galerkin's method. The resulting system of equations is solved using 
the perturbation method in multiple time scales (Refs. 7 and 8). For some 
cases, results were also obtained by direct numerical integration. 

The effect of forward flight is studied with the additional require- 
ment that the helicopter should be in trim, i.e., the thrust coefficient is 
kept constant during the increase in the advance ratio; thus the effect of 
forward flight on a fixed configuration can be obtained. The effects asso- 
ciated with the trim requirement were disregarded in the other works dealing 
with forward flight, except, possibly, in Ref. 4. 
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SECTION 2 


BASIC ASSUMPTIONS 


The rotor blade can be considered to be a thin flexible beam attached 
to the hub at its inboard end and free at its outboard end. Figure 1 describes 
the geometry of the problem. , 

The large deflections will have only a small effect on the tension in 
the blade due to elastic effects, because one of its ends is free (in fact it 
can be shown that this is a third-order effect in terms of deflections) ; there- 
fore, a linear treatment of the elastic restoring forces can be considered ade- 
quate. Such a theory has been derived by Houbolt (Ref. 10). It is very similar 
to the usual engineering beam theory. 


It is assumed that the blade is initially straight, with its elastic 
axis coincident with the x-axis. The blade is torsionally rigid. It has a 
pitch setting of 0 and it can bend in both the y and the z directions. The 
cross section of the blade is assumed to be symmetrical about the major prin- 
cipal axis. 


A quasi-steady two-dimensional aerodynamic strip theory is used as 
justified by Miller (Ref. 11) and apparent mass effects are neglected. This 
means that in the usual unsteady aerodynamic expressions, Theodorsen's lift 
deficiency function C(k) = 1. Stall, compressibility, and reversed flow ef- 
fects are neglected; thus, the aerodynamic load is applicable to moderately 
large advance ratios. 


Periodicity and spanwise variation of rotor inflow is restricted to the 
first cosine component: 


'~X- + 'Xic J(X) COS f 


( 2 . 1 ) 


In deriving the equations of motion, an x, y, z coordinate system (see 
Fig. 1) rotating with the shaft of the helicopter and attached to the blade is 
used. 
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SECTION 3 


THE EQUATIONS OF MOTION 


3.1 Brief Derivation of the Equations of Motion 

A brief derivation of the equations of motion is given below. A more 
complete and detailed derivation is given in Sec. 8. The equations of dynamic 
equilibrium of a blade undergoing only bending’ in flap and lag can be taken 
from Ref. 10: 


J[(£l)yCos z e + (Ej)^c n z eJ^ f [( £ l)i ~( £ l)j] 

. r&s^ej- i f t(x 0 ) i&i 7=/> z 


5 Cr>& Ceo & d v 
dx£ 


fL / [( El)z ' ( El)y 1 ^ j- [( Fl) SCn'e + 

2Xo L L ° * 

(fl), coc^Jg - Te, cose l . £ [Tlx. )£ ]- > 


(3.1) 


where p ^ and p^ , which are given in Eqs . 3.3 and 3.4, include the pertinent 
aerodynamic and inertial forces. 

For present purposes, it is assumed that- the elastic axis, area cen- 
troid, aerodynamic center, center of gravity, and feathering axis of the blade 
are all coincident. The offset between these points in the cross section of 
the blade is important only when one considers the torsional degree of freedom 
also. 

For moderate angles of pitch setting, the elastic coupling due to pitch 

setting was neglected. This effect was included in the coupled flap-lag pitch equa 

tions treated in the last part of this study. Reference 3 illustrates, in detail, 

the importance of this effect, which is strongly stabilizing. As pointed out in 

Appendix L, a consistent treatment of this effect requires a more sophisticated 

treatment of the elastic mode shape than the one used in this study. Therefore, 
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the elastic coupling will be neglected in order to be able to study conveni- 
ently the destabilizing aerodynamic and inertial flap-lag coupling terms. In 
a practical sense, this assumption is a limiting case of reduced elastic 
coupling which can occur for hingeless rotor blades with nonuniform stiffness 
distribution. 


~dX> 






(3.2) 


The loading terms in the z and y direction, with nonlinearities accurate up to 
third order in displacements can be written as 


Pi s Li - 


** 


Jl m W 


l. 


* _ 

w Jl 


SF 


(3.3) 


Py~ Lu - tn Jl 1 [ V -(co + v) +2.U J - Jl V 


(3.4) 


p=- II * mJL z (Xo+ei + T-v) 

* 3Xo 


(3.5) 


The last terms in Eqs. 3.3 and 3.4 represent viscous-type of structural damping. 

The boundary conditions for a hingeless blade can be taken as 

.3 


vv(0,t)= {Q,i) _ d W ( i, i ) _ 2 W (l,t) 

dx 0 L Tlx* 


v(o,t}- 'dv(o l t) _ ^ V ( 1,1) _ lPv( i,t) 


(3.6) 




ixf 


Next, in order to apply Galerkin's method, the displacements u, v, w of the 
beam are expressed in terms of the normal modes of the beam in flap and lag 
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V- - l Y (x») hmti) 

u * -i 

Equation 3.9 represents the shortening effect or inboard movement of a 
mass point on the blade due to bending, under the assumption that the blade is 
inextensible. 

Note, that whenever repeated indices are used in this report, the summa- 
tion convention is always implied, unless otherwise stated. 

The substitution of Eqs. 3.7 through 3.9 into Eqs. 3.3 through 3.5 and 
3.2, together with the application of Galerkin's method on the resulting system 
of equations, yields after a considerable amount of algebraic manipulation 



(3.8) 


(3.9) 


(3.10) 

kc -+ g bLf - kc + Mu a j u - [l Scmr - z(Mr) c„ r ] 

- 2 f 'V ^ 2* ~ e ° " ^rj b h ^ Jx ° 

^ 2j M 

(3.11) 

In obtaining Eqs. 3.10 and 3.11, the boundary conditions, Eqs. 3.6, and 
the orthogonality relations for rotating beams, given in Appendix A have been 
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used. The various quantities M^, ••• etc., used above are defined in 

Appendix B and represent generalized masses. 

Next, the aerodynamic loading terms will be evaluated. The loading 
term in the z-direction can be obtained from Miller and Ellis (Ref. 11) or 
from Bisplinghoff (Ref. 12) . 

Li - a,fbRV T (U T 6-U p ) (3.12) 

The aerodynamic load in the y-direction (see Fig. 2) 

L r - f atK[V p (V T e-U,)+ C_£ Ur 1 ] (3 . 13 , 


where the velocities U , U can be written as 

p T 

Up- uiR cosy/^y- wSl 

''bXo 


(3.14) 


l/ T = dlR (jUSinf +X )-+Jlv (3.15) 

The expressions for the aerodynamic loads given above represent aero- 
dynamic loads, for moderate angles of collective pitch 0. The nonlinearities 
originating from these expressions are due to the retention of the second- 

order terms due to large deflections in the appropriate relations for U U , 

2 2 T P 
U T » and U p . 

The reference plane used in evaluating relations in Eqs. 3.12 through 
3.15, for a hingeless blade, is the hub plane. Therefore, in an exact formu- 
lation, the terms associated with the cyclic pitch variation should also appear 
in these equations. Keeping in mind that the primary aim of this report is 

the evaluation of the nonlinear effects associated with large displacements, 

* 

and the determination of the respective roles of parametric excitation and 
forcing due to forward flight, it was decided to treat the effect of cyclic 
pitch variation only in a later part of this investigative effort which will 
also include the torsional effects. 


■*Note that in this study, the term parametric excitation, stands for the effect 
of the time dependent coefficients in the equations of motion. 
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The substitution of Eqs. 3.7 through 3.9, 3.14 and 3.15 into Eqs. 3.12 
and 3.13 will yield the appropriate expressions for and L^. A further inte- 
gration will give the generalized aerodynamic loads required for Eqs. 3.10 and 


3.11: 


JL 


4M, 


* ~a 




ci Xo 


(3.16) 


A 


Li = 


- I 1 
di l u 



(3.17) 


The complete expressions for A^ . , A are given in Appendix C where the 

1 20 1 19 

various flap coefficients F . . . F and lag coefficients L . . . L are also 

defined. 


Using Eqs. 3.10, 3.11, 3.16, and 3.17, the final form of the general 
equations of motion, for an arbitrary number of modes, can be written in 
compact form 

^ r +f iK t + ^ ^ & - 2. p itn & k* * a* 

i = 1, 2, . . . , N (3.18) 


MlC he + 2su hi + Mu <*> Lc hi = [ 2. S Cfn r - IfMjfhmrJh^ hr 

- x fa fa ' e ° Cc + Au ' 


i - 1, 2, 


M 


(3.19) 


3.2 Equations of Motion for the Two- Mode Case 

The solution for the system of equations derived in the previous section 
for an arbitrary number of modes is very difficult and complex. Therefore, in 
this report only one elastic mode in each degree will be considered. 

The static equilibrium condition in hover (y = 0, with all time deriva- 

o o 

tives equal to zero) denoted by g^, h , will be used as the natural equilibrium 
state of the system. 
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(3.20) 




Ml, 



+ Qa 

cl 


L']-\C.J 


(3.21) 


In the following, the lower indices on the various L 


i 


and F 1 


coefficients 


will be dropped. By expressing 

% * X* 

hi~ L, + ^2 


(3.22) 


and considering only the case A- — =0, Eqs. 3.18, 3.19 and 3.20 through 3.22 

X c 

become 

x t 4 ($ b1 4 C-f )-£, 4 [^Fi F c lJ* <-Oif(F - £ F %* )J^1~ 

[(l,eF i -> l ,F 3 )-l(L'jF"e£ z J 4 
U[C*Z,+X„ +C 1 ^F' i t i + C , (j=) F ' Sx ,] 

c z [ f^e y-_ySF 3 e coszf + £ S’-nZH- J (3 ' 23> 


K + t x z 4 (u~X X 2 = jzc e L z 6 X„ Ur,f 4 /+c t (foxJfeL 10 - 
1 XoL") c.s^ 4 X f [ Yq -C, t X 4 4 C}jU,6L 5c‘m^_2 {zC+L (X i 4£°jcoZ'f' 
+ eL 13 K )]- 

r C ^Z [(£**<) SL'Xosf ] + 

C<z £ djx L 17 sin 2-f ( <$°+ Xj ) +ZCdo JU.L 1 Sin^ J 


(3.24) 
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In the lag equations, terms proportional to (C^ /a)jj have been neglected 
because they are small quantities. The various quantifies used in the last two 
equations are defined below. 


* JL 

2 Mp, 

_jL— . 

ZMfi 









( 3 . 25 ) 


C,rj = - , 

f 5„i - (ftf),,/ 

l 

Ml, 


$D2 ~ 1 ^7 Oc L 4- Cry L \ o 0 

Ml, * 

X H = c,^-c,(zeF'°- f"L) 

Y^- -c„2°+c 7 (L 7 e-zKL e ) 
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The equations above are a coupled nonlinear system of Hill-type equa- 
tions under the influence of periodic forcing. The periodicity of the coef- 
ficients of the equations and the periodic forcing are due to forward flight. 


The terms g and g represent the damping present in the system, 
damping is partly of aerodynamic origin and partly from structural damping 
The damping in the lag degree of freedom is very small; therefore this is 
the potentially unstable mode. 


The 
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SECTION 4 


SOLUTION FOR THE LINEAR SYSTEM IN HOVER 


For the case of hovering flight p = 0, when Eqs. 3.23 and 3.24 are 
linearized, the resulting system of equations is simply: 


** * _ z v * 

X, + $ Dl t-i ~~ y\ h = ° 

X 2 + % z X * + ~ ^ ^ = ° 


The critical condition for the complete system (Eqs. 3.23 and 3.24) 
will be given by the flutter or critical condition of linearized system 
(Eq. 4.1) . The flutter condition is characterized by the existence of a small 
amplitude oscillation for Eqs. 4.1. Assuming the solution in the form 

X 1 = A/** ^ ? X 2 = Az*'^ (4.2) 

Substitution of Eq. 4.2 into Eq. 4.1 yields the following characteristic equa- 
tion 

=° (4 . 3 , 

For a small value of 9, the root of Eq. 4.3 has Real (p) < 0 and the 

solution is stable. For the critical value 0 = 0^, the system is neutrally 

stable. For 9 > 0^, at least one of the roots of Eq. 4.3 has Real (p) > 0 
and the system is unstable. 

At 0 = 0^, there are two solutions to Eq. 4.3 such that p is imaginary 

P~ ~ L(O c 
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Then by setting equal to zero the real and imaginary parts of Eq. 4.3, the 
following relations are obtained 


col = % D> ( ^±l * CUf ' 

+ %D1 


and 


(~ ^ ) (~ " Co c 2 y D1 4 ^c^-H ^ = ° 


(4.5) 


Equations similar to Eq. 4.4 and 4.5 have also been obtained by Ormiston 
(Ref. 3) and Tong (Ref. 7). 


From Eqs. 3.25 it can be seen that g , X , Y are all functions of 0 

D2 H G 

and A , while A is also a function of 0. As shown in Appendix D, the rela- 
o o 

tion between the inflow and collective pitch is given by 



(4.6) 


where 


b e 3 IQ 

C a. 


(4.7) 


Therefore, from Eqs. 4.4 and 4.5 the critical value 0 can be determined. 

c 


46 



SECTION 5 


PERTURBATION SOLUTION IN THE NEIGHBORHOOD OF THE CRITICAL CONDITION 


5.1 The Multiple Time-Scaling Technique 

The multiple time-scaling technique is an extension of the two variable 
expansions introduced by Cole (Ref. 6) , Cole and Kevorkian (Ref. 13) and 
Kevorkian (Ref. 14) . 

In the use of a perturbation method, one may encounter terms of type 
t n sin t, t° cos t, which means that the solution will be unbounded for t -*• 00 . 
These terms (which show the singular nature of the problem) are called secular 
terms: from a physical point of view one does not expect unbounded solutions; 
from a mathematical point of view, it means that the perturbation series break 
down for large t. By introducing an additional time scale, t = £t, where £ 
is some small positive parameter, the original differential equation is changed 
into what is formally a partial differential equation. This can be exploited 
in various ways to examine how the solutions grow and how to suppress the 
secular terms. The multiple time-scaling technique is a further generaliza- 
tion of the two time-scale expansions (Ref. 15) by introducing additional 
time scales t^ = £°t to achieve greater flexibility. In particular, all pos- 
sible secular terms can be avoided and a uniformly valid asymptotic expansion 
in the time domain is obtained. 

The multiple time-scaling technique has been extensively used in treat- 
ing the nonlinear panel flutter problem (Refs. 16 through 18) . It has been 
applied first to the nonlinear helicopter blade-flutter problem by Tong (Refs. 7 
and 8) . 

5.2 The Perturbation Expansion 

Let £ > 0 be some small parameter representing the perturbed state of 
the system relative to the critical condition denoted by subscript "c". The 
following expansion can be established in the vicinity of the critical con- 
dition. Let 


47 



v 2 y z 

X.J- £ a 0 + 6a, v- 6 a 2 + . . . 

< , V 2 / 

X z = £ /> 0 -t-£kyt+ £ h x +.... (5.1) 


be the expansion of the dependent variables. 

The various other parameters of the problem can be expanded as 


c / 4 

y£i. = yro £ -/■ £. +■ - 

e = q c + e 4 t + e 2 s Vz +.. . 

i/ z 

^a e ^oc £ ^2^ *'■ • 


&, = &,„ n„z + 7.,i i + -- 

..o r /^h. 

9 = 9+9 6+ 9 £ 

°4 d 0 dot doz 

X H * X H 0 + X , t + X 2 ^... 


(5.2) 


— 2 _ 2 
“ ^F/O 

_ 2 _ 2 
Cu Z/ ~ (dj LIO 

1/2 

This particular expansion scheme has been chosen so that the 0(e ) per- 

turbed equations will be a system of ordinary differential equations with con- 
stant coefficients. 
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In addition, the multiple time scales will be introduced, using the 
following relations 

'm m = 0, 1, 2, ... (5.3) 


where the \b 1 s are chosen as the new independent variables. Then 
m 


( ) = 2__ 4 ..... 


* * j 2. -s 2 - 

/ ) = _ d + 2.S. — + 

oty' 2 7>% L 


(5.4) 


Some of the perturbed parameters in Eq. 5.2 are not independent, being 
related to 0^ or Q^. The various pertinent relations are given below. 

)„ j ]}-- 

BC 

- Qi (Ta J Qi+boc) _ ( /^. yi<o >oc ) ]■ 

o c l 8 1 ~c^ru ( 5 * 6 ) 

wkei~e 

b 0c = 3ZOc 

CCL 
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§o^J±-(f='ec'F Z ) \ec) 

CvFio 

L-- ^ ( F ' e ‘- F ^‘) 

<ont 

X h0 ~ C H%~ c <( 2 6 c F '°- F" \ oc ) 

= Ch%, -C, (IQ, F'°- F>'\) 

\o~~ - c «3 0 + c i ( L\-z\ oc if) 

( 5 . 7 ) 

Y , : ~ c » 3 0l + c i ( -2^1 L J 

3dzo = % ^cJzcjo L '\ l' 1 \oc & c ) 

Ml, 1 4 ‘ 

= Cj L. [\,Qc -t^oc^ 1 ) 

n i * 


In a similar manner, expressions for X , , g^' y 2 ' ^22' etc "' cou ^ 

also be evaluated. These expressions would be important only if one would 

3/2 

consider expansions of order higher than 0(£ ) and will not be required in 

this report. 

It is important now to associate the mathematical perturbation parameter 
e, which represents basically the ratio between the two time scales used in the 
problem, with physical parameters affecting the rotor blade dynamics. For this 
purpose, a distinction between the two cases must be made: 
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Case 1 : Hovering Flight — p = 0 (i.e., p^ = 0) 

From Eqs. 5.2 through 5.7, it can be seen that 8^ has been undetermined 
and is, in fact, a free parameter. Without loss of generality, it can be taken 
as 0^ = + 1, so that 


e-e*i !5 .a 

i.e., for this case, the perturbation parameter £ is equal to the absolute 

value of the difference between the actual collective pitch setting of the 

blade in the flight condition being investigated, and the critical value of 

the collective pitch 0 . 

c 

Case 2 : Forward Flight — p / 0 (i.e., p 0 / 0, or p, f 0) 

For this case, one may have £ associated withy*.. Depending on the 

value of 0 ) , one can have either: (a) p = 1 and p n = 0 or (b) p = 1 and 
c oil 

Po = 0. Thus 


(a) p 


.1/2 


or 

(b) p = e 

which will be discussed in detail later. 


For the cases with forward flight, Eq. 5.6 is no longer valid. For 
these cases, the inflow is evaluated from the usual relation given in Gesow 
and Myers (Ref. 19) 


Z o = yU. o + 


C T 


while and 6^ are evaluated from 


(5.9) 




I - 



(5.10) 


e-Qc 

t 
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5.3 The Perturbed Equations of Motion 


Substituting Eqs. 5.1 through 5.4 into Eqs. 3.23 and 3.24 and requiring 
terms of the same order in £ to satisfy Eqs. 3.23 and 3.24 separately yields 
the somewhat lengthy system of equations written out in detail in Appendix G, 
Eqs. G.l through G.3. This general system of equations is specialized below 
according to the various cases mentioned in the previous section. 


For hovering flight, p = 0 (i.e., p 

o 

reduce to : 


p^ = 0) , Eqs. G.l through G.3 


1/2 

Equations of 0 (£ ) : 

l— i ( Q'O'bo) ~ 0 

l—l (etc, bo) - 0 

Equations of 0 (e) : 

L i ( 

L z ( bi) ~ ( Q-O, bo) 

3/2 

Equation of 0(e ) : 


(5.11) 


(5.12) 


L 4 ( *i,bz) - Qz ( O-o i boy&ib /) 

L 2 ( bz) - bJ 2 ( do, bo, Clij bi) 


(5.13) 


where the operators L^, and the expressions Q^, N^, Q 2 and are defined 
in Eqs. G.4 through G.9 of Appendix G. 


are : 


For the case of forward flight with p = P q £ 


1/2 

Equations of 0 (e ) : 

L i ( do, bo) - /*■<> to (to, &c) 

L x {ao,bo)= >0 ( to,&') 

Equations of 0(£): 


1/2 


(p^ = 0) , the equations 


(5.14) 


L, ( a «,&,)= A o ^ (to, &c , ciO'bo) ■+ Uo J,( to, &c) + Qi (<*o,bo) 

/ (5.15) 

L Z yLoHl (to, Ki (to } ^c) ~b (do, bo) . 
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3/2 

Equations of 0(£ ) : 

L 1 (& l, b% ) - J^O f~1 ('foflc. J 0*1 5 bl ) f~ JU.0 F 0 (6l) *fb) + yito F^ ( ; io ) 

^ f*o )Q-°) + Qi (cLO,bo ) Q.l,bij 

L z fa Zj />i)^o H i (^ 0/ 9 C) a 4l b 1 ) +JL 0 Ho(% B ')+^o H z ('k&f 0 , 1 *) (5 ' 16) 

+ K z (^ 0j a 0 ) + N z (a 0 ,i>o jd/jbi) 

where the expressions F , H , F , J , h , K, , F, j , H , K are defined in 

0011112222 

Eqs. G.10 through G.18 of Appendix G. 


are : 


For the case of forward flight with p = y^e (p^ = 0) , the equations 


1/2 

Equations of 0(£ ) : 

L | {do, bo) - 0 

1—2. (do, bo) -0 
Equations of 0(e) : 


(5.17) 


L 1 (a 4 , M = JU, F 0 1%j9c) + Qi ( a 0 ,b o ) 

L 2 (a Ho ( % t 9c) ~b ( a °i ^°) 

3/2 

Equations of 0(e ): 

1 fyo d-obo) ~f~Qz (do/bo, a,, bi) 

Lx (d-i'bz) = Ja, Hi { < H 0) 0c )4o,k 0 ) + M 2 (dojbo) bi) 


(5.18) 


(5.19) 


In Eqs. 5.11 through 5.19, the terms F , H , J, , K, represent forcing 

o o 1 1 

functions} F^, F^, H^, represent parametric excitation functions, while 
£> 1 # N^, Q 2 , N 2 represent nonlinear coupling. 


5,4 Solution of the Perturbed Equations and the Conditions 
for Suppressing the Secular Terms 

In this section, the equations derived previously will be treated in 
detail. For convenience, the various cases will be classified in the follow- 
ing manner: 
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Case A — Hovering flight y = 0 (i.e. , y Q = y^ = 0) 

1/2 

Case B — General forward flight case , y = y Q C 

(i.e., y. = 0) and w f 1/2 + £V or u; / 1 + £V 
1 c c 

Case C — Forward flight case, with u) = 1/2 + ev, for this case 

c 

y = y e (i.e. , y =0) 

1 o 

Case D — Forward flight case, with u) = 1 + ev, for this case 

1/2 C 
y = y £ (i.e. , y = 0} 

o l 

In these relations v is a quantity of order one or smaller , which will 
be specified in the analysis. The reason for this classification will become 
apparent within the context of this section. 

5.4.1 Solution for Case A 

This case represents hovering flight (y = 0) . The equations which must 

be solved are Eqs. 5.11 through 5.13. First, the solution to the equation of 
1/2 

0(£ ) will be obtained. As represented by Eq. 5.11, these equations are 

identical to the linear system (Eq. 4.1), treated previously in Section 4. 
Since the damped branches are stable, only the harmonically oscillating 
branches must be considered, which can be represented by 



(5.20) 

where ( )* denotes the complex conjugate; A is a function of ip ,ip , ... , etc., 

o 1-2 

but not of W . In analogy to Section 4, the values of w and 0 are determined 
o c c 

from 





%P40 


(5.21) 
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(5.22) 


$oio $Q2Q (^Fio ~b>uo ) 

+ $020 J 


(% 


%I0 $, 


olo 


<** = CU cX H o Y go 


040 


The quantity v , determined from Eqs. 5.11 is given by 


V, 


.X 


HO L ^>c 


(5.23) 


-<*£ + fa 0 tCOt+ ^F'o 

This quantity represents the ratio between the flap and lag amplitudes. It 
can be shown that 


1^1= J 9 t>2Q Xho 

' Solo Ygo 


(5.24) 


Next, the equations of 0(£) must be solved. The solution of Eqs. 5.12 
is composed of two parts: the solution of the homogeneous system and the par- 
ticular solution which is obtained by the method of undetermined coefficients. 
The homogeneous solution would have importance only if the equations of 0 (£ ) 
were also considered. This being not the case, the particular solution 
would be sufficient. Then 

C 2<*>, t f 0 

, I "Z L P 
'o 


:y z ^( a °IZ u 



(5.25) 


where the quantities a 2 , 0 2> 0 3 are defined in Eqs. G.19 and G.21. 

3/2 

Finally, the 0(e ) equations are treated, substituting Eqs. 5.20 and 

5.25 into Eq. 5.13 and using Eq. 5.26 

* ( 8 t'**)*z teal [AS 


2 Reajt(Az Cu}t ) 2 Reed 

. j <■ (o(-i o)t 

-f (A)*et 


] 


(5.26) 


yield the following relations 
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(5.27) 



Pis (Ao) A 0 J 2 


+ 


+ ( 



} 


L- Z ( a z,bz) = 


2 Real /> 2 Ms +fa A 0 +f> 46 (Ac)*AZ]e. 

^ V/ 



(5.28) 


where the quantities p^ p 2 , p 3 , p 4> p 15 , and p 1& are defined in Eqs. G.21 of 

Appendix G. The terms p ic , p..,. in Eqs. 5.27 and 5.28 represent the nonlinear 

15 lb 

coupling terms . 


5.4.2 Solution for Case B 

1/2 

This case represents the general forward flight case y = y^£ and 

u) / 1/2 + £V or w / 1 + £V. 
c c 

1/2 

The solution to Eq. 5.14 (0(£ ) equations) is given by Eq. 5.20 with 

an additional part representing the particular solution. Thus, 


a ° =2 Real 

bo) 




(5.29) 


where a , b are defined in Appendix G. 

The solution to Eq. 5.15 (0(£) equations) is given by 
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/ 



Z Reat 



<-2 uj t'f'o 

e 




C(l-iVc) % 

e 


~i~ JUi o Ag 



C(1+ioc)h 




+ 



(5.30) 

where the a l3 , b^, a , b.^, a , b , a 3 , 3 3 , a , b 1& are given in Eqs. G.31 
through G.35 of Appendix G. 


to (e 


3/2 


Finally, when Eqs. 5.29 and 5.30 are substituted into Eqs. 5.16 
) equations] lengthy expressions occur which can be written as 


L-i ( ft l>2 ) - Z Real j' £ jo 4 D_ Ao +p^A p + 

P21 A ' 1 Ao + pis (^«)* A c J 4 

( ) t^+ (...) + 


( \ Zwj'k -) 

{ .... J u -f ... Coh*>Tant terms C 


(5.31) 
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L 2 ( 4i,U) = 2 Real j [ p 2 H Ao + Ao + 


+ 


f>46 


(Aoftf] 


T>% 

c.cutSfc / \ Cc^ g f \ 2 3V^ 


+ 


(...) e' 1 % (- -) i 


c(i-^UJc)% 


+ (..■) t Ciu ^ +(...) 


, > CU-SL (oc)% , , 2 

-f- Con St- Terms r 


(5.32) 


where the quantities p^r P 22 are defined in Eqs. G.44 and G.45 of Appendix G. 
5.4.3 Solution for Case C 

This is a forward flight case with to = 1/2 + ev and y / 0 where V is 

c 

a quantity of 0(1) or smaller which will be specified later. From Eq. 5.30 
it can be seen that the expressions for a^, b^ (and their conjugates) be- 
come very large for u) = 1/2 + ev, because M(1 - u) ) of Eq. G.32 is near 

c r^> c 

singular. This situation can be corrected by requiring that ji Q = 0, then 
y = y^e. For this case, the solution is obtained by solving Eqs. 5.17 through 
5.19. Proceeding analogously to the previous cases, the solution to the 
Ofe 1 ^) equation is again given by Eq. 5.20, while the solution to the 0(e) 
equations, Eq. 5.10 is given by 



where the quantities are defined in Eq. G.46. 

The substitution of Eqs. 5.20 and 5.33 into Eq. 5.19 yields 
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L 1 ( a t >i) r Z Rtcd J [fy Df\o ± ^ An 4 F> 1 s (Ao) Ao] e ^ 

, h \* c fy-o )c')H' 0 c( j+^c) % 

-h (Ao) Q 4- (. j £ / 

(...)/ w » + ... } 

(5.34) 


^~i(a. x fa)~ ZRealJ^ fa lAe -(■ fa A 0 + f> /(> (Ao)* Ao^] t 


'-’U'c'Yo 


(/ J ®) fa e 


i{i-02 c )% 


{....) * 


cfy-fcy c ) % 


(...}****• + 


(5.35) 


The quantities p and p defined in Eqs. G.47 and G.48 are associated 
b 6 

with parametric excitation. 

5.4.4 Solution for Case D 

1/2 

This is a case with forward flight y = y e (y, = ° ) and 

o 1 

u> c = 1 + ev. From the solution of the general forward flight case, Eq. 5.29, 
it can be see that the expressions a^ arid will be nearly singular. This 

situation can be corrected by assuming that the expressions associated with 
the forcing F 0 (^ 0 ' ® c ) and 0^) , Eqs. G.10 are mathematically quantities 

of 0(e) which physically means that at the flight condition being investigated 
0^ and g^ are small quantities or in other words flight at low values of C^. 

Then, from Eqs. G.10 


FoWo^Oc )= e ( Foi seni'o - £ F o (io, fa) 


(5.36) 
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(5.37) 


£ (HoiSin+o + Hozcco% ) - £ Ho ('h,6c) 


where F , F , H . , H are defined in Eqs. G.49 through G.52, and it is re- 
ol 02 ol o2 

quired that these quantities be all of order one. 


The equations which must be solved for this case are: 

1/2 

(a) Equations of 0(e ) as given by Eq. 5.11 


(b) Equations of 0(e) as given by Eq. 5.15 

3/2 

(c) Equations of 0(e ) as given by Eq. 5.16, except that in this 

case, the term y F (ib , 8 ) replaces y F (0, , \b ) and y H (ib , © ) replaces 
oooc ool o oooc 

y h (ij) , 9. ) 

o o o 1 

The solution to this system of equations is given by Eq. 5.20 and 



(5.38) 


where the quantities a,,, b , a, ., b. , a,_, b n , a. c , b. are equal to the 

1J 1 J 14 14 lb lb lb lb 

same quantities without the bar, when a „ = b = 0. For convenience, the 
_ _ _ o2 _ 02 

quantities p^ through p^ and p through p are given by Eqs. Q.53 through 
G. 60. 
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The substitution of Eqs. 5.20 and 5.38 into the modified equations 
(Eq. 5.16) yields: 


"fMc? ± b + /*2 / A* f 
3 


L I ( aa,i>z)= ZRwt {[ pi 

p^(f\o)*f)o Je'^ 0 + [/**/>$ 

c (Z- OOc) % 


+/?(*•)* I>„ e 


-h 


n , , * duic'H, 

hoMih^e +(-••) e * 


f-J e 


C (&<*>*) <f e 


* (-...) e 




-A • * ' - 


k Mv +p^/jo-f- p^j.o z A 0 -h 

Pl 6 (fio) $0 J 6 -f- £ Jj.t> p)0 4-^p 3 piy + 


Opo 


A 0 fab] ^ ^ (fa )* f> 2g e 






t fa-^c) % 
Ciu>t % 


+ (...)€ 


:|a+w t )f 0 \ L(i + Zioc) % 

^ (...) e. + 


C\P- 
<2 ° 


(5.39) 


(5.40) 

where p , p 22 are equal to p 21< p 22 when a Q2 = b Q2 = 0, and the quantities p 23 
through p^ 0 are given in Eqs. G.61 through G.68, 
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5.4.5 Conditions for Suppressing the Secular Terms 
for the Various Cases Considered 

In this section, a general condition for suppressing secular terms of 
type ip e 1Wc ^o will be derived. As shown below, the condition for suppressing 
secular terms can be obtained from a relatively simple consideration. 


A general case in which secular terms would appear , can be formally 
written as 


L t ( 0-i, - Iff a 




L 2. ( Ul, bx) = H 2 ^ 


i.'COc'h 


The solution to Eqs. 5.41 can be written as 

— i— clot 

OUz = II 40 e 

, __ icOc^i 

»zo C 


(5.41) 


(5.42) 


From Eqs. 5.41 and 5.42 


2 T- 

-c -f Cj Dio C (U C + U)p l0 


- uC ^c 


Y(xo c 


CrO v CaJc 


Ujc + C( *>c°j BZo 4 ’ (jJ uo 


-1 

s’ "■ 


C ' 


// 40 


If* 


H -Z0 

- ' *< 

> 

r i 

J 



_ 


It has been previously shown in Section 4 , that the matrix of the coef- 
ficients in the last equation vanishes. Thus, these equations are not inde- 
pendent of each other and the first row of the matrix must be proportional to 
its second row, denoting the constant of proportionality by C^: 


c 



z 2 

- °bc +$dio<'Wc+Wfio 

- Y(tO 



c LOc 


- U) c L jCct) c %j,i 0 -l- (xi L ,p 


(5.43) 


From Eqs. 5.43 and 5.23 



y; 




<ko 




X. 


HO 


- o 


(5.44) 
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As shown below, taking the appropriate combination for tt^ and tt^, 

Eq- 5.44 will give the required condition for suppressing the secular terms 
for the various cases considered. 


Cases A and B 

Case B which is the general forward flight case with u) / 1 + VE or 

1/2 c 

ui / 1/2 + ve and y = y £ , will also include Case A as a particular case, 

co 

From Eqs . 5.31 and 5.32 


Hi = pi + p 3 Ao + / Ao 2 y D 1/ Ao + PlS fo 0 ) ^ o 

(5.45) 

"it 2 , 

= A + puftof P -0 P zl Ao+ p/6 (A-o) Ao 

W, 


From Eqs. 5.44 and 5.45, the condition for suppressing the secular term 
can be written as 


24s. =(K. Z + L? 

D+, 7 


(5.46) 


For the case of hovering flight. Case A, y = 0, (y ■= 0 

o 1 


) and 


Eq. 5.46 reduces to 


Ma.= A : Z A 0 + H, RoiAo)* 

1>% 3 


(5.47) 


Evidently, Eq. 5.47 is not subject to the flutter frequency limitation 
of Eq. 5.46. 

The quantities K , K , K , k are defined in Eqs. G. 69 through G.72. 

1 2 3 o 

Case C 

For this case, u) = 1/2 + ev and y = y,e. Where V is of order one or 
c 1 

smaller, it can be shown that 
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(5.48) 


For this case 

o y* 

Tf<= P* + P 3 + P< s' +J**Ps(. bo)* e 


11 z" 


yD 2 Mfi. + p H bo+ pH, Ro Mo) 4- J^ 4 fo (fto)*e 


_CW, 


(5.49) 


and from Eqs . 



5.44: and 5.49, the condition for suppressing the secular term is 

Ao + K 3 ho ( Pof +^-4 e, 1 (5.50) 


! where is given by Eq. G.73. 

Equation 5.50 can be rewritten in a more convenient manner after chang- 
ing the dependent variable. Let 


r * c0 ^ 

Ao-fi<> e (5.51) 


Then Eqs. 5.50 and 5.51 yield: 

14a = k x /■ K A> ( &°) +}*■< (A?) 

■at, 7 


where 


K z = AC 2 f<-V 


(5.52) 


(5.53) 


It should be mentioned that a case to = 2 + £V could also have sane 

c 

practical importance, when the harmonic components of the inflow are considered 
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together with cyclic pitch variations. In the present report, these effects 
will not be considered. A partial treatment of this case can be found in 
Ref. 8. 


Case D 


For this case, the flutter frequency is close to the forcing frequency 


U) = 


1 + ev and y = y £ 
o 


1/2 


It can be shown that 


^ i / (%! ) %'0 > % 1>10 > toFM, <"uo) 


(5.54) 


From Eqs. 5.39 and 5.40, it is clear that many terms will contribute 

ilk-, i(2-UL.)U> -i(l-2u) ) ill 

secular terms. In addition to e <J , terms of type e c ° and e c r o 

will also yield secular terms. For this case 


///-= />, Mfl + f> A 0 + p z /c.„M 0 +P- Ro (hof + 


r / p 

L Z* 0 A 0 />2Sj e 

, z , . . z , * COp 

+ h* P27 e + / h° (hZ) ^ 


^ + j + Z A a + p 1b (Poj A 0 -t 

-t-£/*of>J0 +^o p ZH + Jj.of\o Mo)*/>26 J e - 


z , < %■* z 4 C\)p, 

+J +0 (Ao) p z7 Z 4-Ju.oPo (pty) £ 


(5.55) 

From Eqs. 5.44 and 5.55, the condition for suppressing the secular term 
can be written as 
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Ih. 


K 6 ) A 0 + K>i (\ 0 (Ao) -f- 

[ fj'Q Ky 4-j*.o + juo A 0 (Ao) Kg J 



+ j^otfioY K< 


e 




•/o 



(5.56) 


where the quantities K^, K^, K^, tc^ , and/t^ are given by Eqs. G. 74 through 

G.78, while < is equivalent to < with a _ = = 0. 

6 6 o2 02 

Equation 5.56 can be rewritten in a more convenient manner by changing 

the dependent variable. 


Let 



(5.57) 


Then Eqs. 5.56 and 5.57 yield 


HAo = ( KSfZ- K( ) Ao + *1 V (Ao)*-f 

h % 

ky ^7 + J* 0 rtf ( Ao) Ao] 


+ 


where 


~f~ JU o (Ao) k S ^ J^-o Ao Kf 0 


(5.58) 


/Cj.» K X +CJ 


(5.59) 
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5.5 Stability Analysis and Blade Response Amplitudes 
for the Various Cases Considered 

In this section, the amplitudes of blade response will be determined, 
and the stability of the motion will be investigated for the various cases 
considered in the preceding sections. The stability will be investigated 
using the conditions for suppressing the secular terms, while the amplitudes 
of blade response will be determined by taking appropriate combinations of 
the solutions to the perturbed equations of motion. 

5.5.1 Case A, Hovering Flight 

For this case, y = 0 and the condition for suppressing the secular term 
is given by Eq. 5.47. Equation 5.47 can be solved in closed form, taking the 
solution in the form of (Refs. 17 and 7) 

Ao = f e ^ (5.60) 


where both p and <j> are real. 

Substitution of Eq. 5.60 into Eq. 5.47 yields 

— £ - Kt. r p+ 

Irjr ' + 


(5.61) 


(5.62) 


where the subscripts R and I denote, respectively, the real and imaginary part 
of the appropriate quantity. 

The solution of Eq. 5.61 can be written as: 

'-<)y (5.63) 

From Eqs. 5.62 and 5.63, the solution for <{> can be written as 
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to 


+ K.Z I % - ^31 


ZK-. 


m 


^ f J ^ 3 * ( j _ 

1 ° Kza 


2 Kir fa 


(5.64) 


where the quantities p^, <p° are the values of p and <J> at = 0. 

For the case of hovering flight, the solution is given by Eqs. 5.20, 
5.25, 5.1, 5.63, and 5.64. 



- 2 RtaJL 



l'COl'J'o 

e 





£ 


-h 


f £Ao(A.)* 



(5.65) 


For the case of hovering, from Eq. 5.8, the perturbation parameter e is 

e = 1 0 — 0 I , so that the blade response depends only upon the collective 
c 

pitch setting. 

The stability of the blade response can be investigated using Eq. 5.63. 
The following four cases can occur. 

Case 1: K_ >0 and > 0 

2R 3R 

At lit *= 0, p = p , at some finite time, the denominator of Eq. 5.63 
1 o 

approaches zero. That is p -*■ 00 . This clearly represents an unstable case. 

2 

Case 2: If both < 0 and < 0, p 4 0 as f, ^ This clearly 

■ " JR X 

represents a stable situation. 

Case 3: If > 0 and K„„ < 0, the exponential terms in Eq. 5.69 tend 

2R jR 

to zero and 
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(5.66) 


f. 





This case represents a stable limit cycle oscillation. 


Case 4: If < 0 and > 0, 

_____ - — jR 


Denoting 


A.*. = f ~ ’ 

1 k 3K 


Eq. 5.63 can be rewritten as 


J_ - L . 4 


-Zt„% . 

e ■ 1 


i 2 U‘ 


2. 

u 


-h 


Su 


from the last expression for 

1 - JL -CO 
S/ fee 

The system is unstable and for 


_L.~_L.yo 

f. z fu 


the system is stable. 


From Eqs. 5.66 and 5.68 for 


s: > - 


K. ZR 

K 


IR 


the system is unstable, and from Eqs. 5.66 and 5.69 for 


f X 4 - Kz* 
*'0 

^3# 


the system is stable. 


(5.67) 


(5.68) 


(5.69) 


Thus, in this case, the stability is conditionally dependent upon the 
relative magnitude of the initial value of p when compared to p. 

O X/ • C • 

From Eqs. G.69, G.70, G.24, G.25, and 5.7, it is clear that K is 

'V 

linearly proportional to 9. . That is, the origin (A 't 0) is unstable above 

1 o 

the critical condition (9^ > 0) and stable below the critical condition (9^ < 0) . 
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5.5.2 Case B — The General Forward Flight Case 

1/2 

In the general case with forward flight, y = y Q £ (y^ = Q) without 

any loss of generality y can be taken as y =1. Thus, in this case, the 

° 2 ° 

perturbation parameter is E = y . The flutter frequency for this case can have 

any value except u> ^ 1/2 + ev or w / 1 + ev. 

c c 

The condition for suppressing the secular terms is given by Eq. 5.46. 

The solution of this equation is identical to the solution of Eq. 5.47 with 

2 2 
6 + y 6. replacing 6_. In analogy to Eq. 5.63, p can be written as 

z o 6 2 



[fclli+jAo fyn ) ^kf 





( Kx* +/- o 



(5.70) 


The various cases for which the stability of the response was investigated 
in the previous section remain unchanged and are summarized below: 


Case 


K 2R 

+ 

2 . 
KT y : 
6R o 

Case 

_2: 

K 2R 

+ 

v h R • 





2 

Case 

_3: 

K 2R 

+ 

M o K 6R 


3R 

*3R < 
K 3R 


with 


f ec = -\l +/*'?*■** 


(5.71) 


Case 4 : + y Q Kg R < 0 and <^ R > 0 — this system is conditionally 


stable 


Thus for 


fo * fu. 

and for f* > 


the system is stable 
the system is unstable 
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From these relations it is clear that the forward flight can influence 

the stability through the magnitude and sign of In addition to this 

bR 

effect, forward flight can also affect the stability of the system by influenc- 
ing K . 

* 2R 

From Eqs. G.22, G.23, G.24, and G.70, it is clear that K is dependent 
only upon the critical conditions and (with 0^ = 0) is composed of two parts, 
one proportional to 9^ and one proportional to A . 

^2 .= (*164)61 + (k 2 \,)\i 


where 

as 


K 201 


and K 


2Al 


are given in Eqs. G.80 and G.81. 


Thus K can be written 
ZR 


“ K 


(5.73) 


In forward flight 0.^ must be evaluated from Eq. 5.10. The flight condi- 
tion is determined by a fixed value of C T and a known value of p. For a given 

flight condition, the values of 0 and a are evaluated using an approximate 

R 

trim calculation described in Appendix F. These determine 0^ by Eq. 5.10. 

Thus, a change in (J , at a fixed value of C T> results in a change of 0^ and A^, 
which has a considerable influence on both the stability and response of the blade. 

From Eqs. 5.1, 5.29, and 5.30, the blade response for this case can 
be written as 



J ' [bit J 
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As can be seen from this equation for the lower order terms (multiplied 

by y) , only the nonlinearities and the forcing have an effect, while in the 

2 

higher order terms (multiplied by y ) , the nonlinearities, the parametric ex- 
citation, and the forcing will all have an influence on the blade response. 

5.5.3 Case C — Forward Flight with Dominant Parametric Excitation 

For this case, the flutter frequency = 1/2 + £V. As shown before 

in this case, y = y,£ (y = o) and y, can be taken as y, = 1. Therefore, 

1 o ' 1 1 

the perturbation parameter £ = y. The condition for suppressing the secular 
term for this case is given by Eq. 5.50 or Eq. 5.52 which cannot be solved in 
closed form. Conclusions regarding the behavior of A^ can be reached without 
actually solving Eq. 5.52 by using a phase plane analysis as described in 
Refs. 21 and' 22. This method was applied to a similar problem by Tong (Ref. 7) 
and Kevorkian (Ref. 14) . Suppose A q can be written as 


A.= i+cr[ ; f“'7 

IAJ 1 

From the last relations and Eq. 5.52 

- (k x r fciV? ■'■P 3 ^3* ) J - ( Kxx-U, fCs-i J-f K’3i)' r l~ 

^2 = (K iZ +JU< Ks-l-h f’XizH*' ( *** Ks* + f'k**) *1 ) 


A singular point occurs when both expressions in Eq. 5.76 vanish 
simultaneously; thus a singular point occurs at £ = r| = 0, then p = 0. 

According to Minorsky (Ref. 21), the investigation of the behavior of 
the equations near the origin can be limited to the linear system which can 
be written as 


~~ - ( ) f - ( kn 
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- ( Kh+Lli Ksi) | + -/,***) 

5 'ij ' 


(5.77) 


A proof which justifies the neglection of the higher order terms in 
Eq. 5.76 can be found in Minorsky (Ref. 21). As shown in Ref. 21, the 
stability of the system in the vicinity of the singular (or equilibrium point) 
can be qualitatively determined by using the characteristic equation 

Km fyUj - /Cxi +Ju 4 K, r 

Kxr + *s-r Km -ju 4 - S 

The solution of this quadratic equation for S yields 

^,Z S + Kzn ~ 'Jkir'- l^ij +-/*i iKsj 


(5.78) 


which is also equivalent to 

S 4)1 = Km - ]J juflKs! - Kzl 

According to Ref. 21, for distinct roots of Eq. 
cases can occur: 


(5.79) 

5.79, the following 


(I) 


S^, S 2 real of the same sign. The singular point is a node 
which is stable if S^, S 2 are negative and unstable if they 
are positive. 

The roots are real if y > < 21 ^ I K 5 1 • They will have the 
same sign if 



KzL 

Ik, I 1 
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or 


fX* ^IK.-i'Kxx = 1^2./ 

IK S I Z /kfl 2 - 

From the last two relations, the singular point is a node when 

— , z 


- 2 

Jb£. 

lKsl 1 


+ L Z L. l!£±! 

^ l*st 2 


(5.80) 


that 


The node is stable when < 2R < 0, and is unstable if < 2R > 0. Recall 

K 2 = 


Therefore the node is stable when < 2R ^ 0 and is unstable when < 2R ^ 

(II) S , S are of opposite sign. The singular point is a saddle 
1 2 

point and the equilibrium is always unstable. This will 
occur when 


K'ZR’*' \j /*** - 1^21 I 
Y I A*/ 7 - 


or when 


A 1 "> 


IfCzl 1 

l« 5 l Z 


(5.81) 


(III) If S^, S 2 are complex conjugates, the singular point is a 
focus. For this case, the stability of the singular point 
is determined by Re(s). If Re(s) < 0, the focus is stable; 
if Re(s) >0, it is unstable. 

The roots S^, S 2 are complex conjugates if 


A 


2 ^ j^L 

IK;I Z 


(5.82) 


For this case, Eq. 5.79 can be rewritten as 


5^,2.= K-zr - c ^$7 r ji C zZ J ~ 77 

V /KW 1 ^ 


(5.83) 
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The stability depends upon the real part of S : 


If (C < 0, the focus is stable 

If K > 0, the focus is unstable 

2R 


It has been shown by Tong (Ref. 7) that additional singular points be- 
sides the origin can exist. These will have a relatively minor influence on 
the considerations which will be given below; therefore, these additional 
singular points will not be written out in detail. 

In the analysis given above, the behavior of Eq. 5.52 has been analyzed 
in the vicinity of the singular (or equilibrium point) . Another item of major 
importance is the determination of the conditions for the existence of closed 
trajectory curves (or limit cycles) of Eq. 5.52 in the £,f|-plane. For this 
purpose, the Poincare-Bendixson theory can be used (Ref. 21) . In order to 
apply this theory the following quantity must be defined. 


V = 


'IX + 'SL 

dj ^ 


(5.84) 


From Eqs. 5.76 

— ~ + 3, I" hc 2 n - -if ^ ^jr 


(5.85) 


and 

v = Z (Kzr + 2- 


From Green's theorem 

J V ^(Xdrj-Yd^ 


(5.86) 


(5.87) 
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According to the Negative Criterion of Bendixson, if V does not change 
sign (or vanish identically) with a region D of the £,r|-plane, no closed 
trajectory can exist in D. 


Therefore, the following possibilities exist: 


- 2 

(a) If > 0 (they have the same sign, p > 0) 

ZR jR 

This means that V cannot change sign or vanish and the 
integral (Eq. 5.87) cannot be zero. Therefore no solu- 
tion of Eq. 5.52 can form a closed trajectory in the 
£ ,r)-plane. 

(b) If K K < 0 (i.e., they have different signs) 

2R JR 


From Eq . 5.86 


= ^ ^ 3 * f 

2-K-bA. 


+s 


') 


(5.88) 


for p < — (K ) / (2K ) , V will not vanish and Eq. 5.52 can 

2R JR 

have no solution which forms a closed trajectory in this 

region. Therefore, if there exists (the negative criterion 

of Bendixson does not quarantee existence) a solution of 

Eq. 5.52 which does form a closed trajectory in the £,r|--plane, 

. . 2 

its minimum radial distance from the origin, p . must 


mm 


satisfy 


f . 


- tj-ZK 


(5.89) 


The upper bound of p . can be estimated when rewritting 

min 

Eq. 5.52 in polar coordinates. 

_ ^ 

Let 

and K s = /fC r / 

from these and Eq. 5.52 
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(5.90) 


+/L 1 lK s lcos(lS-(l>)+K iR f t ’]f 

7 

~[Kxl -/ulltsl^n(2f-t)+ K 3 Zf l J 


The upper bound will be estimated without making use of Bendixson's 
negative criterion. The quantity p is always positive. Therefore, the right- 
hand side of Eq. 5.90 will be monotonic and will not change sign outside of 
the range 


me tx (o t 


K*R l*±! 




(5.91) 


__ f ( K-sl 

Km ’ iKs/tl 1 J fCm / '~ 1 jK. : 3K ( 

From this relation, if Eq. 5.90 has a solution in the form of a closed tra- 
jectory in the £,r)-plane, its radial distance from the origin must be bounded 
by Eq. 5.91. 


The last bounds, Eqs. 5.91 and 5.89, can be combined to give a bound 
on the limit cycle amplitude, if it exists. 


max 


(- -u,LM -ox-i-h 

2 *3* r /KmI' 


Ik si 
I Kir. I 


(5.92) 


Note that in both cases considered above, p is a monotonic function of 




For large amplitudes (i.e., large values of p) , the behavior of p can 
be approximated by 

3 


1L 

1% 31 


(5.93) 


From Eq. 5.93, if K < 0, p will monotonically decrease in the region 

3R 

of large amplitudes; i.e., the solution will converge to a limit cycle if it 
exists , or to a stable equilibrium point (or singular point) . Therefore , this 
represents a stable situation. 


If < > 0, will monotonically increase in the region of large ampli- 

3R 

tudes; i.e., the blade is always unstable, if the disturbance is large enough. 
From Eq. G.71, it is clear that depends mainly upon the nonlinear coupling 
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of the system, so that the nonlinear coupling is the decisive factor at large 
disturbances. From Eq. G.73, it can be seen that depends upon both the 
parametric excitation and nonlinear coupling with response due to forward 
flight. Therefore, on the basis of this approximate analysis, it can be con- 
cluded that the location of the singular points and the magnitude of the 

limit cycle oscillation (if it exists) will depend upon the parametric 
excitation and nonlinear coupling. 

In addition to the effects discussed, a change in the forward flight 

velocity with C fixed will affect K by changing 8 , and A according to Eq.5.72 
T 2R 1 J. 

discussed in the previous section. Therefore forward flight will considerably 
influence the stability of the system and it will also affect all of the bounds 
obtained for the limit cycle amplitude in this section. 


The amplitudes of blade response can be obtained from Eqs. 5.1, 
and 5.33. 



5.20, 


(5.94) 


Due to the convenient form of the solution (Eq. 5.94) , it is possible to 
rewrite Eq. 5.52 in terms of the actual physical quantities associated with the 
problem in a manner described in detail in Subsection 5.5.4. 


Let 


'h ~ 

% t --£ K-- 






Then Eq. 5.52 can be rewritten as 

K s (X-t) + l 


(5.52a) 
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When considering only the linear part of this equation which is equiva- 
lent to Eq. 5.77, the last term in Eq. 5.52a can be neglected. Consequently, 
the characteristic equation (Eq. 5.79) can also be rewritten in terms of the 
physical quantities as 


(fo & c)+ (K'i\,)0'o-h°c)- \j fo lk s l - ^(6 -6 C ) (c*) c - ) J 


(5.79a) 


The solution of the linearized system is the same as the one derived in 
Appendix H for Eqs. H.2 and H.3, and can be written as 

7 n o^o 

A- 2 - D 4 £ -b p 2 e. 

where D and D are equivalent to D , D used in Appendix H when K = 0 and 
2 1 2 12 6 

y K g is replaced by y< 5 in the appropriate relations. As is shown in 
Appendix H, A for the linearized system can be written as 


A a = *2 & 




A' 


>a 
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l + ev 


5.5.4 Case D, Forward Flight with a) = 

c 

This is a case in which the flutter frequency is close to the forcing 
frequency. Since the system is being excited at its resonance frequency, the 
amplitudes of blade response can be expected to occur at moderate levels only 
when the excitation is assumed to be small (see Section 5.4.4). The excita- 
tion is represented by pF (ijj ,0 ) ;pH (t l) ,0 ) where F , H are given by 

oocooc oo 

Eqs. 5.36 and 5.37. In most cases considered, cannot be assumed to be too 
small (see Section 7). Therefore, it will be necessary to assume that the 
solution derived in this section is valid for small values of the advance 
ratio p and small values of C T . 

For this case, a complete phase-plane analysis similar to the one per- 
formed in the previous section is possible, but cumbersome; therefore a some- 
what different approach will be used. Using Eqs. 5.1, 5.20 and 5.38, the solu- 
tion for this case can be written as 



(5.94) 
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Due to the convenient form of the solution for this case, it is possible 
to rewrite Eq. 5.58 in terms of the actual physical quantities associated with 
the problem. For this purpose, it is convenient to define a quantity 



(5.95) 


where A q is given by Eq. 5.57. Thus, the first-order solution for the lag 
motion, according to Eq. 5.94 is given by 


% 2 = 2 Reat ( X 2 ) = Zt! ,x fieae(Ao) 


(5.96) 


Using Eqs. 5.59, 5.72, and 5.73 together with 
bJc~ ■/■/■£]? 

the following relation can be written 


^2 = faze, ) ( 0 ~ 6c ) (K 2 ^ ( A o- \o c) + L CCo r -l) 


(5.97) 


- 1 / 2 . 


Using Eqs. 5.95 and 5.97, together with y = y Q £ and ^ Eq. 5.58 

can be rewritten in terms of the physical quantities 

~~ = [ ~^ot\+ t '(u> c -l) 1^2 ] 

OL Tc 


+ (%z)* + ju, [iCtlXz)** K io x z] X 2 + 

yU + /U? Kq 


(5.98) 

2 — * 

In the last relation, the term K x? (X) I s due to the nonlinear coupling 
of the flutter mode, while the term y[K Q (X->) + K iri X-J X-> I s due to the non- 
linear coupling of the flutter mode and the parametric excitation. The solu- 
tion of Eq. 5.98 in closed form is required in order to investigate the non- 
linear response of the blade. Due to the complex form, Eq. 5.98, this would 
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be very complicated; therefore, the linearized form of Eq. 5.98 as given by 
Eq. 5.99 will be treated first.' 

i£k r f K Z& J0-9 C ) + (VAoc)*/^^^ 2 

L 

+ > Z /C f ^2) 7 

(5.99) 

Let — 

X, ^ f -/-to y 

3 1 (5.100) 

then the real and imaginary parts °£ Eg. 5.99 yield a system of two first- 
order differential equations: 

41 * I (n aa , \le-Bc)* (Ki\X 

♦(wc - v ) +/t l t< r + yU z [^3i ^vl] 

-f ^ r 7K 

(5.101) 

4a » l(K !e ,) K (e-«c */ku] 7 + 

d% l 

[(, k u ,jj (e-ty+fUi * + 

y(x,^- «« n\*/* k Kn + y Kni 

(5.102) 
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The solution to Eqs. 5.101 and 5.102 can be found in Ref. 22, and it 
is analogous to the treatment given to Eqs. 5:77 in Subsection 5.5.3. For the 
sake of completeness, the solution of these equations is given in Appendix H. 

The solution of the homogeneous linearized system is given by Eqs. H.6 
and its stability depends upon the roots of the characteristic equation: 

5 «,i = ) + ( K 'z\ l ) 1 ^ 0 '^ oc ) %z] 

(5.103) 

The linearized solution is stable if Re(S^) and RefS^) < 0, and the stability 
boundary is given by 

Re(S,) = 0 (5.104) 

The solution to the complete nonhomogeneous equations is given by 
Eq. H.14. From this relation, an approximate relation for A^ can be written as 


A = ' 

n Q - — 




h e 




* c %] 


(5.105) 


It is important to note that due to the nonhomogeneous form of Eqs. 5.101 and 
5.102, the flutter mode for this case will always be excited. This is also 
evident from Eq. 5.105. 

As mentioned in Subsection 5.5.2, for trimmed forward flight with a 

fixed value of C T , the variation of y will strongly affect the quantity 

(0 - 0 ) and (X - X ). Thus, the stability boundary represented by Eq. 5.104 
c o oc 

will be strongly dependent upon the advance ratio. 
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SECTION 6 


UNIFORMLY VALID SOLUTIONS 


6.1 Uniformly Valid Expansion 

From the numerical results obtained, which will be discussed in detail 

in the following section, it is clear that the cases of practical interest will 

occur at values of collective pitch 0 < 0 . Therefore, it is important to ex- 

c 

tend the validity of the solutions obtained in the previous sections (which 
are valid near the critical condition) to a region which extends below the 
critical condition. 


This problem is discussed in detail in Ref. 8 where a conventional 
method (Ref. 6) for matching asymptotic solutions near and below the critical 
condition is used. 

% 

The matching is achieved by using an additional small parameter £ which 
has the following properties: 


1 » £ ><> and JL << i 

e 

e is related to the parameters of the problem by Eq. 6.1: 

Q-6c = - £ 

A/ 

£ 


( 6 . 1 ) 


The last equation can be also rewritten as 


£ = - 


e-e c 


_ ± 

Qi 


( 6 . 2 ) 


In the cases with forward flight, the perturbation parameter e can be 


associated either with £ = y (Cases B and D) or with £ = y (Case C) . 


From the numerical results presented in Section 7, it is clear that in 
trimmed forward flight at fixed C T , the range of variation for | 0 — 0^ | is 

0 < I &-0cl ± • IZ 

while the range of variation for y can be taken as 
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Then for Cases B and D 


£ = 




(6.3) 


and the requirement that E < 1 may not be satisfied for certain combinations 
of C T and li. 


For Case C 



(6.4) 


Therefore, for this case, the uniformly valid expansion will be correct for 
y < 0.15, approximately. 

If the expansion near the critical condition and below the critical 
condition are denoted by subscripts "n" and "b" , respectively, the requirement 
for matching the solutions is (Ref. 8) : 


X, 

Xz 




a 


<x 


o 


(6.5) 


'V 

as e, E -*-0 with fixed, for a = 1/2, 1, 3/2 and the uniformly valid ex- 
pansion can be written as 




( 6 . 6 ) 


where the last term in Eq. 6.6 is the common part in the expansions 



6.2 Uniformly Valid Solutions for the Various Cases Considered 


The uniformly valid solutions for the four cases considered in 
Section 5 can be easily obtained. 

Case A, Hovering Flight (M = 0) 

For this case, the solution below the critical condition is 


M (*. 

% Z J~° while jx ) •‘• s < 3^ ven by Eqs. 5.65. 


According to Eq. 5.72, is linearly proportional to 0^, below the 
critical condition 0^ = 0 - 0^ and is a negative quantity. For this case, 
the last term in Eq. 5.47 can be neglected and its solution can be written as 

A -A- 

n 0 ~ At e (6.7) 


where A. is a constant determined from the initial conditions. The solution 
i 

will be exponentially decaying below 0^ and Eq. 6.5 is satisfied. 

Case B, General Forward Flight Case 

1/2 

As mentioned before, for this case h = e and the flutter frequency 

can be arbitrary, except to = 1/2 + ev or u> = 1 + ev. According to Subsection 

c c 

5.5.2 for 0 < 0 c and p < 0.3, Eq. 5.46 represents an exponentially decaying 
function. For this case A q can be approximately taken as 


A,- Ac « 


k z 


( 6 . 8 ) 


The solution for / ^ £ is given by Eq. 5.74. The solution below the critical 

I 

condition, with Eq. 6.8, will be composed only of the parts of Eq. 5.74 which 
are independent of A q . These parts will match and Eq. 6.5 will be satisfied. 
From Eqs. 6.6 and 5.74, the uniformly valid asymptotic expansion can be 
written as 



= Z Hod I UA J v j y* 


a °* l t*‘ 


+ lx ho f °^7 


CllOl'f-', 


•f A z ^ 7 Q ' 3 ] e 

' ( ^43 J 




+jdAo j a <H 1 l 


a «u 




4 L z {L)(Ao f It 3 }*/ 1 
/y I 


(6.9) 

Equation 6.9 is similar to Eq. 5.74. The only difference is that for 
Eq. 6.9, the various quantities a q2 , b Q2 , a 2 ' &2 r a 13 X Ho Y Go' ’**' etc * ' 

must be evaluated at the actual value of 0 and \ as determined by the trim 

o 

calculation, instead of evaluating them at the critical condition. 

Case C, Forward Flight with Dominant Parametric Excitation 

For this case, the flutter frequency is = 1/2 + ev and p = e. 

Below the critical condition A q can be taken as given by the last equation 
on page 46. The flutter mode is not excited and will be a decaying oscillation. 
Thus, the solution below the critical condition will be independent of A q and 


will match the appropriate part of 


tel 


, satisfying Eq. 6.5. From Eqs. 5.94 


and 6.6, the uniformly valid expansion can be written as 
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( 6 . 10 ) 

As pointed out in Subsection 6.1, Eq. 6.10 is valid for 0 < y < 0.15. 

In Eq. 6.10, the various quantities , 3' 2 » C t^> ••• » etc. must 

be evaluated at the values of 6 and X as determined from the trim calculation. 

o 

Case D, Forward Flight with 0 ) «* 1 + ev * 

c 

1/2 

For this case y = e y and u =1. Below the critical condition A 

o o o 

can be approximately determined from Eqs. H.14 and H.15 given in Appendix H 
and the solution for is given by Eq. 5.94. As pointed out previously 

(*ifl 

in Subsection 5.5.4, for this case, due to the nonhomogeneous form of Eqs. 

5.101 and 5.102, the flutter mode will always be excited; therefore, the 
matching procedure is no longer simple. The solution JH l can be 

L *iJh. 

obtained only by solving Eq. 5.98 which is difficult to obtain in closed 

form; therefore / J cannot be evaluated. It can be concluded, 

%zJ common 

therefore, that a uniformly valid expansion for this case is difficult to. 
obtain analytically. Therefore, this case can be handled most conveniently 
by direct numerical integration. As mentioned in Subsection 5.5.4, the analyti 
cal treatment for this case is correct only for small values of the harmonic 
forcing which implies low values of y(i.e., 0 < y < 0.20 ) Therefore, direct 
numerical treatment of this case will have the additional advantage of not 
being limited to small values of y and C T * 


88 



SECTION 7 


RESULTS AND DISCUSSION 

7.1 Numerical Quantities Used in the Calculation 

In this section, the numerical values of the various quantities used in 
the computation of the results are given: 

(1) Mass distribution was taken as a constant along 
the blade span 

(2) c = 0.010 in all of the calculations 

d 

o 

(3) Lift-curve slope, a = 2TT 

2 

(4) Equivalent flat-plate area was taken as f/irR = 1/100 
unless otherwise stated. This value of f results in 

C = 0.01. This quantity is typical for modern well- 
designed helicopters 

(5) The mode shape in flap is the same as the mode shape in lag; 
both were approximated by the first nonrotating mode shape 
which can be approximately written as 

Vti) 

This relation was taken from Ref. 12. 

Equation 7.1 satisfies all of the boundary conditions of the 
problem. As pointed out by Bramwell (Ref. 20) , to get the correct 
first rotating flapping mode shape of a hingeless blade, as many 
as six or more expressions of type Eq. 7.2 should be combined: 

] (7 . 2 , 

Approximation (Eq. 7.1) to the mode would have influenced the re- 
sults to a certain degree if the rotating flap and lag frequencies 

were calculated by using the (El) , (El) , and m distributions. 

Y z 
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were 


In fact, in all of the calculations, the values of u) , , u) , 

FI LI 

selected so as to give a certain w , or 6 and the mode shape 

c c 

was kept constant during all of the calculations. 

(6) The flap and lag coefficients, F 1 , L 1 defined in Appendix C (and 
mass quantities defined in Appendix B) were numerically evaluated 
using a seven-point Gaussian integration. This integration scheme 
is accurate up to polynomials of degree 13 (Ref. 25) . The weights 
for Gaussian integration were taken from Ref. 26. The numerical 
values for these coefficients for SL/R =1.0 and A = 0, B = 1 are 
given in Table 1. 

*\# 'V 

(7) The quantities C , D^, F^ , G^ , ri^(l) used in the trim calculation 
described in Appendix F were evaluated using the approximate equa- 
tions given for them in Ref. 20. The approximations of these rela- 
tions in Bramwell (Ref. 20) are based upon a concept of an equiva- 
lent rotor with an elastic hinge offset and are sufficient for 
trim calculations. The concept of equivalent hinge offset employed 
by Bramwell is not similar to Young's more widely used concept of 
equivalent blade because its value is associated only with the fre- 
quency (Ref. 27) . 

(8) The range of thrust coefficients used in the calculations is 
0.005 < C T < 0.015. The practical range is 0.005 < C T < 0.01. 

The range 0.005 to 0.015 was selected so as to include cases which 
could occur during violent maneuvers or gusts. 

(9) The range of frequencies for flap was taken between 

1.05 < uL, < 1.6 
FI 

For lag, the range was selected as 

0.8 < u < 1.5 

Li 

For most cases devoted to studying the flap-lag-type of 
instability, the lag frequency was taken as greater than one to 
avoid the possibility of air or ground resonance (Ref. 28). A 
limited number of cases with lag frequencies below one were 
studied in Ref. 9. 

(10) The range of Locke numbers was taken as 
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( 11 ) 


S-tf-Io 

All calculations were performed for either y = 5, or y = 10. 

The Locke numbers more representative of a hingeless helicopter 
blade are 7 < y < 10 (Ref. 28) . 

Unless otherwise stated, the solidity ratio in all of the calcu- 
lations was taken as 


a = 0.05 

(12) The range for the structural damping in the calculations was 
taken as 

0 < n . n < 0.02 

SF1 SLl 


This range is sufficient to clearly illustrate the effect of 
structural damping. In some modern elastomeric bearings, the 
value of the structural damping can be as high as 0.05. 


7.2 Description of the Methods Used for Obtaining Numerical Results 

The results presented in this section were obtained by two distinct 
methods : 

(a) The expressions obtained by solving Eqs. 3.23 and 3.24 by 
the perturbation method were programmed on a computer, to 
obtain the analytical, solution in numerical form. 

(b) For a certain number of cases, the solutions were obtained 
by direct numerical integration of Eqs. 3.23 and 3.24, using 
the predictor-corrector method (Ref. 29) . 

The reason for numerical integration was twofold. First, it served as 
a convenient way to check the solutions obtained analytically. Second, in some 

cases, the two methods were complementary to each other; i.e., by knowing the 
solutions from numerical integration, it was easier to derive them in analytical 
form. 

It should be mentioned that the a priori knowledge of the limit cycle 
amplitude as obtained from Eq. 5.66 was extremely helpful in obtaining the limit 
cycle by numerical integration. Specifically, the low damping in lag and the 


91 




the time scale = £i p Q in which the amplitude growth is represented, combine 

so that the time required to reach a limit cycle can be of order ip = 600, if 

o 

the initial conditions are not close to the limit cycle value of the amplitudes. 
This can be shown analytically, also. This effect would preclude the use of 
numerical integration for obtaining limit-cycle amplitudes in hover without 
a priori knowledge of the value of the limit-cycle amplitudes. 

Finally, it is of interest to note that the predictor-corrector method 
is equivalent in accuracy and computing time to the fourth order Runge Kutta 
method. At first glance, it appears that the predictor-corrector method should 
be more efficient in terms of computing time because it requires one half the 
number of points per interval of integration as the Runge Kutta method. As 
pointed out by Lapidus (Ref. 30) , the numerical stability of the predictor- 
corrector method for linear problems is worse than that of the Runge Kutta 
method by a factor of two (approximately) . Therefore, in order to obtain the 
same accuracy, twice as many intervals are required in the predictor-corrector 
method. This effect will cancel the apparent advantage of the predictor-corrector 
method due to the smaller number of points per interval. For a nonlinear analy- 
sis, this situation will be somewhat modified, but numerical experiments per- 
formed indicate that the two methods are essentially equivalent in computing 
time when requiring the same amount of accuracy. 

7.3 Results for Hovering Flight 

7.3.1 Stability Boundaries in Hover Without Structural Damping 

From the solution to the linearized problem, described in Section 4, 

with no structural damping (n = ri T1 = 0), stability boundaries resembling 

SFX SL1 

ellipses can be drawn. The value of 9^ specified on the curve is the value 
of collective pitch above which the linear system becomes unstable in hover. 

The unstable combinations of and are given by the area inside of 
the curve. 

Figures 3 and 4 show the stability boundaries with the inflow evaluated 
from Eq. 5.5 for hingeless elastic blades with F 1 , L 1 as given in Table 1. The 
limit cycle amplitudes which would occur when crossing into the post-critical 
range 0 > © c were calculated using the equations and the criteria given in 
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Subsection 5.5.1. The dotted part of the stability boundaries represents combi- 
nations of flap and lag frequencies such that unstable limit cycles occur. In 
order to check upon this prediction of stable and unstable limit cycles, three 
points along the stability boundary (Fig. 4, 0^ = 0.20), denoted by a, b, c 
in Fig. 4 were also checked by numerical integration. The results are shown 
in Figs. 5, 6, and 7. Figures 5 and 6 indicate the stable limit cycles, while 
Fig. 7 shows an unstable limit cycle. 

From Figs. 3 and 4, it is clear that by increasing y, the unstable 
areas enclosed are considerably increased.* The inner curve 0^ = 0.175 in 
Fig. 4 shows approximately the minimum value of 0 c for y = 10 below which 
no instability can occur. 

Figures 8 and 10 show the appropriate stability boundaries for the 
elastic blade where the inflow as represented by Eq. 5.5 is replaced by the 
inflow calculated from Eq. 7.3. 


>o 




2V 0 
a C 



(7.3) 


This inflow relation is equivalent to taking the induced velocity at 3/4R 
of the blade as representative of the constant induced velocity over the whole 
disc. As pointed out in Appendix D, this assumption gives values of C T which 
are not in agreement with those given by momentum theory. The use of Eq. 7.3 
in the computation of the stability boundaries decreases the size of the un- 
stable areas enclosed by the stability boundaries. The physical explanation 
for this effect is clear if the effective angle of attack defined as 

o( - £ 0 — £ ^ o 

$2 


l,= (t\ L jLHF‘ ; 

Vfil tf Fl (Rl ir FI 

is considered, as the physically meaningful quantity. 


(7.4) 


* 

Stability boundaries given for Y = 5 and y = 10 can also be considered to be 
representative of the lift deficiency function of C(k) = 0.5 (with respect to 
Y = 10). (C(k) for this case is not frequency dependent (Ref. 31)). 
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Equation 7.3 yields higher values of X than Eq. 4.6 for the same value 

o 

of 0; therefore, from Eq. 7.4 it is clear that the use of this inflow relation 
yields lower values of a, resulting in a system which is more stable than in 
reality. 

Finally, if Eq. 7.3 is replaced by the assumption that constant induced 
velocity is that which would be obtained by taking the angle of inflow at 
3/4R as representative of the blade, the following relation for X is obtained 



Of the various inflow relations (i.e., Eqs. 5.5, 7.3, and 7.5) the highest 
value of inflow, for a given value of pitch setting, is obtained from Eq. 7.5. 

A typical stability boundary obtained using Eq. 7.5 is shown in Fig. 9. Com- 
parison of Figs. 4, 8, and 9 shows that this assumption results in a further re- 
duction in the unstable area inside the stability boundary. 

Figure 10 shows the comparison of the stability boundaries obtained by 
considering the elastic blade as modeled in this report and comparing it to the 
centrally-hinged, spring-restrained blade for which the stability boundaries 

were obtained by Ormiston and Hodges (Ref. 3). Only the case y = 5, 0 =0.20 

c 

is considered. As seen, the elastic modeling of the blade results in a slight 
increase in the unstable area inside of the stability boundary, while at the 
same time the location of the whole ellipse is shifted in the u> F1 , ^ L1 plane. 

Comparison of the two stability boundaries given in Fig. 10 reveals a 
significant effect of the mode shape in shifting the stability boundaries. 

This result implies that a hinge offset should be included in the centrally- 
hinged, spring restrained rigid blade model of the elastic hingeless blade. 

7.3.2 The Effect of Structural Damping on the 
Stability Boundaries 

Figures 11 and 12 show the effect of structural damping in lag on the 
value of 0^ for the following four cases: 
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03 = 

F10 

1.2 

03 — 1.03861 

LJ.U 

n SFl = 

0.005 

Y = 10 

W F10 

1.2625 

“lio = 1 - 1 

n SFl = 

0.005 

Y = 10 

W F10 

1.200 

“lio = 1 - 1375 

n SFl 

0.005 

Y = 10 

W F10 

1.175 

“lio = 1 * 33319 

n SFl 

0.005 

Y = 10 


As seen from these curves, the structural damping has a very strong 

effect on the value of 0 . The increase in 0 due to the addition of r) , is 

c c SLl 

very strong for small additions of structural damping and levels off around 
% 

n , = 0.015. This is mainly due to the low value of C. /a. 

SLl d 

o 

The effect of the structural damping on the stability boundaries 
is very important because it raises the values of 8^ beyond practical 
values of collective pitch, and consequently (see Fig. 4) it raises the mini- 
mum value of 0 below which no instability in the linear sense can occur, 
c 

The amount of structural damping in flap has no effect on the value of 
0^, except that it changes the third significant figure in 0^ , within 

° < n SFl < °* 02 * 

7.3.3 Limit Cycle Amplitudes 

For the cases where stable limit cycles exist, see Figs. 3 and 4. The 
limit-cycle amplitudes can be obtained from Eq. 5.66. 

From Eq. 5.65, for values of e, e/0^ < 1, the maximum values of X 2 
at their limit cycle value can be obtained approximately from 

M® I ' S(.c. 

where p„ is given by Eq. 5.66 and 
• c • 



(7.6) 
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A typical limit-cycle amplitude- response curve as obtained from Eq. 7.6 
is given in Fig. 13 and is indicated by the full lines. 


In order to check the results obtained from the perturbation method, 
these results were also obtained by numerical integration. As seen from Fig. 13, 
the agreement between the results is quite good. 


Figure 14 shows the effect of the structural damping on p. for a 

Jc.c. 

typical case. As can be seen, the decrease in p. starts to level off 

^ Xt.C. 

around n = 0.015. 

SLi 


To illustrate the effect of the structural damping on the location 
and steepness of the limit cycle, amplitude-response curves, the points cor- 
responding to n. * 0, n . = 0.0025 and n Tl = 0.00625 of Fig. 14 are plotted 
in Fig. 15. These curves were calculated using Eq. 7.6. As seen from Fig. 15, 
the structural damping changes drastically the location of the limit-cycle re- 
sponse curves. Increasing the value of n T . tends also to reduce the steepness 

SLI 

of the limit-cycle-response curves. 


The steepness of the limit-cycle-response curves is an indication of the 
stabilizing effect of the nonlinearities of the system. It also indicates how 
far 0^ can be exceeded before the amplitudes of response become too large to be 
of any practical value. From the results presented in this section, it can be 
seen that the limit-cycle amplitude response curves are quite steep. This 
means that the nonlinearities in the system are weak and their stabilizing effect 
is not strong enough to reduce the amplitudes of limit-cycle response to practi- 
cal levels once the critical condition is exceeded. 


7.4 Results for Forward Flight 
7.4.1 Trim Curves 

As pointed out previously in this report, the effect of forward flight 
can be correctly investigated only when considering the behavior of the rotor 
at a fixed value of C T while varying M . This can be accomplished by requiring 
that the rotor be in a trimmed condition. Using the trim procedure described 
in Appendix F, a subroutine which calculates the trim conditions has been in- 
corporated in the two computer programs (one using numerical integration and 
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one using the analytical expressions of the perturbation method) . 

A typical set of trim curves for = 1.20 and two typical values of 
are shown in Figs. 16 and 17. It is clearly evident from these figures 
that an increase in y is always first accompanied by a decrease in 0 and X 

o 

up to ^ y = 0.2, after which 0 increases quite rapidly. 

7.4.2 Effect of Forward Flight from Numerical Integration 

Before starting to discuss the results obtained from the use of the 
perturbation method, it is instructive to look at a typical case in forward 
flight which has been solved by numerical integration. 


The results showing the amplitude response for trimmed flight at 
C T = 0.00,6 are given in Fig. 18 for various values of y. The range applicability 
of the analysis done in the present report is 

0 < JU, -4 . 30 

For the case considered, the 0 < 0 c ; therefore, the system is below 
critical and the amplitude grows with increasing y. The quantity (X^ av 
plotted in the amplitude response curves is defined by 



X 


imaX 


Xt rnCn 


z 


(7.7) 


It is necessary to use Eq. 7.7 because the natural equilibrium condition 
defined by Eqs. 3.20 and 3.21 is defined for y = 0. Therefore, the equilibrium 
position about which the blade will oscillate in reality will be a function of 
y. This effect is very small in lag (therefore only the max. value of x 2 
plotted) , but considerable in flap. As seen from Fig. 18, the increase in for- 
ward flight speed increases the amplitudes of response. Below the values of y 
the appropriate values of 0 are given for trimmed flight at fixed C T . At 

y = 0.1, 0 is much lower than 0 and the blade response is moderate, and 

c 

‘Xl’av =• X 2 - 

For y = 0.4, 0 = 0.2936 which is only 5% less than 0 c? therefore, the 
typical behavior of the lag degree of freedom near the critical condition mani- 
fests itself by a sudden growth in \ . 
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7.4.3 Effect of Forward Flight on Case B 

1/2 

This case represents the general forward and flight case with y = e 

and a) ^ 1 + ev or w = 1/2 + £V. For this case two possibilities exist: 
c c 

(1) The combination of C T , y and the structural damping in 
lag n SL1 result in a flight condition for which the col- 
lective pitch is below critical (i.e., 0 < 0 ). The blade 

c 

response is given by Eqs. 6.8 and 6.9. This case repre- 
sents essentially the forced response of the system be- 
cause the flutter mode is a decaying oscillation. The 
flutter mode is a decaying oscillation because 

= feft). —T + (M, 

M' M> 

is usually a relatively large negative quantity. 

(2) The combination of C T ,y and low or zero structural 
damping can result in a flight condition for which 
6 > 6 c - In this case, the blade response is de- 
termined from Eqs. 5.74 and 5.70. The stability in 
this case will be determined by < 2R , < 3R , and ic 6R . 

As shown in Subsection 5.5.2, the quantity which de- 
termines the stability is 

Stable limit-cycle oscillations can occur when 
2 

K on +y K > 0 and k < 0. The quantity in Eq. 7.9 

Zj\. O DR JR 

is strongly dependent upon the trim condition through 

2 

the variation of 0, = (0 - 0 )/y with forward flight. 

1 c 

The amplitude growth in this case occurs in the time 

2 

scale <|> = eip = ip. The quantity A governing the 

I o o 

postcritical amplitude response can be written in the 
following functional form as 



(7.8) 


(7.9) 
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The two possibilities mentioned above are illustrated by the numerical 
results given below. A representative case was chosen with the following 
properties . 


“VlO = “fi 1,175 


to = a) = 1.33319 
L10 LI 


0 J = 1.32641 

c 


0 = 0.357523 

c 


y = io, a = 0.05, 


C = 0.012 

Dp 

n SFl = n SLl 

C = 0.01 
T 


0.005 



0 . 01 , 


a = 2TT 


For this case, due to the presence of the structural damping 0 < 0 , 

c 

Figure 19 illustrates the blade response at fixed C T as a function of y. 

At four values of y(y = 0.1; 0.2; 0.25; 0.32) the values of the collective 

pitch 0 as obtained from the trim calculation are also given. The dotted 

lines in Fig. 19 represent the results obtained by using Eqs. 6.8 and 6.9. 

The flap and lag amplitudes were taken as the maximum values occuring between 

250 < ip < 350, the quantity (X, ) is defined by Eq. 7.7. Due to the re- 

1 av 

marks concerning the validity of the uniformly valid expansion made in 
Subsection 6.1, the same curves were also calculated by direct numerical in- 
tegration. The results from the numerical integration are given by the full 
lines in Fig. 19. The agreement between the two sets of curves is quite rea- 
sonable . 


In order to show that Eqs. 6.7 and 6.8 represent correctly the time 
history of the blade response, the time history for the blade response at 
y = 0.25 was obtained by numerical integration and is given in Fig. 20. 

Figure 21 represents the same blade response history as given by Eqs. 6.7 and 
6.8. The two sets of curves are quite similar.* 

Figure 22 shows the blade response amplitudes as evaluated from Eqs. 

6.7 and 6.8 for C T = 0.0078. As can be seen, decreasing reduces the blade 
response levels. 

* 

From a mathematical point of view, the results of the numerical integration should 
be harmonically analyzed, and the coefficients for a Q2 , b^, a i 3 ' b i 3 ' ••• • etc * < 

should be obtained. These should then be compared to the same coefficients as ob- 
tained from the perturbation method. 
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Considering the same case for which the relevant quantities were 
given above, with structural dancing set equal to zero (n = n oT , = 0) , the 
blade response in the postcritical range (0 > 0^) can be considered. Removal 
of the structured, damping results in a new value for the critical value of 
the collective pitch 0 c = 0.20 (see also Fig. 12). The blade response ampli- 
tudes are given in Fig. 23 for C T = 0.0078. The curves given in Fig. 23 were 
evaluated using Eqs. 5.70 and 5.74. As can be seen, the postcritical region 
is characterized by large amplitudes in lag. Although this is a stable limit 
cycle oscillation, the amplitudes in lag are so large that the results do not 
have any practical significance. By comparing the blade response at 0 < @ c 
with the response at 0 > 0^ (Figs. 19 and 23) , it is interesting to note that 
below the critical condition the amplitudes in flap are usually larger than 
in lag. While in the vicinity of 0^ and above it, the lag amplitudes are 
much larger than the flap amplitudes . 

Comparing Fig. 22 to typical blade response in hover. Fig. 13, it is 
seen that forward flight tends to reduce ( |x 2 I )/( IXj I )- 


The effect of trim at fixed C T on the stability of the blade and the 
amplitudes of the blade response are illustrated by Figs. 24 and 25. 

Figure 24 shows a plot of the quantity defined by Eq. 7.9. As can be seen, 
this quantity is considerably affected by y. Figure 25 shows a plot of 
p g as calculated from Eq. 5.71. This quantity which determines the magni- 
tude of the blade response is also considerably affected by the y through the 
requirement of trim at fixed C^. 

7.4.4 Effect of Forward Flight on Case C 

The flutter frequency for this case is o> = 1/2 + £V and y = e. 

c 

In this case, parametric excitation will be the predominant effect. 


From the stability boundaries given for y = 5.0 and y = 10 (Figs. 3 

and 4) , it is clear that the value of 0 c » even with n spl = n gL1 = 0 for this 

case will be so high as to have almost no practical value. Adding a small 

amount of structural damping (ri = 0.005) will increase 0 even further. 

SL1 C 

Thus it appears that the cases of interest would be 0 < 0 , which would re- 

c 

quire the use of the uniformly valid expansion. As mentioned in Subsection 6.1, 
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the validity of this expansion is good for y < 0.15. For such low values of 
forward flight, it is doubtful whether parametric excitation could have a 
significant effect. 

As shown in Subsection 5.5.3, the solution for this case is given by 
Eq. 5.94. It has been shown that in the vicinity of the origin, the stability 
is governed by the real part of s^ where s^ is given by Eq. 5.79a, page 46. 
When Real (s^) < 0, the system is stable and for Real (s^) > 0 the system is 
unstable. Also due to the homogeneous form of the linear part of Eq. 5.52, 
the flutter mode below the critical condition will not be excited, it is also 
shown in Subsection 5.5.3, that if a closed trajectory (limit cycle) exists 
its distance from the origin will be bounded by Eq. 5.92. Also, for stability 
in the nonlinear range, it is required that K < 0. It has also been shown, 
by applying the perturbation method in the neighborhood of the critical con- 
dition, that the location of the equilibrium point and the magnitude of the 
limit-cycle amplitude (if it exists) will depend upon the parametric excita- 
tion. The blade response up to the first order can be approximately written 
as 


X 2 = 2 Rea?. [A,( ykfo, ilii , 

(7.10) 

Due to the large values of 9 for this case, no numerical results are 

c 

given since they would be impractical. 

7,4,5 Effect of Forward Flight on Case D 

1/2 

The flutter frequency for this case is ^ = 1 + ev and y = E 
In this case, parametric excitation has an effect, but the forcing will be pre- 
dominant because the system is excited very close to its resonant frequency. 

Near the critical condition, the flutter mode is determined by Eq. 5.58 
or Eq. 5.98. The stability of the system near the origin depends upon the sign 
of Re(s 1 ), for Re(s 1 ) > 0 the system is unstable and when Re(s 1 ) < 0, the system 
is stable. The linearized system is nonhomogeneous , Eqs. 5.101 and 5.102. 
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Therefore the flutter mode is always excited the same is true for Eq. 5.98. 
The solution for this case is given by Eq. 5.94, and up to the first order, 
the lag motion can be approximated by 

X 2 r 2 u. Real f Aof &-(k Xp-Xoc 

A 2 -. ’ ■ 

It should be noted that the analysis for this case was performed with the 
assumption that C T and y are small. 



For numerical calculations, a representative case with the following 
properties was chosen: 


0 ) 

F10 

= “Vl 

= 1.20 

0) 

= u 

= 1.03861 

L10 

LI 


U) = 1.04146 0 = 0.20 


c 

y = 10, a = 0.05, 


c 

C J = 0.01 a = 2ir 
d 


Sp = °* 012 


n SFl n SLl 
C T = 0.0080 


0.0 


For the value of C chosen above 0 > 0 . 

T c 

Numerical results for this case can be obtained either by integrating 
numerically Eqs. 5.58 or 5.98, and then making use of Eq. 5.94, or by using 
the numerical integration program. For convenience, the second method was 
used. The results showing the blade response amplitudes are given in Fig. 26. 

As can be seen, this case is characterized by large amplitudes of blade re- 
sponse which become very high even for moderate values of forward flight. 
Therefore, flutter frequencies U) c close to u> c = 1 should be avoided in practice. 

Numerical integrations done for this case indicate that the special ex- 
pansion used for Case D, and physically characterized by the beating phenomenon 

as represented by the term e lV,J in Eq. 5.105 (or in Eq. 5.57) , is valid up 
to the advance ratios of y < 0.2. 
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SECTION 8 


GENERAL EQUATIONS FOR NONLINEAR, COUPLED FLAP-LAG-PITCH 
MOTION OF HINGELESS HELICOPTER BLADES 


8.1 Introduction 

Although research in the field of aeroelastic problems of helicopter 
blades, both hingeless and articulated, has been going on for a considerable 
amount of time (Ref. 32) , a complete set of consistently-derived equations 
for nonlinear, coupled flap-lag-pitch motion of hingeless (or even articu- 
lated) helicopter blades is not available in the literature. The author be- 
lieves that the lack of these equations has acted as an obstacle in the de- 
velopment of reliable analytical methods for predicting rotor blade stability 
in the range where the amplitude of the blade motion is large, because most 
of the work in this field was restricted to solving a special system of equa- 
tions derived under a particular set of simplifying assumptions. Thus, the 
basis for a meaningful comparison of the results obtained by various re- 
searchers was not available and the comparisons were usually of a qualitative 
nature. 

It is instructive, therefore, to review the various systems of equa- 
tions available in this field. In this review, both the linear and nonlinear 
equations of motion will be considered. 

The best system of equations is the one derived by Miller and Ellis 
(Ref. 11). This system of equations, which was derived in order to treat 
linear blade response problems, identified all of the important physical 
effects associated with the coupled flap, pitch and lag motion of the blade. 
Another advantage of this derivation is that it also considered distributed 
torsional effects and some finite displacement effects, like effects due to 
large coning angles. Another useful device used in this derivation was the 
modeling of the blade by a concentrated mass at the blade c.g. offset by a 
distance x i from the elastic axis. Finally, the last advantage of this deriva- 
tion is its clear treatment of the derivation of the aerodynamic loads. 

Chronologically, the next system of equations was derived by Houbolt 
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and Brooks (Kef. 10) . The purpose of the authors in Ref. 10 was to derive 
a linear system of equations for coupled flap- lag-pitch motion; as such, 
their derivation was quite successful. They succeeded in deriving an engineer 
ing beam-type theory for treating the elastic deformation in flap-lag and 
torsion of thin elastic rotating beams in a consistent manner. The linear 
inertia loads were also consistently and systematically derived. However, 
the Coriolis inertia load, due to the elastic shortening effect, was not in- 
cluded in their derivation. An additional point in favor of the Houbolt 
and Brooks (Ref. 10) derivation is that it contains some of the elements re- 
quired for extending the derivation into the nonlinear range. Their treatment 
includes finite, spanwise varying, ‘built-in twist and offset between the blade 
elastic axis and the line of cross-sectional centers of gravity. However, in 
these equations no provision was made for treating built-in coning and the 
aerodynamic forces were not treated. The system of equations derived in 
Ref. 10 has been used in a number of papers dealing with blade flutter. 

The equations of motion derived in Ref. 10 have been extended by 
Lemnios in Ref. 33 to include the Coriolis force in the chordwise direction 
due to flapwise bending (elastic shortening effect, see Eq. 3.9). In addi- 
tion, the aerodynamic terms of the equations of motion were treated in detail 
in Ref. 33. The only nonlinear effect taken into account by Lemnios is the 
elastic shortening effect, while many other second-order effects associated 
with the inertia-loading terms have not been considered; thus, his equations 
are not consistent. In Ref. 33, only the offset between the elastic axis and 
the blade center of gravity was considered. Built-in twist and distributed 
torsion were both treated. 

A further improvement of the equations of motion derived by Houbolt 
and Brooks (Ref. 10) was made by Pizialli (Ref. 34) . This system of equations 
although still limited to the linear range, had the following improvements; 

1. Built-in twist as in Ref. 10 

2. Noncoincident and nonstraight cross section; e.g. , 
centroid and elastic axes which do not pass through 
the rotor axis of rotation. 
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3. Details of articulation such as to include: 

(a) Root elastic boundary conditions 

(b) Radial and chordwise hinge offsets 

(c) Hinge inclinations 

(d) Pitch axis offset and inclination 

Another system of equations of motion has been derived by Bielawa 
(Ref. 35) . These were derived for the purpose of investigating the higher- 
order effects due to both elastic root torsion and distributed elastic torsion 
in the spanwise direction. In order to simplify matters, the location of the 
blade elastic axis, the blade center of gravity, and the blade aerodynamic 
center were assumed to be coincident in each cross section of the blade; 
thus a large number of second-order effects vanish. Large elastic displace- 
ments in flap and lag were assumed.* On the other hand, the Coriolis loads 
due to the blade shortening effect were neglected. The effect of large dis- 
placements on the aerodynamic loads was not carefully treated. Thus, no 
clear picture of the various approximations involved could be obtained from 
Ref. 35. The main contribution of Ref. 35 is a careful and detailed treat- 
ment of the various torsional effects. Thus, this work, too, suffers from 
a certain amount of lack of consistency. 

Finally, another detailed linearized version of the coupled flap-lag- 
pitch equations of motion has been derived by Arcidiacono (Ref, 36). The 
equations of motion were derived for linearly twisted rotor blades . The 
motions include flapping and lagging for the articulated blade, as well as 
flatwise, edgewise, and torsional deformations for the articulated and non- 
articulated blade. Fully-coupled aerodynamic forcing functions were also 
derived based on quasi-steady aerodynamic theory. These differential equa- 
tions of motion were also expanded in terms of the uncoupled vibratory nodes 
of the blade. 

In this section, an attempt will be made to derive, consistently and 
systematically, a system of equations of motion for coupled nonlinear flap- 
lag-pitch motion of hingeless helicopter blades. 

The equations will be carefully derived so that the various approxima- 
tions involved in obtaining the elastic, inertia and aerodynamic loads will 
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be made clear. 


It is of importance to note that for a hingeless blade, the equations 
of motion are more complicated than for an articulated one. The reason is 
that for a hingeless blade, the only physically meaningful reference plane 
is the hub plane. In this reference plane, the angle of pitch setting cannot 
be taken as a fixed quantity in time, but will be a time -dependent quantity 
due to the presence of the cyclic pitch. Due to this effect, additional 
inertia terms will appear in the equations of motion. 

For the particular case of an articulated blade, where for convenience, 
blade stability analyses are usually performed in the no feathering plane , 
the general equations derived in this section will still be applicable by 

* ffk 

setting 0=0=0 and replacing the elastic hingeless mode shapes by appropriate 
mode shapes for an articulated blade. 

In the equations which will be derived, only elastic root torsion will 
be treated and no provision for built-in twist will be made. Small angles of 
built-in coning will be included in the analysis. The arguments for neglect- 
ing higher order nonlinear effects will be stated and consistently applied to 
obtain the final form of the equations of motion in a general -form. 

8.2 Basic Assumptions 

The assumption used in the derivation of the equations of motion will 
be given below. For convenience, they are divided into three groups: 

A. Geometrical Assunptions 

These are assumptions regarding the geometry of the blade and its 
attachment to the hub. 

(1) The elastic blade is attached to the hub at some offset 
distance e 1 from the axis of rotation (see Fig. 27) 

(2) At its root, the blade can have a built-in coning 
angle 3^. In addition, the feathering axis can have 
an inclination 3 (measured in a vertical plane) 
with respect to the hub plane. It is assumed that 
these angles are sufficiently small so that the 
assumptions sin3=3, cos3=l (and the same for 3 ) 
are valid, unless otherwise stated. 
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(3) The feathering hinge is inboard of the "virtual" 
flap and lag hinge (see Figs. 27a and 27b) 

(4) The blade is initially straight. There is no angle 
of built-in twist. 

(5) The cross sections of the blade are assumed to be 
symmetrical about the major principal axis (see 
Fig. 28). 

(6) The equations will be derived in such a manner as to 
enable the following choice in the location of the 
elastic axis (E.A.), aerodynamic center (A.C.), axis 
of twist [ (A. T . ) or feathering axis] and blade cross 
section C.G. (C.G.): 

(a) E.A. , A.C. , A.T., coincident with C.G. offset. 

(b) E.A. and A.T. coincident, both offset with 
respect to A.C. and C.G. which are coincident. 

B. Elastic Assumptions 

. These are assumptions regarding the elastic properties of 
the blade and its deformations. 

(7) During the deformations, cross sections are assumed to 
remain plane and normal to the elastic axis. 

(8) Shear is neglected. 

(9) The Blade is a thin flexible blade attached to the hub. 

(10) The blade can bend in two directions normal to the 
elastic axis. 

(11) The blade is torsionally rigid, except at the hub. 

The twisting of the blade is represented by a root 
torsion denoted by $. 

(12) The deflections of the blade are moderately small so 

2 

that terms of 0(2^) can be neglected when compared 
to 1. (e^ being the order of magnitude of the de- 

flection.) With this assumption, it is sufficient 
to retain only the linear and the second-order non- 
linear terms in the equations of motion. All third 
order terms can be neglected. 
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(13) Only a linear treatment of the elastic restoring forces 
will be considered; i.e., large deflections will have 
only a small effect on the tension due to elastic ef-^- 

f ects in the blade since one of its ends is free . This 
assumption is consistent with neglecting terms of O(e^). 

C. Aerodynamic Assumptions 

These cure the assumptions made in calculating the airloads acting on 
the blade. 

(14) Two-dimensional quasi-steady aerodynamic loads 
are used. 

(15) Apparent mass effects in the aerodynamic loads 
are neglected. 

(16) Stall, compressibility and reversed flow effects are 
neglected; thus, the aerodynamic loads are valid for 
V < 0.3. 

8.3 Displacements, Coordinate Systems, and Coordinate 
Transformations 

8.3.1 Sample Displacement Fields 

In the present analysis, the feathering axis of the blade will be 
assumed to have an orientation given by an angle 3 with respect to the hub 
plane (see Fig. 27a). Therefore, it is important to define sample dis- 
placement fields which can be used latter with the general equations of 
motion, carefully and in detail. 

A. First Displacement Field 

Consider first the case where there is no elastic root torsion. The 
geometry for this case is shown in Fig. 27a. The required quantities are the 
displacements of a point A (A' in deformed state) located on the elastic axis 
of the blade. For this purpose, only the hub plane coordinate system, attached 
to the blade and rotating with angular velocity (2 is required. According to 
assumption (4) Subsection 8.2, the initial position of point A on the elastic 
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axis coincides with the x-axis. Also note that in its undeformed state, the 
blade cross section at A has a collective pitch setting of 0. From the 
geometry of Fig. 27a* 


VO = X 0 -Strip + We COS ji 


( 8 . 1 ) 


V = 


% 


( 8 . 2 ) 


A 

U = - 


w« - *. ( ( I - C°sp) - i / [ ( * (Jgj J JL,, (8. 


3) 


Now assumption (2), Subseciton 3.2, will be used. Thus 


W = x 0 p +• w € 


(8.1a) 


A 

V = Vg 


(8.2a) 


U = -V)*p 




(8.3a) 


In Eq. 8.3, all quantities are small second-order quantities, there- 
fore the approximation cos 8 = 1 is not permissible. The last term in Eq. 8.3 

is the elastic shortening effect, also given in Eq. 3.9. The quantity w is 

e 

measured perpendicular to the real position of the undeformed blade (see 
Fig. 27a) . 


Finally, the displacements due to the root torsion will be taken into 
account. The assumption will be made that the torsional displacement occurs 
after the elastic flap and lag displacements (i.e., w^, v^) have occurred. 

If the torsional displacement around the feathering axis is $ (see Fig. 27a) , 
then the displacements due to torsion are (4> is a small angle) 


Vi T = <1 

v T - - 


(8.4) 


Note that the inextensibility assumption is used in Eq. 8.3. In general, this 
is not absolutely necessary. 
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Thus , the total displacements of a point on the elastic axis of the blade can 
be written as 

W = W + W T = f V/ € + 

(8 * 5) 

V' = v -*■ -vr T - V e - Wg 4* 



It is important to realize that if the elastic flapping and lagging 
displacements are not assumed to occur before the root torsional elastic dis- 
placement, then the additional displacements due to the coupling between 
the steady-state elastic flap and lag with $ will have a different form. 

B. Second Displacement Field 

Another possible displacement field is one in which the feathering axis 
has an orientation determined by the angle 3 with respect to the hub plane, 
while in addition, the blade has a built-in preconing angle given by 3^ with 
respect to the feathering axis. This displacement field is schematically 
shown in Fig. 27b. 

From considerations identical to those applied in the previous sample 
displacement field 


u 


- v<e 
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(/» V*-) 


A f 
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3 i vi fr 


Ir = 



W 


X a (^+p p ) +- Vie 


( 8 . 6 ) 


Again, the displacements due to torsion are given by 
- -X «f P § COS& - V) c § 

(8.7) 

VI T ~ 


Thus, the total displacements of a point on the elastic axis of the blade 
is given by 


U - U - 


- W. r ) - £ (f */3 r ) 


nr = v •+ v. 


— We § - Xo^f, ( fr+|c 05 0- J 


( 8 . 8 ) 


V\J 


X 0 (^3+| 3 f ) + We + :^ + wr 7 
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In Eqs. 8.7 and 8.8, the relation 
St»(9-+$) - + $corfr 

has been used. For the particular case of 3^ ~ 0, Eqs. 8.8 reduce to Eqs. 8.5 


8.3.2 Coordinate Systems and Coordinate Transformations 


In order to handle, in a convenient manner, the various inertia and 
aerodynamic loads derived in this section, the vector method will be used. 
As will be seen below, six different coordinate systems are required to 
describe the various quantities. These are: 

1. The hub plane coordinate system, shown in Fig. 27a rotating 
at constant angular velocity with the blade. The z-axis 
coincides with the axis of rotation; the x-axis coincides 
with the assumed initial position of the elastic axis 
and the y-axis is perpendicular to the x-z plane. The 
unit vectors of this coordinate system are denoted by i, j, k. 


2 . 


3. 


4. 


A second coordinate system shown in Fig. 29a. Its unit 
vectors are denoted by i , j , Jr. The i, unit vector 
coincides with the orientation of the feathering axis , 
the j, j axes are parallel and k. is perpendicular to 

-V <vl 


ii and Ji- 

A third coordinate system with unit vectors i , j , k 

*v 2 ~2 ~2 

is shown in Fig. 29b. It is obtained by rotating the 
i-T / j.» k system by an angle 0 + ‘P around the i, axis. 

'V X -vX /vX /v X 

A fourth coordinate system with unit vectors I , J , K 

~1 ~1 ~1 

is defined in Subsection 1.1, Appendix I. 


5. A fifth coordinate system will be attached to the cross 

section of the blade at its elastic axis. It moves with 

the cross section as the blade deflects. The unit vectors 

for this system are: I - tangential to the deformed elastic 

z 

axis of the blade; J , normal to the deformed elastic axis 

2 

and parallel to the hub plane, and K normal to I and J . 

*v 2 /v2 2 
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The projections of these unit vectors, at a point P 
on the elastic axis, are schematically shown in Fig. 27. 

6. A sixth coordinate system similar to the previous one 

and rotated about the £ axis by an angle 6+0. The unit 

vectors are I _ coincident with 1^; Jo coincident with 
3 2 

the blade chord, and K„ normal to I, and J„. These unit 

~ 3 ~ 3 ~ 3 

vectors are schematically shown in Figs. 27 and 28. 

According to the proof of equivalence of rotations given 
in Appendix I, within the approximations inherent in the 
various coordinate transformations, the unit vector 

~3 

can be considered a good approximation to the blade-chord 
orientation in space after the deformations have occurred. 

The relations between the various coordinate systems are given in 
detail in Appendix I. 

8.4 Derivation of the Inertia Loads and Moments 

8.4.1 Derivation of the Inertia Loads in the x,y,z 
Directions 

The inertial loads will be derived using the blade model shown in 
Fig. 28. The blade is assumed to be represented by its mass-per-unit 
span m, concentrated at the blade cross section C.G., and offset by a dis- 
tance Xj. from the elastic axis. This model, used first by Miller (Ref. 11), 
is much more convenient than the distributed mass used by Houbolt (Ref. 10). 
On the other hand, care must be taken to correct the propeller moment and 
the rotational inertial term in q^, by the missing cross-sectional polar 
moment of inertia. 

Let x, y, and z denote the undeformed position of the cross-sectional 
blade C.G., while its position in the deformed state will be given by x^, y^, 
and z^. Then the following relations can be written: 
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x =■ x„ + e, 

^ = x x COS & 

T: - X x & 



(8.9) 


X. = X. 


e, +- u - 


52T 

2* 


(*-*> 


- — (?,-w) 

3X„ ' ' 7 


( 8 . 10 ) 


= 3T + 


x l co5[e- ■*■$) = v ♦ (x i cqt^-x j § sit) ©) 


( 8 . 11 ) 


2, = w +■ 


x s»S)(< 9 + §) = vl + x l + $cos&) 


( 8 . 12 ) 


Note that the blade cross section is not in the y-z plane, but is 
normal to the elastic axis. Only the E.A. point is in the y-z plane. This 
is the reason for the last two terms of Eq. 8.10. In addition, it will 
always be assumed that 4> is small, so that 


Cos(fr-t$) = cos©- -§sr»>©- 

sf* (©-+$) = SVT> &■ $ COS S' 


(8.13) 


From Eqs. 8.10 through 8.12, the position vector, the blade center of 
gravity in the deformed position is given by 


R 

/V 


A 


t X. 




(8.14) 


From elementary mechanics (Ref. 37) , the acceleration in inertial space 
can be written as 


Ct = fL„ . +• 2. U) * It R 

/%/ H <v ^ 


U)»R r 4-CJ » ( U> * ( 8 -15) 
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Due to steady rotation, the third term of Eq. 8.15 is zero. From the 

geometry 3/3x = 3/3x , using Eqs. 8.14 and 8.15, 

o 


= ^(x,- x ,n i-zij'Si) +2X.SI) 


(8.15a) 


From Eqs. 8.15 and 8.10 through Eq. 8.12: 

+ —■ x i cos(9- + f )(©■ + 4) Z + x x sm(e-f £) 


- n$l* x z si7 '( & + $) “ 2 lot «>*(* + §)(& + $) 

+ X sin (&+$)( 0 + $) 2 - UtL X cos(& +$)( & + $) 


- n 2 [x 0 + e i + o - x i cos(fr+§)- x r s)n(e-+$)J 

- z il [ v - x x sin ($- + $)(£+$) J ^ 

+ j > v- - x j cos( ©■+$)(§■ + §) - x t sin (§'+•<!)(& +-$ ) 

/v t 

- Jl x £ V + X r C 0 £ (& + #) J + 2 51 [u x x cos(©-+$) 

+ XSin(& + §) ($• + §) - X- Si"n( fr + §) 
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3W 

3X, 


X t COs(0' + §J( & + § ) J 


w-x i sm(e-i-'|j(&+$) X + X x c0i(fr +$)( & + $) 


} 


(8.16) 


The components of the acceleration in the x,y,z directions, respec- 
tively, can be obtained from Eq. 8. 16. After substituting Eqs. 8.13, the 
complete expressions for a^, a^, a^ are given in Eqs. J.l through J.3 of 
Appendix J. 

From the form of these expressions, it is clear that a considerable 
number of terms are negligible. The considerations for neglecting terms are 
given below. 

Denoting the order of magnitude of the displacements by 0{e^): 

1. v,w and their derivatives are of 0( e D ) 

2. u from Eq. 8.6 is usually a second-order 

2 

quantity; thus, u ^ 0 (e d ) 

3. <J> is usually a small quantity (0 < $ < 3°) ; thus is 

usually 4> < 3v/3x , dw/3x 

3 ° 2 °2 * 2 ** 

Therefore CKe^) < <P < 0(e D ). Thus, $$,$ 

can be considered negligible when compared to other 

second-order quantities 

4. It is reasonable to assume that the order of magnitude of 

„ Or 

sin 0 = 0(e D ) . Thus 

o(£.)Mn e- S' 0(6,*) 
o(e B )e$ = 0 ( e^) 
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5. 


From the first part of this report it was found that the 

cyclic pitch is usually bounded by 0 < 0.5 6 , where 0 , 

o o 

is the constant collective pitch setting. Thus, it is 

2 

reasonable to say that the order of magnitude of 0 (3w/3x) , 
0^w is given by: 

e l o(e D )zo (£ b 3 ) 

6. From the first part of this report, it was found that in 
the vicinity of the critical region (0 = 0^) , the lag 

displacements were 3 to 5 times larger than the flap 
displacements. Therefore, when neglecting terms in the 
following equations, it will be usually assumed that 

Vih &■ = vi cos O' 
and lr Cos & > VI si'-b ©- 


All quantities, which according to considerations 1 through 6 

3 

above are of CKe^), will be neglected, unless otherwise stated. 

Neglecting almost all O(e^) in Eq. J.l, the inertia load in the x- 
direction will be given by 


p -ft 

XI * 

[Because L (see Eq. 8.56) is assumed to be zero] 


(8.17) 


From Eqs. J.l and 8.17, using the nondimensional time derivative, 
Eq. 8.18 is obtained 


_ d\r_ 




— 


siv>e- + ( C0S& -t 2^ SW O') 4- 2f • 

\3X 0 dX 0 / J 


-'vnilY f $ - a!w CO£&+ (^.co^ - |3r J 


117 



+ 2$ sii.9 i. cos&)+# fit s«9- - 


fall si^e- - 2!i!_ cos©) l +V(;|y»i’>fr 


- 2 51 X t | 2. 




— 3^ cos 8^ + 2 +5>cos©)J' 


(8.18) 

In the last expression, the first two brackets [ ] and { } represent 
the linear inertia load which was also obtained by Houbolt (Ref. 10) . The 
second { } bracket represents the nonlinear terms, while the last { Jbracket 
represents the additional inertia load due to the presence of cyclic pitch. 

It should be emphasized that a considerable number of terms in Eq. 8.18 are 
still negligible. These are retained for the sake of generality. 

In a similar manner, neglecting the 0 (E^) terms of Eq. J.2, the inertia 
load in the y-direction will be given by 

T^i = ' 7 ” a > ( 8 -i 9 > 

From Eqs. 8.19 and J.2: 

P = - nmSl 3 ' ( Tr + 2.1/ 

T ‘ 

— ^ cos &■ - $ si^ 


- >) - mXjfl* Si*) fr 

- 2 ( Cos* -*?*Siv>e) 1 

\?x 0 iy M 0 QV J J 
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+ j~-^*+ 2 S’!) co^©- - eT ( simfr + § C.OS 9 -) + 2 (| 7; S, ' T) ^’ 


— l^cosO")® 4 c os ©-1 

^ X o ' - 1 

( 8 . 20 ) 

Again, the first group of terms in Eq. 8.20 are the linear terms 
identical to Houbolt’s (Ref. 10) . The other two groups are the nonlinear 
terms and those due to cyclic pitch. 

The total load per unit span in the y-direction is given by 


r, m L a . 


3 u ' 'V 

Similarly, using Eq. J.3 and 


l 


- -ona. 


the inertia load in the z-direction is obtained 

^ = - mnSl ^ w +• $ cosfr') - yr\X x Sl & (cos&- $ sm &) 

The total load per unit length in the z-direction is given by 

t ■ L .r*r„ 


( 8 . 21 ) 


(8.22) 


(8.23) 


(8.24) 


It should be emphasized that the total loads are always assumed to be 
acting at the elastic axis of the blade cross section. In reality, the 
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inertia loads act at the blade cross section center of gravity. Therefore, 
at the elastic axis, these loads will also generate a system of inertia 
moments . 

8.4.2 Derivation of the Feathering Moment at 
the Blade Root 

In deriving the feathering moment due to the complete system of loads 
acting on the blade, the various moment arms must be carefully considered. 
According to assumption (3), Subsection 8.2, the pitch bearing (or the 
feathering hinge) is inboard of the "virtual" flap or lag hinges. Geometri- 
cally , its location will be taken at x = e^ (or x^ = 0) . 

In writting the moments, the i, , j,, k, coordinate system defined in 

>vl ~1 ^1 

Subsection 8.3 is useful. From Eq. 1.1, assuming that 3 is small 

i sr L, - M, A 

«V ( 

/ o on 

M. * a +• 

~ — f -V, 1 

The moment vector about the elastic axis in the hub plane coordinate 
system is 

? = i + (8 * 26) 

'V ^ 

From Eqs. 8.25 and 8.26 




as 


or 


? - (8 - 27> 
The total loading per unit length at the elastic axis can be written 


r = ci + n 
t + r,i. 

~ »■ I /V 


(8.28) 


(8.29) 
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The position vector of a point on the elastic axis of the deformed 
blade, with respect to the pitch bearing, can be written as (see Fig. 27a) 


H, 


= jl (x D +u) +jv + 'kvj 


(8.30) 


Or using Eqs. 8.25 


R - t, 

e ~ 


Rx^u) ^wj ■hj.'lr + W-f X # fU)A J 


(8.31) 


The moment of the loads given by Eq. 8.29 about the point represented 
by the pitch bearing is given by 


jr ' £•*/> v.(-) 


(8.32) 


The required moment about the feathering axis at the blade root is 
given by the i^ component of Eq. 8.32 


~p O'."')] (8 - 33) 

The total feathering moment will be given by 


M 


T.F.fi 




where 


(8.34) 


h = l + n* (8 - 35 > 

In order to evaluate the feathering moment, the quantities M p ^, q^ 
are required. The derivation of these quantities, together with the various 
approximations involved, is given in Appendix K. 


From Eqs. K.6, K.10, K.14, 8.35 and 8.34, M 

its final form. Note that between 0 < x < e , L , 

1 z T 

are all zero. Thus 


T.F.A. 

L , u, 

y r 


can be written in 


v, w, k , x, , and m 
m I 
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tf TFA -j ♦itfvc.se-) 

+ fsm»(«-px.) -f%m»(w-px.) 1 2(w-p.)(^X- C0 S » 


£0L sm»)l 
3X.9Y 


Ax„ '5^1,4 - Sl'If ( sm S cos S' *Jcoi2&J 


WXt ^( smfr»§«s&) v + m x i Sl'[Tr t (siv»*$CK») 


-*v3co s9 * 2.l(s;~mcos») - Y‘«9 
- a ^ si ” 1& ' 2$ 3^Jt WS ^J 

+ mx >M <»**■ * |r ^'^b ~ wsi<9 / 3 


si^e-cose-al \ <A-x 0 -IfStd- - ( I, -I )fr SI* siv&cosfr 

1 J J t W v 


*4- 


onil l x :$>cos& - mx*Sl'[-2 (jj'/ 1 '* 9 ' 


- c°s&) frsi'*) &] kx K 


(8.36) 
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where 1^ is the feathering moment of inertia of the blade 


x. 

( k* + Xj ^ rr\cL )( 0 


and 


A, 

I. - 

ft 


Thus, 


(8. 37) 


(8.38) 


X t 'matx 0 = 1^ ~ I „ (8.39) 


8.5 Derivation of the Aerodynamic Loads 

8.5.1 Approximations Made in Deriving the Airloads 

The purpose of this subsection is to show the various approximations 
involved in evaluating the aerodynamic loads acting on a section of the 
blade. In calculating these loads, assumptions (13) through (15) of Sub- 
section 8.2 will be used. 

Consider the deformed position of the blade in the hub plane coordi- 
nate system (x, y, z system, see Figs. 27 and 28). The position of the 
elastic axis is given by the position vector 

+ W/k (8.40) 

^ c v ' /V 0 rV 

-N/ 

For air-load calculation, the displacement u in the axial direction 
(which is a second-order quantity, anyway) is unimportant and will be neglected 
in this treatment. Thus 

jje = (*. +*•)}, + + (8 - 41) 
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The velocity of this point (the elastic axis) , can be written as 

V EA = fl E + W x fl e (8.42) 

<v /w 'V ^ 

Or, using Eq. 8.41 

= -IrJK * £il(X 0 + € «) + ^ (8.43) 

nJ 

From the geometry shown in Fig. 2, the velocity of the air in the hub 
plane coordinate system is given by 


V 


pSlRcosYi — h&R srn'yj- 


/V 


(Note that in this report the induced velocity is positive down.) 


(8.44) 


From Eqs. 8.43 and 8.44, the flow relative to the blade at the elastic 
axis is given by 


U = - L ^SlRcos'f + Tt-Sl J ~ j- £^JlRsiviV 

+ ir + si (x a + e,)J wj 


(8.45) 


In order to transform the flow relative to the blade, as given by 

Eq. 8.45, into a physically meaningful form, it is convenient to use the 

I_, J_, K„ attached to the elastic axis of the blade and defined previously 
~2 ~2 

in Subsection 8.3.2. 


Note that it is important to use this coordinate system because accord- 
ing to assumption (7), Subsection 8.2, the cross section of the blade is 
normal to the elastic axis. In this plane, the instantaneous geometric angle 

of pitch is given by the angle between the vectors J and J. or (0+$) . 

~2 ~3 

From the relation between the various unit vectors, Eqs. 1.5, the 
following relation can be written. 
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(8.46) 


* - *I‘ 


k - ' sW 




+ 1C 

J. a. % 


From Eqs. 8.45 and 8.46, the flow relative to the blade can be expressed 
in terms of the I , J , K system. Thus 

(8.47) 

or 


w ■ + J* + \ ^ 

For convenience, a picture of the blade and flow geometry resulting 
from these considerations is given in Fig. 30. 


(8.48) 


The total load on the blade can be written symbolically as 


L. L~ kv + 


(8.49) 


From two-dimensional quasi-steady aerodynamics, it is well known 
(Ref. 11) that the quasi-steady lift at the aerodynamic center (c/4 point) 
is due to the angle of attack at the rear neutral point (3/4 c) . At the rear 
neutral point the induced velocity can be written as 
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Thus 


L n s jX* k *[(•**>- ^ + — ■ J <8 ‘ 


50) 




(8.51) 


Where the following approximations were made in writting Eqs. 8.50 and 8.51 

L»-' (i^) ? 

V 1 U Y. 

i.e., the angle of inflow is small and 

both of which are reasonable approximations. Also, the component of the 
velocity parallel to the blade deformed elastic axis will have no effect on 
the aerodynamic loads. 

Prom Eqs . 8.49 and 1.5, L_ can be rewritten in the i , , j , k coordinate 

A ~ 1 -v- 

system. Thus 


L* -~J:[ L z^wf + J + <f L it + 


(8.52) 


According to the symbol convention of this report, L can also be 

A 


written as 


La - i L *r^ L 3r + 4 L 


Et 


(8.53) 


Comparing the last two equations, it is clear that 
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(0.54) 


L 






(8.55) 


L 


XT 


I 

2 i 3X L y r iX 


(8.56) 


It should be noted that Eqs. 8.54 through 8.56 are a direct result of 
the approximations inherent in the coordinate trams formation, Eq. 1.7. 

In the present report, the quantity L will be neglected in the calcu- 

XT 

lation of p because it is a small quantity when compared with the centrifugal 
force acting in the opposite direction. This is a nonconservative assumption 
because the L^load per unit length tends to reduce the stabilizing effect 
of the centrifugal force. 

Finally, the aerodynamic moment along the I direction can be symboli- 
cally written as 


h, = n x I t (8.57) 

From Eqs. 8.57 and 8.46 

These two moments in the j, k direction will introduce additional small 
bending moments in the flapwise and chordwise directions. These will also be 
neglected in the present report. Thus, the aerodynamic moment can be written 
as 


M 


6 



(8.59) 


Finally, approximate expressions for the velocity components, as given 
by Eq. 8.47, are required. Consider Eq. 8.47: u^ is not required in this 
analysis and will be left unchanged. For low values of y [ assumption (16) 
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can be con- 


of Subsection 8.2] , the first two terms in the expression for u 
sidered small when compared with the other terms. Thus 

and 

<< &(*•*<) 

so that 

I 

* 

Similarly, in the expression for u^* v(3w/3x)i2 < vQ and will be 
neglected so that 

n S , 8 . 61 ) 

3X ' 

Comparing Eqs. 8.60 and 8.61 with Eqs. 3.14 and 3.15, it can be seen 
that they are identical. Thus 


Up ~ u Zz 

(8.62) 

U T = " U Yx 


8.5.2 Derivation of the Aerodynamic Loads and 
Moments per Unit Span 

Originally, it was intended to derive the air-laods per unit span in 
a general form and use the quasi-steady assumption (assumption (14) , Sub- 
section 8.2) C(k) = 1 only in the actual numerical calculation. Unfortunately, 
the presence of the cyclic pitch, constant pitch, and time dependent motion at 
flutter frequency would require a formal splitting of the aerodynamic loading 
into these three distinct groups. Thus, in order to remove this unncessary 
complication, it was decided to invoke the quasi-steady assumption right at 
the beginning. 


In order to clarify the various assuinptions associated with the 
quasi-steady approximation and neglecting the apparent mass effects, the 
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results of two-dimensional unsteady aerodynamics will be given below. For the 
case of a two-dimensional airfoil, in unsteady motion, determined by a time- 
dependent angle of attack a and a downward displacement of the elastic axis h, 
(see Fig. 8.1), the unsteady load and moment per unit span can be written as 
(Ref. 12) 



FIG. 8.1 


L = L„ c(4) + f + V, A -(x„ - 'A J 


(8.63) 


where is the quasi-steady lift given by 


L q r [X* ■*" ^ x *) * J 


(8.64) 
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and 


li = c(k)L Q x„+$ a(bR) 3 (x,-.5-tR)[i;'-(x,-.rtR)iiJ 

-j.f V ^L U -*J - S A £ ( u /^ 


(8.65) 


Here the assumption has been made that the aerodynamic center is at C/4. 
Setting C(k) = 1, as justified in Ref. 11, and neglecting the apparent mass 
terms associated with h and a yields 


L = aj. () V,UR[v (| oL + b + (|-X,)bRiJ 


( 8 . 66 ) 


M * *f,,V ( t [v,* + K J + <XJ. V, ( tR)^ I - X, ) Lr i ( X, - . 


(8.67) 


where x. = x/bR. 
A A 


Replacing h = - u , V = U , a = 8 + <J> and a = 0 + $ 

Jr A 1 

Eqs. 8.66 and 8.67 become (for small angles of inflow) 


l2(0 + I) 


*f, u T -% +(| -x,)ua(»* i)J 

tv = ^rHl \ (>(*♦*) -'v] + ( , -\)M lRft (^J 


( 8 . 68 ) 


(8.69) 


Finally, the loading in the y-direction can be obtained from Fig. 2 

L=-^L - ptRC, U* 

y. * *• T (8.70) 


In Eq. 8.70 it is implied that 
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and 


-I Uf> *v l If, 

v; ' u T 

co S [t«. *’(■£)] = 1 

Prom Eqs. 8.68 and 8.70 

Lj. •‘-S» a - h *{' / r L U r ( 9 + ^- U fJ 

+ u t] (8.71) 

8.6 The Equations of Motion 

8.6.1 The Elastic Restoring Forces 

According to assumption (13) Subsection 8.2, a linear treatment of the 

elastic restoring forces will be considered sufficient. Such a treatment has 

been derived by Houbolt and Brooks (Ref. 10) . In writing the equations of 

equilibrium for the beam element, another convenient assumption will be made, 

following Ref. 10. It will be assumed that the element is cut by slices 

perpendicular to the hub plane. This assumption is one of convenience (in 

calculating the loads and moments, it was always assumed that the cross 

section is normal to the elastic axis) , and will have a negligible effect on 

the results. Finally, in the beam equilibrium equations in the y and z- 

directions, the terras 3q /3x , 3q /3x will be neglected. Thus the equations 

y o z o 

of equilibrium can be written as 



£*{ [(«),■ *■•->**] 



Si?) 9 
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■ax 1 


[H-( ei )>^c k6 ^i + [( EI ) 3 s^e.(Bi 


~ T e , 


*] -£J>WsfJ = r. 


(8.72) 


'o ~ ^ 

Note, that only the elastic deformation contributes to the elastic 

/ 

energy of the system. Therefore v , w defined in Fig. 27 are used in 

e e 

Egs. 8.72. 

The elastic restoring moment about the feathering axis can be assumed 

r 2 

to be represented by a spring constant = If U) q . 

8.6.2 Blade Equations of Motion 

The equation of motion in the z-direction is obtained by combining 
Eqs. 8.72, 8.23, and 8.68. 

The equation of motion in the y-direction is obtained by combining 
Eqs. 8.72, 8.20 and 8.71. 

The feathering equation of motion can be symbolically written as 




(8.73) 


where M is given by Eqs. 8.36, 8.68, 8.69, and 8.71. 

1 • r • A • 
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SECTION 9 


NONLINEAR AND LINEARIZED EQUATIONS FOR THE COUPLED FLAP-LAG-PITCH 
MOTION OF HINGELESS BLADES IN HOVERING FLIGHT 

9.1 Introduction 

In this section, the general equations of motion derived in the previ- 
ous section will be specialized to the case of coupled flap- lag-pitch motion 
in hovering flight. Using Galerkin's method for the spatial variable, the 
partial differential equations will be reduced to a system of nonlinear ordi- 
nary differential equations. In this process, no attempt will be made to in- 
clude an arbitrary number of elastic modes (as was done previously for the 
case of flap lag) ; thus only one elastic mode for each elastic degree of 
freedom (i.e., flap and lag) will be used. Due to the considerable amount of 
algebraic manipulations involved , this process of reduction is given in a 
most concise form possible. Where details are required, they are given in 
the various appendices. 

Based on physical reasoning, it is shown that the nonlinear effects 
in the feathering equation must be treated in a different manner than those 
associated with the flap and lag equations. Therefore, the treatment of the 
flap and lag equations is separated from the treatment of the feathering equa- 
tion . 

The resulting system of coupled ordinary differential equations is 
linearized about a natural equilibrium position, which is taken to be the 
equivalent linear steady-state position. In this process, various nonlinear 
effects are transformed into coupling effects. Thus the linearized equations 
contain many new terms when compared with those which would have been obtained 
from a purely linear treatment. 

Next, the linearized equations of motion are transformed into a system 
which can be conveniently solved by defining various flutter derivatives 
(Refs. 11 and 39) . 

Finally, by adding the nonlinear terms to the linearized equations of 
motion, the complete coupled nonlinear equations of motion are written in a 
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form suitable for numerical integration. 

9.2 Assumptions Made in Reducing the General Equations 

In order to reduce the number of terms in the general equations , some 
assumptions must be made. These assumptions, given below, are only a matter 
of convenience, and will have little effect on the mechanism of instability. 

(a) In order to simplify the inertia loads, it is assumed that 
the offset between the elastic axis and the blade cross- 
sectional center of gravity is zero, i.e., x^ = 0. 

(b) The distance e between the area centroid of the tensile 

A 

member and the blade elastic axis is also taken as e =0. 

A 

(c) In order to describe the displacements of the blade, dis- 
placement field B, given by Eqs. 8.8 and shown in Fig. 27b, 
will be used. 

(d) It is assumed that the flap and lag motions of the blade 
can be represented with a sufficient degree of accuracy 
by using one elastic mode in each of these degrees of 
freedom; thus 

(9.1) 

(See Figs. 27a or 27b for the definition of w , v .) 

e e 

The mode shapes used in Eqs. 9.1 represent the appropriate 
mode shapes of a rotating beam without the effect of elastic 
coupling, i.e., at 0=0. 

(e) The boundary conditions for v , w axe taken as the usual 

e e 

ones for a hingeless blade, i.e., Eqs. 3.6 are assumed to 

apply to v , w replacing v and w. 
e e 

(f) The angle of preconing 8 and the built-in coning angle 

P 

of the feathering axis 8 cure assumed to be small. 
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9.3 The Equations of Motion in Flap and Lag 

From Eqs. 8.8 and 9.1, the displacement field can be written as 




(9.2) 


= ~M,f t -1*1,3,$ ~ x 0 i^ p (smd + § cosft) 


(9.3) 


w = /x 0 (p+fr) + - XA,y, $ 


(9.4) 


From Eqs. 8.17, 8.18, 8.20, 8.21, 8.23, and 8.24, assuming that x = 0, 
the loads can be written as 


p - - 21 » - 'vn.Sl* flf- ( x„ +e, + **.) - 
r x 0*. L 


(9.5) 


f> - Ly r + 2 U - '*) - $ 

3 


Slv- C 


(9.6) 


Sl 


** 


- L — lmJ2*W - °l -52. VJ e 

OSF 


(9.7) 

— r r i-*L, *v ^ *c 

■ e *r 

where the effect of a viscous type of structural damping has been included in 


the last two relations. Note that only the elastic part of the displacement 
will dissipate energy through structural damping. 

Denoting the elastic coupling effect by Ec^ and Ec^ as defined in 
Appendix L by Eqs. L. 3 and L.4, and using the assumption that e = 0, 

A 

Eqs. 8.72 can be rewritten as 


k 1 f> 1} > + E ‘0 " Ect ™/ = £ 


+ o ( t 

OX. Ox. 


(9.8) 
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J_ f P 9%/e 
9X* 



(9.9) 


Substituting Eqs. 9.2 through 9.7 into Eqs. 9.8 and 9.9 yields a 
system of equations on which Galerkin's method is applied. The application 
of Galerkin's method consists of the following steps: 

** 

(a) The u , u terms occurring in Eq. 9.5 are neglected. 


2 — 

(b) The flap equation is multiplied by n,& dx and the terms 

1 o 

associated with the inertia and elastic forces are inte- 
grated between = 0 and = i, while the aerodynamic 
loading terms are integrated between x q = A and x^ = B, 
where A and B are the tip loss factors. 


(c) 


(d) 


The resulting equation is nondimens ionali zed by 
2 

dividing by I ft . 

b 


2 — 

The lag equation is multiplied by dx Q , then steps 
(b) and (c) , described above, are repeated. 


The algebraic process detailed above is straightforward and elementary 
therefore, the details of the algebraic manipulation will not be given. The 
equations of motion obtained from this process for the flap and lag degrees 
of freedom are given in Eqs. 9.10 and 9.11, respectively. 


M F , I * * 7 sn $ , ♦ ( n F , ^ ♦*«,)},« 

** 

B T ( Q) - +' 



+ A fit 
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nT, K, +2M U w L|1 j k, + (RL = 

* * 

§, E «- b 7 [(j,$)-3,^] - ^2Srh ] 

(9.11) 

- — ^'°Pp (sin 0-+ 3>co.S0-) +• 2 Si-wtH^ +• A L ,r 


In the derivation of Eqs. .9.10 and 9.11, the orthogonality conditions 
for rotating beams, Eqs. A. 2 and A. 4 have been used. It should be mentioned 
that Eq. A. 4 should be modified to account for the effect of the total 
coning angle (6 + 3 ) • As shown in Ref. 23, this effect will increase the 

P —2 2 

rotating lag frequency u) _ by an amount 1/2(3 + 3 ) . In the present treat- 

L1 p 

ment this small correction term will be neglected. 


The quantities E^, E^ 2 appearing in Eqs. 9.10 and 9.11 are the elastic 
coupling effects. These quantities are defined in Eqs. L.16 and L.17 (or L.5 
and L.6) . The quantities B appearing in Eqs. 9.10 and 9.11 are generalized 
mass terms defined in Appendix M. 


The quantities A^ 1T and A^ 
and lag, respectively, defined by 


are generalized aerodynamic forces in flap 


A 


FIT 



r 


s 



A 


(9.12) 
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«... - £ 


'UT 


0*1 


V» i 


ft 

L V \ <Ax c 


(9.13) 


Next, the aerodynamic loading terms and will be evaluated. 

The aerodynamic loads and L are given by Eqs. 8.68 and 8.71, while the 


relations for and are given by Eqs. 8.60, 8.61, and 8.62. 
M = 0, and 


For hover. 


9.12, 


Up = Wil t R> 

u T * &SL + ii-(x 0 + €,) 

Substituting Eqs. 9.2 through 9.4 into Eqs. 9.14, 9.15, 8.68, 
and 9.13 and performing the required integrations yields: 


(9.14) 


(9.15) 

8.71, 


A 


FIT 


* nv 

2 R 


— 1 


•=•(»♦*) - \> f * - 


- [ 2 F ”( N §) 


+ (i) 1 i F 


\5 


f (*»,$) J - 2 j F'p P i 9tos$ 



(9.16) 
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and 


A l»t ' Z\rM 


>.[l' (e + s) - L^.J f ^ l 4 


-*lVJ - f + 


(9.17) 


The quantities L 1 , F 1 used in Eqs. 9.16 and 9.17 are defined in 


Eqs . C . 1 and C . 2 . 


In all of the equations associated with the flap and lag degrees of 

freedom, most third-order terms have been neglected. Also, some small 

terms multiplied by C d /a have been neglected. 

o 


9.4 The Feathering Equation of Motion 

9.4.1 The Higher Order Terms in the Feathering Equation 

In treating the feathering equation, it is important to realize that 
the third-order inertia terms in this equation can be of importance and are 
non-negligible. The reason for this is clear when the general form of the 
feathering equation is considered. From Eqs. 8.36 and 8.73, where x = 0, 
the feathering equation can be simply written as 

£ +J! a + §C.oS2&) + If (a)* $ =■ 

c 

J{v- L zr ~ L 1 T (w-p 0 ) + l\ ~ + nr 



(9.18) 
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The left-hand side of this equation is multiplied by the feathering 
moment of inertia of the blade which for the case of = 0 is a small quantity. 
It can be easily shown that after the integrations, all of the quantities on the 
right-hand side of Eq. 9.18 are multiplied by the flapping (or lagging) moment 
of inertia of the blade 1^, or by similar generalized mass quantities which 
are of the same order of magnitude. The ratio of I^/I^ ^ 1000 for most hinge- 
less blades. Thus, the small nonlinear inertia and aerodynamic terms on the 
right-hand side of Eq. 9.19 are multiplied by a large quantity and their in- 
fluence on the stability of the torsional degree of freedom could be con- 
siderable. This has also been found in Ref. 11. Great care is, therefore, 
taken in treating these terms in this derivation and none of the potentially 
important terms are neglected. On the other hand, all fourth-order terms or 
terms of equivalent magnitude will be neglected. 

9.4.2 Final Form of the Feathering Equation 

The final form of the feathering equation can be obtained by substi- 
tuting Eqs- 9.2 through 9.4, 8.68, 8.69, 8.71, 9.14, and 9.5 into Eq. 9.19. 

After performing the required integrations, which are straightforward, and 

2 

dividing the result by , the following final result is obtained. [Also, 

using sin 0=0, cos 0=1] 


i# + = [bVji* 

¥■* 

** * * 

+ P[S’C ^6’ '*(£$) +p r $J * fl T, 


where A , is given by 
T1 


(9.19) 
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An = i t_i(x„[T 4 (e+$)-T 2 \ o - iF'l. | [zl'(o+§)- C t +[eZ z {jtjL ' 3 

* ($F"il + <k L '(t$)]-z irfoie J + i (^)(y t ,-.s)(T‘i-d,fL i p 

- ^fq^e^k.-LK^-L f,KF' 0 - ifila+VL'U.SJk.k* 

- 1^/' + T'fre(«+H)) + {(l-x,)i,kh, + (F'e-x 0 F') 2< § 

- T%p,e + i" (G)k,L".ifrfj,F'} . 

+ k U k < d ’ Hh**) b * £ L 'J}) 


+ $($)Xl y 'l i r[ T, (o+ 2 )->* 7 l ]+iil i />t f:, (e + 2 )-^° F ‘-] 

- Lp F k,[KL'(e+§) +i c £L']+ , }<K[F z (e+t)-KF 1 ] 

+ Cis[T i fif, + ^F , ]j + {l%l[F e (e + §)- Z X 0 F $ J 

+ kl’^d x +t F 3 ]- { e * h ‘idh F '°+id]}^ 

In writting the last equation, the assumption sin ©= 0 , cos0?1.o 

12 3 

has been used. The quantities T , T , T used in Eq. 9.20 are defined by 
Eqs. M.19 through M. 21. 

9.5 Linearization of the Equations for Coupled Flap-Lag-Pitch Motion 
9.5.1 Determination of the Static Equilibrium Condition 
The complicated system of nonlinear differential equations obtained in 
the previous sections will be linearized about the static, linear equilibrium 
condition. Denoting the static equilibrium position in flap, lag, and torsion 
by g° , h° and & , respectively, the dependent variables can be written as 

1 i r _ 
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(9.21) 


' *> , * h ", * ** 

5 * i.*t 

It is important to note that the quantities g~, h^, although similar 
in nature to those given by Eqs. 3. 22, will have different values and will be 
given by different algebraic relations due to the presence of the torsional 
degree of freedom. 


Setting all the dynamic terms equal to zero in Eqs. 9.10, 9.11, 9.16, 
9.17, 9.19, and 9.20, and substituting Eqs. 9.21 in the resulting equations 
gives (after some algebraic manipulation) the following system of equations, 
which for convenience is written in matrix form 


- — 



\ 

r 


( \ 

^11 ^tz $13 


a. 




c - 

S ZI ^ZZ S X3, 

1 

K 

> = < 

L. SN 

“ + 1 

■■ 

C S S 

^31 3Z 33 

1 



T sh 


Cs 

l— 




t J 


J 


(9.22) 


where the quantities S^_. , F^, L gw » T SN , and are given in Eqs. N.l 
through N . 15 . 


The quantities F SN r L SN * T SN represent the static nonlinear quantities 
in the flap, lag, and feathering equations. The linear static equilibrium 
condition about which the equations of motion are linearized will be defined 
as that obtained from solving Eq. 9.22 with the nonlinear terms set equal to 
zero, i.e. , 



(9.23) 
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9.5.2 The Linearized Equations 


Using the static linear equilibrium conditions, determined in the 
previous section, the linearized equations can be obtained by substituting 
Eqs . 9.21 into Eqs. 9.10, 9.11, 9.16, 9.17, 9.19 and 9.20. In this process 
only the linear terms in X-^' X 2 ' and *t> will be retained. The linear steady 
state part of the equations will vanish due to Eqs. 9.23, and the nonlinear 
part of the steady values will be tacked onto the nonlinear equations, which 
will be treated in the following section. After dividing the flap equation by 
and the lag equations by M the linearized flap and lag equations can 
be written in the following convenient manner. 

*, * + X < = F X ** + F *. 

+ * Fy V - F f f 


+ F 


r'f 


(9.24) 



\ =L x,X +L x, 

+ L >, x ' + L y'f + L: f t + L rT 

(9.25) 

The quantities F**, F* , F , F**, F*, F, and L**, L* , L , L**, L*, 

X 2 X 2 X, 4> <P <$> X 1 X A X x <P <J> 

are flutter derivatives associated with the flap and lag equations, re- 
spectively. These quantities are defined in Eqs. N.16 and N.17. The quanti- 
ties g , 9 D 2 represent equivalent damping terms in the flap and lag equations 
given in Eqs. N.16 and N.17 (the last equation in each group). These are 
partially due to structural damping and partially due to aerodynamic damping. 
Finally, it should be mentioned that in this process of linearization, while 
the quantity $ has been replaced everywhere by <t> +<p , the inflow ratio X has 
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been left unchanged and is calculated from Eqs. 4.6 and 4.7. Physically, 
this approximation means that the linearization has been performed in such a 
manner as to allow small variations in the thrust coefficient. 

Next, the feathering equation is linearized. In this case, due to the 
algebraic form of the equation obtained, it is convenient to divide the equa- 
tion by the quantity I (see Eqs. N.18 and N.19) 

E 


I E = I - T^ 


(9.26) 


Physically, 1^ represents the increased moment of feathering inertia of the 
blade due to the blade elastic axis displacement from the feathering axis. 
Thus the feathering equation can be written in a convenient form by 


t ^ T +K st = V*. + T x, *. + T x, 


+ t *A v* 

(9.27) 


where the quantities T**, T* , T , T**, T* , T represent the flutter 

Xi x 2 x 2 X 2 

derivatives for the feathering equation. While g is the equivalent damping 
in feathering, it is of aerodynamic origin only, and is the equivalent 
spring in the feathering degree of freedom. These quantities are defined 
by Eqs. N.24 and N.41. 


9.6 Nonlinear Equations for Coupled Flap-Lag-Pitch Motion 

The nonlinear equations for coupled flap- lag-pitch motion can be 
easily obtained by using the linearized equations of motion obtained in the 
previous section. The nonlinear parts of Eqs. 9.10, 9.11, 9.16, 9.17, 9.19, 
and 9.20 are extracted, the nonlinear steady-state parts from Eq. 9.22 added, 
and the resulting relations added to the linear system, Eqs. 9.24 through 9.27, 
obtained in the previous section. The final form of these equations is given 
below. 
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+ L nl [x x % X) X x % lf v^v|>) 


(9.29) 


- ( Vi, - T y> )%, - (t^ ♦ Ml - T, f )Y - 

TpM -*f,t + T *A * V' + T S»*» + V 1 


T «.( 


(9.30) 


where the quantities F , L , T , T** , T** T** are defined in 

H NL NL NL NX 1 NX 2 , N<f> 

Eqs. 0.1 through 0.6. 


Note that all second derivatives are written on the left-hand side of 
Eqs. 9.27 through 9.29. This form is necessary for the numerical integration 
of these equations. 


145 



SECTION 10 


STABILITY OF FLAP-LAG-PITCH MOTION IN HOVER 

10.1 Introduction 

10.1,1 Brief Review of Past Work 

In this section, the stability of flap- lag-pitch motion will be in- 
vestigated using the equations obtained in the previous section. Both static 
stability or divergence and dynamic stability, or flutter will be considered. 
When searching in the available literature for previous investigations in 
this area, one finds that the complete three degree-of-freedom problem has 
not been investigated before. 

The most extensively investigated case has been the case of flap-pitch 
motion which is similar to the classical flutter problem associated with 
fixed wing aircraft. For this case, both flutter and divergence in hover 
have been investigated by Miller and Ellis in Ref. 11. They concluded that 
the important parameter of the problem is the offset between the blade cross- 
sectional center of gravity and the aerodynamic center. Both flutter and 
divergence boundaries were plotted as a function of this parameter. The im- 
portant effect of preconing on both flutter and divergence has also been 
pointed out in Ref. 11. 

The flutter boundaries in flap pitch for hovering flight were also 
obtained by Daughaday, DuWald and Gates in Ref. 42. Their results were in 
general agreement with those of Ref. 11. Divergence boundaries were not 
presented in Ref. 42. The flutter boundaries obtained support, the claim 
made in Ref. 11, that quasi-steady aerodynamics yields conservative flutter 
boundaries. The experimental results obtained in Ref. 42 seemed to indicate 
better agreement with quasi-steady aerodynamics. Furthermore, it was found 
in this work that flutter could occur when the mass and aerodynamic centers 
were coincident with the elastic axis; this was found to be due to gyroscopic 
coupling between torsion and flapping. 

It is interesting to note, therefore, that no investigation of the 
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static or dynamic stability of the complete flap-lag-pitch motion in hover 
has been undertaken up to the present. In particular, the various effects 
of feathering axis orientation with respect to the hub plane, preconing, 
collective pitch setting combined with elastic coupling, and inplane stiffness 
of tiie rotor on the divergence boundaries has not been investigated. 

Finally, it is worthwhile to mention briefly a few other papers which 
have a connection with the coupled flap-lag-pitch problem. The pitch-lag 
problem in hovering of an articulated blade has been excellently treated by 

Chou in Ref. 43. Chou identified the mechanism of instability as due to 

pitch changes caused by lag motion; the coupling between these degrees of 
freedom was due to rotor head geometry. Stability criteria for hovering 
rotors were obtained. 

The aeroelastic stability of helicopter rotors in hovering flight was 
studied by Zvara in Ref. 44, both theoretically and experimentally. The main 

purpose of this work was to determine the relative merits of the various aero- 

dynamic theories as applied to different rotor configurations. Cantilevered, 
articulated, and teetering blades were evaluated. Only the flap-pitch 
degrees of freedom were considered. 

Perisho (Ref. 45) treated the flap-pitch motion in forward flight, 
including the effect of reversed flow. Blade response curves were obtained 
using numerical integration. Stability boundaries and divergence boundaries 
were not explicitly obtained, although some of these boundaries are indicated 
on his curves. Bielawa (Ref. 35) treated the complete flap- lag-pitch problem, 
considering mainly second-order effects due to distributed torsion, together 
with the effect of time-dependent coefficients. Stability in hovering flight 
and divergence boundaries were not considered by Bielawa. 

Recent flutter analyses seem to be devoted to considering some effects 
while neglecting other effects, usually of equal importance. Representative 
examples of this trend are found in a recent paper by Stammers (Ref. 46), 
where the effect of the periodicity of the coefficients was investigated using 
a regular perturbation method. The Coriolis effects in the flap-pitch equa- 
tions of motion were neglected, reversed flow effects were also neglected. 
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and inertia terms associated with the feathering equation did not include all 
of the important effects. In this paper divergence boundaries were obtained 
as a function of the blade center-of-gravity elastic axis offset, and flutter 
boundaries were also obtained. Another example of this trend is in Ref. 47 
where the effect of the wake in forward flight on the blade stability is 
treated by extending Loewy's aerodynamic theory (Ref. 48) to forward flight 
while neglecting the effect of periodic coefficients and reversed flow at 
advance ratios of 0.3 < M - 0.8. The inertia characteristics of the blade 
in Ref. 47 were represented by taking the tip cross section of the blade as 
a typical cross section of the problem. In this case, stability boundaries 
for forward flight were obtained. Again, only the flap-pitch degrees of 
freedom were considered. 

10.1.2 Objectives of the Present Study 

In the present study, using the linearized equations of motion, the 
divergence of the coupled flap-lag-pitch motion will be investigated first. 

The effect of preconning, elastic coupling, feathering-axis orientation with 
respect to the hub plane, and inplane stiffness of the rotor on the divergence 
boundaries will be shown. 

Next, the flutter, or dynamic stability of the linearized system of 
equations (Eqs. 9.24, 9.25 and 9.27) will be considered and stability 
boundaries similar to those obtained in Section 7 will be given. 

10 . 2 Divergence Boundaries 

10.2.1 The Linearized or Approximate Divergence Boundary 

A linearized approximate divergence boundary can be immediately ob- 
tained from Eq. 9.23 by requiring that det | S^_. | = 0. The divergence bound- 
aries obtained from this requirement are similar to those of Ref. 11. They 

are characterized by the property that the divergence boundary is independent 

o o 

of the values of the quantities g, , h, , and $ which represent the linear 

II o 

static equilibrium condition of the blade. 

From Eqs. N.l through N.9, it is easy to see that the only element of 

the [S . . ] matrix which contains w explicitly is S„. Therefore from the 
ll o ■* 33 
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requirement that det |s^| = 0# the following relation can be written 

■ i{ Mr + f I - r*.? 


( 10 . 1 ) 


where 


t - 


jjajttjji -SnS^aSia 


(10.2) 


_ ^ '3 C ^Al $31 ~ 1 


L 2 = 


V.. ** 


(10.3) 


where S. . is given by Eqs. N.l through N.9. From Eqs. 10.1 through 10.3, the 
approximate or linearized divergence boundary can be easily calculated. 

10.2.2 The Exact Divergence Boundary 

The exact divergence boundaries are obtained by including the effect 
of the static equilibrium condition of the blade (i.e. , g° , h°, and $ q ) . 
Mathematically, the exact condition for the divergence can be shown to be 
given by the following relation: 


it 










(10.4) 

0 
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This requirement is identical to setting = 0 (where v_^ is derived in the 
next section) . 


10.3 The Flutter Boundaries for the Linearized System 


The flutter boundaries can be obtained from the solution of the 
linearized system as represented by Eqs. 9.24, 9.25, and 9.27. The solution 
to this system of equations is given by 


*,* A,e , ' r 
y = , A,« py 


(10.5) 


Substitution of these relations into the linearized equations of motion yields 
the characteristic equation, Eq. 10.6, given on the next page. This equation 
can be expanded to give a sixth order equation which can be written as 


where the coefficients v , 

6 

P.7. 


+■ ^ |> 3 + V +> 9 , |> * ^ * 0 

. . , V are given in Appendix P , 
o 


(10.7) 

Eqs. P.l through 


At the flutter condition 


jo * l W c (10.8) 

Substituting Eq. 10.8 into Eq. 10.7 gives two equations, one for the 
real and one for the imaginary part of Eq. 10.7. 

From the real part 
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(10.9) 


” 0, (o — 0 u) +■ Q 0 - 0 

b C C 1 fc 

From the imaginary part 


^5- Sf ~ 40 * ^3 ■*" H.^1 ~ ° 
or 

^5-H.^ ~ ''jH + ” 0 (10.10) 


The last equation is a quadratic in 0 ) c which can be easily solved; thus 


K)„ 


% i Vi*- ‘li.Oj ' 




From Eqs . 10.11 and 10.9 the flutter boundaries can be obtained. 


( 10 . 11 ) 


10.4 Flap-Pitch Stability Boundary 


In the presentation of the results, it will be useful to have stability 
boundaries for flap-pitch motion similar to those obtained previously for 
coupled flap-lag in Section 4. Therefore, the stability boundary for flap- 
pitch motion will be derived below. From Eq. 10.6 


P+ r + (<»i, + F *i) 

-[//*. r x J 


- [ p y + F rl 


P %i * 


- o 


( 10 . 12 ) 


In analogy to Section 4, at the flutter condition p = iio^ and the 
imaginary part of the characteristic equation (Eq. 10.12) yields the appropri- 
ate relation for the flutter frequency 


COjS Ae ± ft>3 (top, +F 1 S)- Ff Tj£ Ff 
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(10.13) 



The real part of Eq. 10.12 yields 






(10.14) 


The simultaneous solution of Eqs. 10.13 and 10.14, with an appropriate 
relation for the inflow ratio, yields the critical value of the pitch setting 
0^ at which pitch-flap flutter occurs for a given torsional stiffness U)^. 


153 



SECTION 11 


RESULTS AND DISCUSSION 


11.1 Numerical Quantities Used in the Calculation 

Most numerical quantities used in the calculations are identical to 
those given in Section 7.1 unless otherwise stated. For all flap- lag-pitch 
calculations x = 0. 

In the calculation of the divergence boundaries , the following geo- 
metrical blade properties were used: 

Y = 8; a = 0.08; b = 0.0313; I = 0.0013 

These properties are close to those of the Boelkow M105 4-bladed hingeless 
rotor . 

In the calculation of the effect of the torsional degree of freedom on 
the flap-lag- type of instability, the previous values for y and a, used in 
the flap-lag calculations had to be used again. Thus, for this case 

Y = 10; a = 0.05; b = 0.025; I = 0.001 

In the calculation of the effect of lag on the flap-pitch instability, 

the same values of Y, a, b and I were used as for the divergence boundary cal- 

culations. 

For convenience , the inflow for all cases calculated was taken from 
Eq . 7.3.* 

1 1.2 Static Stability Boundaries 

11.2.1 Approximate Divergence Boundaries 

For convenience in numerical calculations, only the linearized or 
approximate divergence boundaries, defined in Subsection 10.2.1, Eq. 10.1, 
were evaluated. Conceptually, these approximate stability boundaries are 
similar to those evaluated by Miller (Ref. 11) . Equation 10.11 includes up 
to the first order, the destabilizing feathering moment due to the drag 

*In all coupled flap- lag-pitch calculations the structural damping in the flap 
and lag degress of freedom * The sructural damping in the 

feathering degree of freedom was always assumed to be zero. 
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acting through the initial deflection (static) in flap together with the 
stabilizing feathering moment due to the lift acting through the initial 
deflection (static) in lag. 

Figure 31 is a typical divergence boundary of this type. The blade 
cross sectional center of gravity and the elastic axis are taken to be coinci- 
dent- 

Thus, x^ represents the aerodynamic center-blade cross sectional 

center-of-gravity offset. As can be seen from Fig. 31, the rotating-lag 

frequency (or chordwise stiffness of the blade) has a considerable effect on 

the divergence boundary because it changes the effective moment arm through 

which the lift produces a stabilizing feathering moment. As can be seen, 

changes in co between 1 and 2 have the greatest effect. It is clear from 
LI 

Fig. 31 that the soft inplane hingeless blade has superior divergence charac- 
teristics. 


Figure 32 shows the effect of collective pitch setting 0 on the approxi- 
mate divergence boundaries. These plots indicate that at lower collective 
pitch setting the blade is statically more unstable. Physically, it expresses 

the fact that the stabilizing moment due to lift is approximately proportional 
2 

to 9 » while the destabilizing moment due to drag is proportional to a power 
of 0 somewhere between 1.5 and 1.2. Calculations made for pitch flap with the 
exact divergence boundary (nonlinear) show that increasing G could be de- 
stabilizing . 


Figure 33 shows the destabilizing effect of the preconing, which has 
also been indicated in Ref. 11. From this plot it is clear that preconing 
has a strongly destabilizing effect on static stability. 


spect 
Again , 


Figure 34 shows the effect of feathering axis orientation $ with re- 
to the hub plane, measured in a plane perpendicular to the hub plane, 
this is a destabilizing effect, similar to 3 , but not as strong. 


The effect of elastic coupling is stabilizing. Figure 33 was recalcu- 
lated, including elastic coupling; it was found to be stabilizing, but only to 
a minor degree. The effect of elastic coupling is dependent upon the combi- 
nation of flap and lag frequencies. Therefore, for some other combination of 
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these quantities, the effect could be stronger. 


11.2.2 The Exact Diverge nce B ound ary f or F l ap Pitch 

In order to illustrate the considerable difference between the approxi- 
mate and exact divergence boundary, these two are given in Fig. 35. The exact 
divergence boundary for flap pitch is treated in Appendix R, and the appropri- 
ate curve was calculated using Eq. R.5. 


As can be seen from Fig. 35, the exact divergence boundary is approxi- 
mately 33% higher than the approximate linear divergence boundary. Additional 
cases for different values of collective pitch setting 0 were also computed 
and the same difference between approximate and exact boundaries was observed. 

o 

It is of interest to compare the initial static values of g, ,'i> evalu- 

. 1 o 

ated at the exact divergence boundary and compare them with the values result- 
ing from Eq. 9.23. 


For the x 

A 


0. 04 point of Fig. 35 


($> ) 

o exact 


<£) . . 

° linear 


0.04919 

0.0623 


* g i^ exact 
‘^l 1 linear 


o.ioii 

0.1117 


The flutter calculations are dependent upon the initial values g^,$ Q ; 
thus, the exact calculation of these quantities, can have an effect on the re- 
sult of the flutter calculation. 


1 1. 3 Flutter Boundaries 

11.3.1 General 

The purpose of this section is twofold. First, it will show the 
effect of the addition of the torsional degree of freedom on the flap-lag- 
type of instability treated in the first part of this work. Second, it will 
show the effect of the addition of the lag degree of freedom on the flap-pitch 
type of instability as evaluated from Eqs. 10.13 and 10.14. 

The coupled flap- lag-pitch stability boundaries are evaluated from 
Eqs. 10.9 and 10.11 using a computer program which searches and iterates for 
the flutter points occurring between 0 < 0 < 0.5 for a given combination of 
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flap, lag, and torsional frequencies. This program also continuously checks 
the value of v , evaluated using the linear static equilibrium condition, for 
the occurrence of divergence. 

The coupled flap-pitch stability boundary is evaluated using Eqs. 10.13 
and 10.14 by a separate program which searches and iterates for the value of 
at which flutter occurs at a given collective pitch setting 0. This pro- 
gram evaluates only the upper branch of the stability boundary. 

11.3.2 Effect of Torsional Degre e of Freedom on the 
Flap-Lag-Type of Instability 

* 

In order to xllustrate this effect, the following cases are considered: 


Case 


(a) 

U) F1 = 1.175; 

II 

r-H 

►4 

3 

1.075764 



(b) 

W pl = 1.175; 

(Jj = 

LI 

1.28303 

Case 

2 

(a) 

W F 1 = 25 ' 
= 1.25; 

ii ii 

13 1 3 

1.11966 



(b) 

1.39403 


Note that all cases plotted are calculated without the effect of elastic 
coupling. 

Cases 1 and 2 correspond to the appropriate points on the stability 

boundary given in Fig. 8 for 0 = .20. Point (a) is on the lower branch 

c 

while (b) is always on the upper branch. 

Figure 36 shows the effect of torsion on Case 1(a) . At a high 

value of 10 (li) = 100) , this simulates a torsionally rigid blade; for this 
o o , 

case, 0^ = 0.2068 which compares well with the value 0^ = 0.20 obtained from 

the flap-lag calculation. As the torsional stiffness is gradually decreased, 

— 't 

the value of 0 is slowly increasing. In the vicinity of to = 32, 0 in- 
c o c 

creases sharply. At some very high value of 0 divergence occurs. Thus, the 
effect of addition of the torsional degree of freedom is stabilizing for the 
lower branch of the flap-lag stability boundary shown in Fig. 8. 

When the same process is repeated for the upper branch of the flap-lag 
stability boundary. Case 1(b), it is seen from Fig. 37 that the value of 0^ 

* 

A limited number of additional cases will be considered in Ref. 9. 
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is continuously decreasing and at W q = 5, it reaches a value of 0^ = 0.047. 
Thus, the addition of the torsional degree of freedom is strongly de- 
stabilizing for the upper branch of the flap-lag stability boundary. 

During this process, the flutter frequency remains fixed for both the 
upper and the lower branch of the stability boundary, and is approximately 
equal to the lag frequency. 

For the sake of completeness , two additional points on the flap- lag 
stability boundary were considered. These are Case 2(a) and Case 2(b); the 
results are shown in Figs. 38 and 39. Again, the torsional degree of freedom 
is destabilizing for the upper branch and stabilizing for the lower branch. 

11.3.3 Effect of Lag Degree of Freedom on the 
Flap-Pitch-Type of Instability 

In order to evaluate this effect, a pitch-flap stability boundary is 
required; this boundary is shown in Fig. 40. As can be seen, the pitch 
setting is strongly destabilizing. 

In Fig. 41, the flutter frequencies associated with the stability 
boundary shown in Fig. 40 are given. It is of interest to note that the 
flutter frequencies in pitch-flap are usually lower than cu^. 

A representative point on the flutter boundary (Fig. 40) with 

6 = 0.20 w = 7.8125 uL, = 1.2 

c o FI 

will be used to investigate the effect of adding the lag degree of freedom. 

The results of adding the lag degree of freedom are shown in Fig. 42. 

For 0) =20, which represents a blade with the lag degree of freedom 

LJ. 

effectively suppressed, 0 = 0.197, which is quite close to 0 =0.20. De- 

c c 

creasing the lag stiffness of the blade gradually reduces the value of 0 . 

c 

The flutter frequency during this process is shown in Fig. 43. It is in- 
teresting to note that during this gradual decrease, the flutter frequency is 
essentially the flap-pitch flutter frequency. 

In the vicinity of to = 5.0, where the rapid dip in the value of 0 

LJ. C 

occurs, the nature of the instability changes; it becomes a flap-lag-type of 
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instability and the flutter frequency is very close to the lag frequency. 

— 'Xi 

At u) = 2.25, 0 attains a minimal value of 0 =0.04. This region is 

LI c c 

obviously identical to the region of low 0 's shown in Figs. 37 and 39. 

- n, c 

In the vicinity of u)_ = 1.1, a rapid growth in the value of 0 occurs, and 

LI c 

no flutter occurs below oj £ 1.0. This seems to indicate that the soft in- 
plane hingeless blade seems to have the best coupled flap-lag-pitch flutter 
characteristics. This apparent advantage of the soft inplane hingeless blade 
should be considered within the limitation of the analysis performed in the 
present study. The various other instabilities associated with this type of 
blade such as ground resonance, air resonance , and instabilities due to positive 
preconing, were not considered and could be major design problems in the con- 
struction of a stable hingeless rotor system. 

For a few cases, the effect of elastic coupling on these stability 

boundaries was investigated. It seems that the elastic coupling is sufficient 

to eliminate most of the unstable regions, except the unstable region in the 

vicinity of u) = 2.5 shown in Fig. 42. 

LI 
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SECTION 12 


SUMMARY OF RESULTS AND CONCLUSIONS 

12.1 Conclusions for Coupled Flap- Lag Motion 

Conclusions for Hovering Flight 

1. In the vicinity of the critical condition (0 0^) , when flutter 

occurs, the amplitude of flap motion is much smaller than that 
for lag. This is due to the fact that the lag mode is the poten- 
tially unstable mode, due to its low damping. In this region the 
ratio of lag over flap amplitudes can be between 5 to 10. 

2. The structural damping in flap has no effect on the stability of 
the system. On the other hand, a small amount of structural 
damping in lag (n. T = -5% of critical damping) is sufficient to 
stabilize the potentially unstable lag mode by raising the values 
of the critical collective pitch setting 0^ above values which 
could occur in practice. According to Ref. 3, the elastic coupling 
effect is also sufficient to eliminate most unstable areas inside 
the flap-lag stability boundary. In an actual hingeless rotor, these 
two effects will always coexist, indicating that the flap- lag sta- 
bility problem may be readily avoided in an actual rotor design. 

3. The limit-cycle-amplitude-response curves in lag are steep and the 
limit-cycle amplitudes are large. This means that the nonlinearities 
in the system are weak and they cannot stabilize (in a practical 
sense) the response once the linear-stability boundary is exceeded. 

4. Combinations of flap and lag frequencies, corresponding to the 
region of unstable limit cycles in the stability boundaries (Figs. 3 
and 4) , should be avoided in the design of hingeless helicopter 
blades, because for this case a blade stable from a linear point of 
view can become unstable, if it encounters a disturbance which is 
large enough. 



5. The structural damping in lag tends to reduce the limit-cycle 
amplitudes in lag and it also reduces the slope 0x2/39) of the 
limit-cycle-amplitude-response curves in lag (see Figs. 14 and 
15) . 


6. In the neighborhood of the critical condition (0 > 0 ) , the 

c 

value of K , defined in Eq. G.71, determines the behavior of 
JR 

the blade motion in the large amplitude range because the sign 

of 3|A o |/3ij^ is the same as for ic^ R when |a q | is sufficiently 

large. Therefore, a large positive K represents a strongly 

destabilizing nonlinear effect. It is undesirable to have a 

blade with a positive < because in this case instability can 

3R 

always be excited, if the disturbance is sufficiently large. 

In violent maneuvers or in strong gusts, such large distur- 
bances may occur in practice. The value of K is independent 

3R 

of the advance ratio and is a complicated function of • 


(A) .v 

L10' r 


7. 


etc. 

In hovering, for 0 > 0 and K. 


c 3R 

for lag motion is approximately given by 

\'k 


< 0, the limit-cycle amplitude 


1 

l /Cw J 


kiR 

while the limit-cycle amplitude in flap is given by 


2 . 1*1 



ft 


8. Figure 10 shows that the effect of using the approximate 
centrally-hinged, spring-restrained rigid blade mode of 
the elastic hingeless blade is to significantly shift the 
region of flap-lag instability to lower values of blade 
frequency. This indicates that the use of the exact ro- 
tating mode shape could be of importance in the calcula- 
tion of flap-lag stability boundaries. 

9. The flutter frequency of flap-lag oscillations is very close 

to the lag frequency u) . Therefore, for practical purposes, 

LI 

the lag frequency can be taken as equivalent to the flutter 
frequency. 
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Conclusions for the Forward-Flight Cases 

10. When the rotating lag frequency of the blade is 
not in the vicinity of the 1/2 or 1, the t im e de- 
pendent coefficients (or forward flight, y) in the 
equations of motion has only a forcing function type 
of effect. Thus, the time dependent coefficients 
will not have any influence on the stability of the 
system for this case and for the range of advance 
ratios considered in this study (0 < y < .3). 

11. Below the critical condition (0 < 0 ) for all cases, 

c 

the flap amplitude X 1 * in forward flight, is larger 

than the lag amplitude X 2 * At approximately 0.10^ 

below critical the lag amplitude starts to grow, and 

for 0 > 0 , the lag amplitudes are much larger than 
c 

the flap amplitudes . In order to eliminate large oscil- 
lations in lag, it is reasonable to design the blades in 
such a manner that no angle of collective pitch inside 
the flight envelope should ever exceed 0.70^. 

12. Due to the large amplitudes of vibrations in lag which can 
occur in both hover and forward flight, the fatigue life 

of the blade for inplane vibration could possibly be an im- 
portant criterion for hingeless blade design. 

I 

1/2 

13. For the general forward flight case. Case B (y = e , y =1, 

o 

U) / 1/2 + ev or 0 ) / 1 + EV) , the blade response and stability 

c c 

2 2 

in the neighborhood of 0 is determined by y + K_ 1(0 - 0 )/y , 

C O OK 2R C 

2 

(A - A )/y ] and ic_ . When the first of these quantities is 
o oc 3R 

positive and (C^ R < 0, stable limit-cycle oscillations occur. 

The flutter mode in this case has the functional form of 
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14. Case C (y = £, U, = 1 , y = 0, w = 1/2 + £V) . In this case, 

1 o c 

the parametric excitation is predominant. Up to the first 
order, the amplitude of the lag motion, near the critical con- 
dition, can be approximated by 



Near the origin (A =0), the stability is determined by Re(s,) 

o 1 

(s^ is given by Eq. 5.79a on page 46). For Refs^) > 0, the 
system is unstable and for Re(s^) < 0, the system is stable pro- 
vided k < 0. 

jK 

1/2 

15. Case D (y *» e , y = 1 , y, = 0,u> = 1 + ev) . The parametric 

ole 

excitation has an effect but forcing is predominant because the 
system is excited very close to its resonant frequency. Up to • 
the first order terms, the lag motion can be approximated by 



The stability near the origin (A =0) depends upon the sign of 

o 

Re(s 1 >. Where s^ is given by Eq. 5.104, the system is stable for 
Refs^) < 0 and is unstable when Rets^) > 0. For y > 0.25, this 
case is usually characterized by large amplitudes of blade re- 
sponse. Therefore, the lag frequency range of 0.88 < w <1.12 
should be avoided in the design of hingeless helicopter blades. 

16. The important nonlinear terms in Eqs. 3.23 and. 3.24 have been 

identified. The strong destabilizing term in Eq. 3.23 is 
13*2 

C^ (£/R) 0F y^ 2 ' When neglecting this term, the regions of un- 
stable limit cycles in Figs. 3, 4, and 8 vanish completely. 

„ 18*2 

The term - C^(J6/R)L ln Eq. 3.24 is stabilizing. When neglecting 
this term the regions of unstable limit cycle grow and cover 
approximately 60% of the stability boundary. Thus, it is danger- 
ous to neglect nonlinear terms when investigating the nonlinear 
blade stability problem. 


163 



17. p or a n forward-flight cases, the stability and the amplitudes 
of blade response are considerably affected by the requirement 
of trimmed flight at fixed C^. Flutter and blade response cal- 
culations, which neglect this effect when evaluating the effect 
of forward flight, have a doubtful validity. 

12.2 Conclusions for Coupled Flap-Lag-Pitch Motion 

1. From the numerical calculations performed it was 
found that the results associated with the 
stability of the coupled flap-lag-pitch system 
are quite sensitive to the numerical values of the 
coefficients in the feathering equation. The 
various quantities F* 12 } > etc * 

are dependent on the assumed mode shape, consequently 
the use of the exact mode shape of the blade (corresponding 
to the actual load distribution) may be important for certain 
flight conditions. Similarly the use of the exact, span- 
wise varying, inflow could also effect the results. These 
effects, can also combine to modify the static equilibrium 
position of the blade, thus affecting both the divergence 
and flutter boundaries. Use of an additional elastic mode 
in the flap degree of freedom could also be important. 
Therefore, the results and conclusions, presented below, 
should be considered as subject to certain limitations due 
to the fact that the various effects mentioned above have 
not been included in the analysis. 


164 



From the approximate static stability boundaries given, it is 
clear that the aerodynamic center-blade elastic axis (c.g. 
in reality) offset is a destabilizing effect. Similarly the 
effect of preconing (3^) and blade feathering axis orienta- 
tion with respect to the hub plane ( 3) are also destabilizing. 
The strongest destabilizing effect seems to be the preconing. 


The lag degree of freedom has an important effect on static 
stability because the deflections in lag generate an effec- 
tive moment arm through which the lift produces a stabiliz- 
ing feathering moment. Thug, a soft inplane hingeless 
blade has much better divergence characteristics than a blade 
which is stiff in the chordwise direction. 


Conceptually, the approximate linear divergence boundary is 
not the correct divergence boundary of the system. For the 
small number of cases considered, the approximate divergence 
boundary seems to be nonconservative. 

The addition of the torsional degree of freedom has a 
stabilizing effect on the lower branch of the flap-lag 
stability boundaries (Figs. 3,4, and 8), while it is 
destabilizing for the upper branch. Thus, the upper 
branch of tiie flap- lag stability boundary should be 
avoided in the actual design of rotor blades. 

The addition of the lag degree of freedom with to > 1.2 

LI 

is destabilizing for the coupled flap-pitch motion. 



With rotating lag frequencies of cj t1 < 1.1, the coupled 
flap- lag-pitch motion is quite stable. Thus, the soft 
inplane rotor with a rotating lag frequency in the 
vicinity of a) T1 = 0.75 (midway between .5 and 1.0) seems to be 
the best from both the divergence and flutter point of view. 
This apparent advantage of the soft inplane hingeless blade 
should be considered within the limitation of the analysis per' 
formed in the present study. The various other instabilities 
typical of this configuration were not considered (see 
page 159) . 

Due to the strong effect of the initial deflections 
g°, h°, on both the static and dynamic stability 
boundaries of the blade, flutter-divergence interactions 
could easily occur. 

When including the third order inertia terms in the feather- 
ing equation it has been found that it is also necessary 
to include a considerable dumber of third order terms , as- 
sociated with the aerodynamic loads (because they are of the 
same order of magnitude) . This is required in order to rep- 
resent correctly the behaviour of the blade. It should be 
pointed out, that this conclusion may be modified by the 
considerations given in the first conclusion of this section. 

Based upon this work, a parametric investigation aimed 
at determining the optimum configuration for a hingeless 
blade seems feasible. 
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1 


■ficients Defined in Appendix C 


Lag Numerical Value 

Coefficient (Nondimensional ) 


1 

It 

0.2888881 

L 2 

0.3999991 

3 

h 

0.0 

It 4 

0.2253962 

L 5 

0.0 

It 6 

0.0 

7 

It 

0.2059954 

L 8 

0.2567891 

I* 9 

0.0 

D 10 

0.3716033 

I* 11 

0.4999981 

12 

1. 

0.0 

L 13 

0.2567891 

L 14 

0.2059954 

15 

F 

0.0 

F 16 

0.3333315 

F 17 

0.3333313 

F 18 

. 0.1901494 

F 19 

0.1901494 


Mass Quantities Defined in App. B 
M fi = .7703672 
= .7703672 
C = 1.199996 


TABLE 2 


NUMERICAL VALUES OF B 1 , T 1 COEFFICIENTS DEFINED IN APPENDIX M 


Coefficient 


I 1 

B 2 

B 3 

1 4 

1 5 

1 6 

B 7 

5* 

I 9 

5 10 

1 11 

1 12 

B 13 

B 14 


B 


■15 

-16 


B 17 

B 18 


Numerical Value 

.8666655 

.9037014 

.7703676 

.7703676 

.8666654 

.9037012 

.7703671 

.7703676 

.9037012 

.8666654 

.8666642 

.8666642 

1.1999960 

1.5000000 

.9037001 

.9037001 

.8148102 

.8148102 

.333333 

.500000 

.250000 
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FIG. 11 EFFECT OF VISCOUS STRUCTURAL DAMPING ON ANGLE OF COLLECTIVE 
PITCH AT CRITICAL CONDITION 
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FIG. 16 TYPICAL TRIM CURVES FOR C T - 0.005 
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FIG. 17 TYPICAL TRIM CURVES FOR C T ** 0.0075 
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^l^av' *2 Note : x^» X 2 are fia P ar } d lag deflections, about the 



FIG. 19 AMPLITUDE RESPONSE CURVES IN FORWARD FLIGHT, COMPARISON 
OF RESULTS FROM PERTURBATION AND NUMERICAL INTEGRATION 
AT SUBCRITICAL CONDITIONS . 
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FIG. 22 EFFECT OF C ON AMPLITUDE RESPONSE IN FORWARD FLIGHT 
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FIG. 23 AMPLITUDE RESPONSE CURVES IN FORWARD FLIGHT, POSTCRITICAL 
REGION (FROM PERTURBATION METHOD) 
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CASE B 



1 1 1 1 1 ® (trim) 

0.243 0.236 0.21 0.204 0.237 

V 

FIG. 25 EFFECT OF FORWARD FLIGHT, WITH TRIM AT FIXED C , 

T' 

ON p. (CALCULATED FROM EQ. 5.71) 

• C • 


198 


£j°* e ; » X 2 are flap and lag deflections, about the equilibrium 

position, at blade tip nondimensionalized with respect to the 
blade length. 



_| 1 1 ( 1 1 +.0 

0.248 0.241 0.215 0.205 0.207 0.219 

FIG. 26 AMPLITUDE RESPONSE CURVES, IN FORWARD FLIGHT, FOR CASE D, 
IN POSTCRITICAL REGION (FROM NUMERICAL INTEGRATION) 
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REAR VIEW 



FIG. 27a DISPLACEMENT FIELD WITHOUT 
PRECONING (8 =0) 



[ASSUMPTION (4) , 
SUBSECTION 8.2] 


DEFORMED BLADE 



TORSION AND WITHOUT 








POSITION OF BLADE CROSS SECTION 



BLADE CROSS SECTION BEFORE DEFORMATION 


FIG. 28 BLADE MODEL AND POSITIONS OF THE CROSS SECTION 
BEFORE AND AFTER THE DEFORMATION 
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Elastic Axis Blade c.g. Coincident 

10 



FIG. 33 EFFECT OF PRECONING B AND PITCH SETTING 0 ON THE 

P 

APPROXIMATE DIVERGENCE BOUNDARIES 
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Clastic Axis and Blade c.cj. coincident 
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FIG. 36 EFFECT OF TORSIONAL DEGREE OF FREEDOM ON FLAP-LAG-TYPE OF 
INSTABILITY 
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FIG. 38 EFFECT OF TORSIONAL DEGREE OF FREEDOM ON FLAP-LAG-TYPE 
OF INSTABILITY 
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TYPE OF INSTABILITY 
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FIG. 41 FLUTTER FREQUENCY AS A FUNCTION OF PITCH SETTING 



FIG. 42 EFFECT OF LAG DEGREE OF FREEDOM ON FLAP-PITCH-TYPE OF INSTABILITY 
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APPENDIX A 


ORTHOGONALITY CONDITIONS 


The orthogonality condition for a rotating beam vibrating out of the 
plane of rotation can be obtained following Bisplinghoff (Ref. 12) 


I 

JlH 


i i 

i dZ ) d% -° 

*9 


(A.l) 


when i / k and 

' j‘hi In ") ) J *° = 


Jl z t 


2 // ■ ■ j _ .3 

= I W*) olXo l 


(A. 2) 




JV- 


when i = k. 

For a beam vibrating in its plane of rotation, the orthogonality con- 
ditions are (Ref. 23) 

jyih-iC £ J '°+ e J ( 1- d *‘) 


JlH 


o 


clXo - 0 


n. 


when m f n and 
/ 


— ( ( E $ Z ^ Y*) **<> + ^J(\ nX i <**' ) ( tc) dX* ~ 

J o o Jj( e 

- jt 3 J my* dxo - ( o->Li j £* J myfclxo 


(A. 3) 


when m = n = i 


(A. 4) 


These relations also determine to . and u 

Fx Lx 


21cS 



APPENDIX B 

COEFFICIENTS ASSOCIATED WITH GENERALIZED MASSES AND DAMPING 

h' 1'f+VK. 

fo- = t 3 J o % i‘ k ( ) clXo /ik 

Mfc= 

Hu* l 1 ]'**?? +*•/%, 

C * ly >*r c <**•/ h 

l mr = ej'rltl ( C "c •>*) **A 

''O 

fa -)i« % n'e d *’) f i “4A 


The nondimensionalized viscous-type structural damping coefficients are 

* - *** ** V 

^ mXo cfXo 

q __ / ^v«~» - = toJL.-tiL-'l'". 

Q SC If 


Jlj mXo L c/Xa 


where n . and n . represent the fraction of critical damping associated with 
SFi SLi 

the ith mode in flap or lag. 
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APPENDIX C 


EXPRESSIONS ASSOCIATED WITH AERODYNAMIC LOADS IN FLAP AND LAG 


Note that all integrals are performed between the lower limit A and 
the upper limit B. A and B represent the tip loss factors. 


F- = J 

F = fir, dx ° 

Fc~ x ° 

Lc-fr^ 

Fh jq. dxo 

L l , =J‘ <jc*i)rc d * ; 


L Fj* z rc d *° 

frjyx-K* 

(c ' 1) L%/xj(5)r^* 

ti-fttcik dr ° 





L* * = J% p d *° 

ti-J dZ 

F** /?*) #i-i < d *’ 

*£'{* vj* 1 

Ci *f*% t c dZ 

*/ W* 

i« d *° 


CC.2) 
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Ck’ji&Wypxi 

C'llLYr^' 

L«fnr n dr ° 

&Uxc dF ° 

L't 

FL’jwJ*° 

La 

"n 

£ 

ii 

L'L^Jny^Fdu 

F n -jxJ*, 

L-‘U‘jK% 7 U d *° 

f 1 - J X* dx° 
F n - ['-* 

Jo 

F Zo = J * ('*[)) 

LcufrAp/*' 

Lcu=Jf. i k yj*° 


The quantities A^, 


A . are defined by the following relations : 


A 


F-C 



Fc Q-F?\o+F? h- coif + /i. (z e F c - 

_\o Fc) S£nf-/±F'?6 cos If - \ic J± s in Z<fJ 
X/ ^ 

- •/*• f f 4 Pel s tn 2 +) & - £ ( P,l + A f a 
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* 

Srn 


” [( X f~Cm & ~ ^°) ~ ^ lc c OS^ 4-XjA,9 F Crn S Cn'f'J L 

+ [(k f 6 Cr^ k r + $)/*£'*+%£. +({j F Zi k *]J 


(C. 3) 


^ u ~ fe- l£j j £ \o( L< & - L?)io) -i Lc hie +Qo ( Lc j + 


( L 5 c 9\,c - x \o >/c L-l )tosf 9 Lc +2Cjv juLc ) SCnf 

~ it ^ 10 ^ ^ Ldo J jlLc ) t'OSZ'f' 4 - &\,c L c ScnX.'f'J 

^ fa (r)[( ^ ~ c ^ L^°) - x\tc Lck coS 'jf' + /*Lci& Sv ” K f'J 

+ [^.(e Lck°-%Lck \o)cQS>f-J^ \ic Lc% lZ ('l-t- CoS w)* 

+ & j££k sc»i+]- c @)[(l£ \J + lCg Let ) + 

/5* ^ ih 

'hie 9 Lei */* -b ^5-? Jin. Lek ^ tri fl- If// 4 ’ (~ J ^ ^ #rI fa 

+ * i 11)™* + (±)*( Lm y t +e LctZ £ C )]J 


(C. 4) 
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APPENDIX D 


INFLOW RELATIONS 


D. 1 Inflow Relation in Hovering Flight 

(Ref. 19) gives the following relation 


(D.l) 

Thus, constant induced velocity over the whole rotor area is compatible, 
with a twist variation inversely proportional to x. In the present report only 
constant values of collective pitch will be considered. Therefore, an ex- 
pression for the constant inflow will be derived in such a manner that it repre- 
sents a weighted integral of the varying inflow, resulting in the correct value . 
for the thrust coefficient. 

From blade element theory with constant induced velocity 


Extended blade element theory 
for the induced velocity 


_£L = r/77 

dlR 4(> [V 


7>l 6 X 
CL<J 



(D.2) 


Also, from blade element theory 


otT- J 


olX 


(D.3) 


From Eqs. D.l and D.3 



(D.4) 


Evaluating the integral in Eq. D.4 yields: 
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(D.5) 


‘~P 


From Eqs. D.2 and D.5 


where 


er 


U-- ne 

CTcb 

The advantage of relation D.6 when compared to (D.7) 


‘-J 


(D.6) 

/ 


(D.7) 


is due to the fact that, for a given C T , Eq. D.6 agrees to within 
't 4% with A q = / C^/2 , while Eq. D.7 gives a discrepancy of 20%. 

D.2 Inflow Relation in Forward Flight 

In forward flight, the constant part of the inflow can be obtained from 
the well-known relation given in Ref. 19 



"tan oil 


-b 



(D.8) 


For a given value of C T and p, A q can be obtained from Eq. D.8 by 
solving a quartic equation. This has been done by using a simple iterative 
process on a digital computer. 

In practice it has been found that it is possible to replace Eq. D.8 
by a simpler relation which provides a reasonable approximation for A q in the 
range of advance ratios (0 < p < 0.3) and thrust coefficients (0.005 < C T < 0*0i5, 
for flutter) considered. 

<\, 'V, 

For low values of p, p o> 0, For the value of p 't 0.3, 

it has been shown in Ref. 24 that 
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i Or\, 0 


■i Ct 
ZjU, 


Therefore, an approximate relation which would be good at the endpoints can 
be written as 




z[^i 


ctu.ot.fi 4 


Cr 

<Z ( 


w 


Using this expression, A q can be approximated by 


(D.9) 


\*U.icvn ot* f Cr (d .10) 

zfVVJ 

Expression D.10 can be considered to be an approximate solution to the 
quartic (D.8) . The maximum error due to the approximation occurs at u = 0.05 
and is 7%, at y = 0.1 the error is 2.5%, at p = 0.15 the error is less than 
0.8% and afterward (D.10) is equal to (D.8) up to 3 significant digits. 


Therefore, Eq. 8.10 was used instead of solving Eq. D.8 by iteration 
for every case considered. 

Note that A in Eq. D.9 is calculated from A = /c /2 . 

H H X 
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APPENDIX E 


THRUST COEFFICIENT FOR A HINGELESS BLADE IN FORWARD FLIGHT 

In this appendix, the thrust coefficient of a hingeless rotor will be 
derived. All quantities are referred to the hub plane . According to 
Bramwell (Ref. 20) , the motion of the blade as represented by one elastic 
mode is sufficient for the derivation of the thrust coefficient. The contri- 
butions due to lagging velocity in U p and will be neglected. 

With these assumptions and Eqs. 3.7, 3.8, 3.14, and 3.15, the velocity 
components U p and U^ can be written as 

+ ]** (£• 1 ) 


U T = dip (E.2) 

From Eq. 3.12 

T=m>_ [ 7 L 2 t<J*olf= 2^(^)apl / / U?(6-<Jp)dX W 

&TT Jo J A ITT Jo J A Ur) (E. 3 ) 

Assuming a cyclic pitch variation given by 

Q-Go-tO^^Cnf + eic^f (E.4) 


the corresponding flapping motion can be represented by 

(E - 5 > 

The substitution of relations E.l, E.2, E.4, E.5 into E.3 yields, 
after a considerable amount of algebraic manipulation. 
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— 'V — % 

For A = 0, B = 1.0, by integration by parts 

<-r s 

3 ^ 

From its definition F ^ 0.5, therefore, the last term of Eq. E.6 
negligible, and the thrust coefficient is given by 

C T - qa-tl\l6. (3F*+i F' ? K+JkF' 7 , 

jL, ( R, / ( 3 2 

where the quantities F 16 , F^, F'*" 8 are defined in Appendix C. 

Equation E.7 agrees with the corresponding equation given in Ref. 



(E.6) 


is 


(E.7) 


20 . 


227 



APPENDIX F 


APPROXIMATE TRIM EQUATIONS FOR A HINGELESS ROTOR 

F.l Assumptions and Basic Equations 

In this appendix, an approximate method for calculating the value of 

collective pitch 0 q and the appropriate values of 0t R , 0 , 0^ s for a given 

flight condition will be described. The flight condition is determined from 

the values of C and u which are assumed to be known. 

T 

The method of calculation is based upon the following assumptions: 

(1) The helicopter is in steady level flight. 

(2) Pitching and rolling moments on the rotor are 
equal to zero. 

(3) The rotor hub and the helicopter center of gravity 
coincide, i.e., the whole aircraft is represented 
by a point mass coinciding with the hub. 

(4) The cyclic pitch variation is assumed to be given 
by Eq. E.4 and the corresponding flapping motion 
is given by Eq. E.5. 

•\» 

(5) The angle a is small. Then sin a = a , cos a = 1. 

R HR R 

The geometry of the problem is given in Fig. F.l. 

The equations of equilibrium, tangential and normal to the flight path, 
can be written as 

m T " ^tot (f>1) 


<F ‘ 

From assumption (2) , the moment equilibrium in pitch and roll yields 


Mmr- 0 

(F . 3) 

Mtr-O 

(F.4) 
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These four equations are sufficient to determine the four unknowns 

0 , a . 0, , and 0, . 
o R lc Is 

From assumption (1) , y = y = 0 V = 0, and Eqs. F.l and F.2 reduce to 
^ToT^° (F.5) 




(F.6) 


From Fig. F.l 


DtoT = -TstnOi R + H CoSoifi-h Df> 


(F.7) 


L == Teas oi R 


(F.8) 


The drag of the helicopter can be represented by using the equivalent 
flat-plate area 



(F.9) 


With these relations and assumption (5), Eq. F.5 can be rewritten in coef- 
ficient form 

- Cr°<-« + C h + (F.10) 

In order to solve the problem, the quantities C , M , M„ must be 

H mr xr 

evaluated. 


F.2 The Pitching and Rolling Moments 

Neglecting the Coriolis forces, Eqs. 3.1 and 3.3 yield 


— — [( E l)y d 1 — — ( ) -+ m, 

DXo' l o -fixS J 'iXe/ 



(F.ll) 
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By integration, the moment at the blade root due to aerodynamic, centrifugal 

and inertia forces can be written as (for Z = R, x = x, and e =0) 

o 1 


Mr. (o,+)=J jT L 2 -nuJL L ^ L ~>rtdl l ^J XodXo 

Using one elastic mode given by Eq. 3.7, Eqs. F.ll and F.12 yield 

M * ( °* ] = i l£[ (£I \ S H (t & )- 4 *° 


(F.12) 


(F .13) 


From the free vibration problem equation (A. 2) 


c/x 


h 7% y & ( T %) = no **'*> 


(F. 14) 


From Eqs. F.13, F.14, and D.15 

/ 

M R lo,+)= 7/ df, (F 15) 

According to Fig. F.2 below, the pitching and rolling moments are ob- 
tained from M by a simple vector decomposition. 

R 



FIG. F.2 BLADE ROOT MOMENT DECOMPOSITION 
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(F.16) 


Mm= Mk (0 7 t) c °s+ 
Mt - W) 


For trim purposes, the average values of these quantities, per revo- 
lution, are required. 

IT 

lefO.'P \C0b4' Ot'f' 

(F .17) 

2/7 «2o 

From.Eqs. F.17 and F.15 


M , ma. - J L / M^o,* ) cobf df 
2J Jo 

Ml*.* -1 / Mr(. 0,+) SCnfdf 

2/T 


He** ± («r F *-4)l s ji\ £ nZ^dx. 


(F.18) 


Equation F.18 represents pitching and rolling moments due to one blade, 
for blades 


M m r = Mm* 

Mi r = Me* n* 


(F.19) 


Thus, Eqs. F.3 and F.4 are equivalent to the requirement that 


tc " U ° 

F.3 Flapping Coefficients 


(F.20) 


The so-called "flapping coefficients" are obtained from the steady-state 
solution of the flap equation when the Coriolis forces, damping and displace- 
ments or velocities, due to lag motion, are neglected. In the analogy to 
previous derivations, only one elastic mode is used. From Eq. 3.18 
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— ** 


Mf 4 f t ■+ Mf/ ff 4 - A Ft 


(F.21) 


With = 0, from Eq. C.3 


/>c,= %(${[ f'9-F ! X*+ 

_ J?_ F i e col IF -u ( F l ’coiF + 

2 7 

M F f i(n 2 F ) % - | ( F *+jh F * u '» t) l J 


(F .22) 


Substituting Eqs. e. 14 and E-15 into Eqs. F.21 and F.22, together with the 
requirement that the coefficients of the constant term, cos ty>-term and 
sin 4»~term should be equal to zero, yields the following equations: 



From the cos y-term: 



(F . 23) 


(F . 24) 
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From the sin ip-t era: 


-±(L) l l [(r'+ 3f>_ 1 )e <5 


(F . 25) 

Using Eq. F.20, taking (£/R) = 1 , A = 0 , B = 1 , m = const-, Eqs. F.23 
through F. 25 can be rewritten so as to be similar to the equations given in 
Bramwell (Ref. 20) 


OJr/ 2 %° - £ [ 0* (c< + { /*& ) - F* Xo 9<s] 


(F . 26) 


where 


[( c, *<$, i y - ) e« £ -y. A° g , ]*° 

-£j>[ ^ ^ ' 6 ~/f J J“° 

C, s 6, ? 

ft '? F/>9 £),=;< 

1 Z/ / ) * W 4 Vf* 


and 


F 6 Z l IF 9 

2 2 


(F .27) 


(F.28) 


(F.29) 
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From Eqs. F.26 through F. 28 


rv 


C, + G, 3 
v 


/ 



(F. 30) 


(F. 31) 



(F. 32) 


F.4 The Horizontal -Force Coefficient 

The horizontal force per unit span is given by the following relation 




(F . 33) 


where <j> = arc tan (U^/U^) = U^/U^. 

Using Eqs. 3.12, 3.13 and 3.7, the average horizontal force per revolu- 
tion can be written as 
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8 /IF 

nb 1 f / f [ VrCcto+a [U T Up6 -U P )J 

xir ^ 1 

- a. (6U T Z - UpUr) ff 4 ^/c eo+Jctfo df 


(F . 34) 


The substitution of Eqs. E.ll, E.12, E.14, E.15 and F.20 yields, after a 
considerable amount of algebraic manipulation: 

c^= i. cr ju Cct° (L j ^ V / A^ r /Cr/7 J 

- i $(i, 


(F.35) 


X8 JL9 

where F and F are defined in Appendix C. 

For consistency with the previous section, let 

ZiL/.o j 4 »o; «= ' 

U / 

and introducing the usual correction for radial flow effects, Eq. F.35 can be 
rewritten as 

Ch- C Cd.o 4 FS: £ 0 o Ac Ja. F l -hfy$ F*^ \oJ 

-3;{F n B' C -/*i;F»)ck 1F - 36 ’ 


, „19 „20 . , . „ 

where F , F are given m Appendix C. 
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F.5 Trim Calculation 


From Eqs. F.10 and F.36 


_ +0-H5 GJa.C4o -/■ 6a £0 o JLlF 1 ^ -f Q. ^ F ,7 ^ oj 

-i' ( F ' ,9 « -Mi fl °) § '■{(■£, j/‘"° 


In Eq. F.37, 0 is given by Eq. F .30 , is given by Eq. E.17 
Xo is given by Eq. D.10. 

This system of equations is solved by an iterative process on , 

computer. First 0, a are obtained and from their value 0, , 0, , and 
c r Is lc 

also evaluated. 


(F.37) 


and 

l digital 

X are 
o 
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APPENDIX G 


DEFINITION OF THE VARIOUS EXPRESSIONS, EQUATIONS, AND QUANTITIES 
USED IN THE PERTURBATION METHOD 


Note that almost all of the expressions given in this section have been checked 
using the FORMAC algebraic manipulative program. 

1/2 

Equation of 0(£ ' ): 


L j ( do, bo) - juuo Fo (Fo, Be) 

L 2 ( ao ibo)-JLcoHo( { Jo l &c] 


(G.l) 


Equations of 0(e): 

L 4 ( <*y, b< ] =• J*o { F 0 ^&c,do,bo) ( 'h t> , ®c) 

■+ J+4 Fo ('f' 0/ &c) -f Qi (&o, bo) 

Ll(a-I, &y) =r y 0 0 do, bo) -h Be } 

-/■ yh,y Ho('4o i Be) + A// (do, bo) 

3/2 

Equations of 0(e ): 

L 4 ( «2, M = J+o F 4 (%, 0 C} -b yUo F 0 [Girfo) + 


(G.2) 


('/o, «ojio)-f ( F’ofdo ) ~b/A'i Ft ('bo l Bc,do,bo) 


F 3 C Z - SinZ^o ^ 0 C l F 7 ^ o ^u, + 

-f C& 0 »L'/' o Qc J*oJ±4 C Z F^ -f Q z (do, bo, d-4 } b^ 


(-1 ( q Hi 6 / 0 , ®c, a i,^l) + JZo Hoi'4o / &i) -+ 

/±o Hz lF'o,0 C) cio, Fo) -f Kz ^%,doYju,H4(%Q c ao 1 b 0 s j 
-h iirt ZFo yAo J*-> F-l L + M z (do, bo, 0-4 J bt) (G.3) 
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(G.4) 


The operators L^, 1*2 are: 


L , fa, 4) = pk + ?JL + 


h'h 






Ll (a ' k) ?»20 i- 4, & (G.5) 

Q, (ac,kc) - c, a<> ^ + cM U F » feit+F* &» *? 1 (G.6, 

1>% \*/ L [ 2)%1 2>% lire J 

(ao,i>o) - - C 44 apdao „C}L [ f'djlo\ ll*+Q c L? 9 Ills, !tks 1 (G.7) 

■»% 7>Yi vft J 

Q, f a*, 4o, a,, 4«) = C 4 a, 34* Cu Go -J-2C, / V\& C F *dhi 34 0 

'aY'o T5’/ / o ‘dVb'Wo 

f C< iF'*[ ^ Mo J.JCI0 M± l+XllL e 

R L 'h'te, 'h'k -avk. J 


« L^% T>% 1% T>% 

-rt) 'd£o,Xu Mo _q Mo 7 - g °- 

L4 M% H ° 3d ‘° **+■ 




^ WoD*, 

-bfo T>% 

Mo di» 1 „ c„ ao Mi + Y, M° 

*£ J w. -n, -3% 


(G.8) 


i., -,M =. C 7 (l)[lL 11 $g ^ + 


'D'fe 

^<<0 _ « Ms. -2 C" V, 

+ \G 0 — $020 - 


c)qc 


(G.9) 
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(G. 10) 


Fo (V C* [ (zF%-FH x )sCnt, -^CosfJ] 

Ft-('k,» c }= C 6 [(zll *£ + \ocScL 2 )scnf 0 -t 

L ' / a 

(e c L , °-2L' , )>oc)$ 0 CoS'{' 0 J 

Note that starting with Eq. G.10, the term L 1 (C d )/a will be neglected i 
all of the calculations. ° 


F, ( < k } &c, <*-o> bo) - f ~ 29 c A Ci F" SC nfo 3 h> _ C,F 9 dao Sint 

R 7>% 

~ C 2 F Q a 0 cos% + ^.cos^ ~l (g. li) 


J< (Foj B c )--1 C z <j 0 F 7 S Cn lt 0 FI C 2 9 C F 3 cos +J - &cC x F 3 (G . 


12 ) 


H i ( 'i’Q t 0 C , a o, bo) = Ctj L 9 C ScVi^ - cHio 


7>% 






— c v L & c Cj 0 Ho cos FC 6 i1°cl 0 Q c CO s% - 
T>% 

2 \ oc C € L" o-o cos 


(G.13) 


AT; (Be, </■„)= l 2 ' C& Q c i [ Sin 2 ^ 

F^aoHo) = a e F ll, C z l d_bo cos ^ 

% 


(G.14) 


(G.15) 
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(G.16) 


J z (io ,flo)= '|c 2 F 7 a, Sen 2.^ 


/6 

H 2.(^0/ Be, a 0 ,bo)= -Zcos%a a CrjL dap _cos^ 0 a 0 c 7 L dbp {G11) 


T>% 


7>% 


K z (ao^o) -A. ®c a ° SCnlfc 


(G. 18) 


°^2 / — /. 

. *{<:<*)% ( 2oj ‘) 
ji 2 


Q (L)^ a) c (F l 3 Qc+ F Vr) 


- Co l r r Z - Crj ( | j tw t f L^r* °c L ' 9 ) J (G. 19 ) 


where 


M f A) - 


c -A z +L‘j t)l0 A - Xwo t A 

^>UO ~A ^^cIdiqA 


- XrO ^ A 


(G. 20) 


to 


I /C«c [ («,f- *>] 4- C, 4 C0c[2F% + F' S (^f)]J 

lulo L R 


- Crj f t\ (xJc 


Hereout- 
Pt = [-(ZCa>c + £ o(0 ) V ' + X- H ° ] 

p2 ~ ~ 2* tw c - ^£20 4- Kio F~ r J 

h-Fn^Vr+X-J^c 

Pn - [ ~*li c'Coc + Yt t-co c tT, J 


(u"te)4 e c L l3 [v, + (*,)*]} 


(C 


(G. 22) 

(G. 23) 
(G. 24) 

(G. 25) 


. 21 ) 
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(G. 26) 


Pis - c h j - ■+2*''cot e ti -t- Z ( v,) * J 

+ c < (^) Co ^ [‘<F l 3 e c &^ 2 ci l F' <r + ZF*fi z (V,)*] 

fa= C„ [i-<r,)*C0> c (-cti) - 2c/ 3 iu c v r J-C 7 (A'ju>o[^ i l 

+ (ir r )*ZA 2 » c L , *+2o< l e c L' 9 ] 




'M* 

(G. 27) 


(G. 28) 


(G. 29) 
(G. 30) 


(G. 31) 


(G. 32) 


(G. 33) 
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/>„= kz(v?fcv c F'*c<(A ) J -h a ol Co c F' !r Cj(tyCa OL cocC t/ + 

+ i b oz (V r f c H +3.Q C i> oz <*> c f= /J - 1 Cu> c j F ,q g o 

+ 1 fyrfcoc, f s - ± (V,f c z F c + ^ f"c 2 f£j e c (G 36 

p l2 - ~G C b oz ("*> ) LOc (A J C 7 L! 9 - e c a oz CV C I'AJ C 7 L 

— (tfr ) ^oc Cg L — J.(tK-) Co c 9 c L Crt -f- i C^F] Co c 4r F Cj 

l ' v° (G.: 

— I dot (tf r ) CV C C,i - t doz (^r) On — 2 &ol iV~r) CO c C>j L 

+ J- Coj c 9 c <j o L 16 c 7 + l(irf 6 C Cg l‘° 


- boz C ° c ^ C ’ — &oz Oo c F Ci ^ Aj-/- t'^oz ^r c ‘{ 

+ C Q,o 2 . < ^ J c c -<j - 2 9 C b oz co c F c, + d ( ' c ° c C z ^AJF J ^ 

(G.3 

~±Co c v,. C f F 3 - 1 v r C z F 6 _ Co c F y/ C 2 f£)&c 
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(G.39) 


jb ~ Q r b 0 J lO r Vy Cri ( t— -f Or ^Oj Cj c ^7 (— 1 
fin [a) [f>] 

$ 

_ V^. \oc Cq L -f- A ^ ^7 ~ t'^c. ^7 L- 

2. 


- 1 cioz Vb c u - L a °i V ~i' c " + 2 a 0 z co c Vy A c 7 L 


IS 


R 


to 


— AC (a)(_ e c L CCq 4- A “tfr Cq 
Z 2. 


f>jr] - - hoz Q-OZ C< f F -f L b B2 a aZ C H - O c hf>Z C I (-g)F 


j3 


Pis 


-h 4 CC Z F a j. J_ L b D i f A) Q —J. Got Cj F 

V ^ l do 2 

a oz c z FUjL Cz e c F*- Kz F"c z | 9 C 

2 Lf R 


jcp 2 2. f& / 

= 9c l>oz Goz c 7 (L\ L - t’<Z el Cl/ 4- Zoz (L \ C 7 C 


U 


_ a 0Z X oc CqL" _ ± i & c % L' 1 -h A Got & C L C 7 

Lf 2 

1 17 i I iG . . lo 

— L CtQl L Crj A °OZ &cg 0 Cry L 4. A ^6 — 


(G. 40) 


(G,41) 


ft 5 = C , (|j F 15 ]~ a 0 z(t,o 2 )4 t>02 (an)*] + 

2-%b„ (4«)*c, /^A) -f Cz/i\ P** [‘hz+fcba) J, G .42) 

10 ? y 2 \R l 

-h Cl /~ 9 [ GLo Z + ( Cioz)* J - ~ Cz F^ ^ a OZ 4-(a.oz) J 

4, Czdcf 3 -f C z f^) 9 c F [boil(l)oi) J F Cij fcboi (ttoi) -f{t boi) Q01J 
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/ 5 20 = 


- ^7 [ (b<n) <%t>2 + ( Cl ° z ) bozj'2 (xoz ) «roi C 7^j 

Cg /. /; A oc [ doz + (q-oz) ] - j- *-7 Q [ Cl o2+ M*] 

“ C 7 L '\ [( Ca °l)+ ( C *n)*] -$*-% L,Q9 * 

-/- 1 c 6 L l0 O c £ a 0 i 4- (aoif J 


if 



(G. 43) 



c 2 ^ j 0, F" [fan)* (4-coc) + 1>4 h(i + <*>‘)] + c <*c*,(, c h -b 
i hoL fan)* c 4 - 0 (^m) c v + 2^j #c ^ Cl J^bo-L ibil) (1-OJc) 


+ (bo Z ) 4 L<t (4+ Wc)J -Cctoz (b 45 )* Ci, (i-toc)4-Lfaoi)*lwC ¥ (i+COc) 

■bCy/f^.) F ~r CLoi ( by,})* [-i-CVc] -f- (O-ol) bji/ (i-bCt> c ) -h 
[fi. I L 

loi (o. a 3 ) (4-C Oc)+(bo2. ) °-w ('fbcoc)j * J- C<F 9 [fanf (I' 00 *) 
+ a 4H (nu>c)J~ Si F^ [(an)* +a^ ] +4- C L ^ F ILi j^t(^ 


~t" C bjtj (j (J'ftOc) + C f ( k 0 i) — boi j ^ t^c f f & 02 ) Ctoz'jJ (G.44) 

/^ 2 = c 7 Jz [(CUc-l) ( b+OOc) Qrjif J + Crj L L ^ ['('f~UJc)( c lil) 

( 1 +(*>c) a <V J -/■ C 7 L [-CL01 ( 1 + (Uc) Fr + (4-COc) Vr 
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where C 


* C,L Z i / <,< ’ 1 -htUlc) («.,.)* Q„] 

+ C 7 L% l [ (Ufa* + 

a °\ - ( i*u c ) ( a 02 fl,„] + C 7 J«2' Sc [ iti-co,) (i, 3 f 

c) -n] + Cj. L Q c V-r-Coilavt +(a„)‘) L J 

+ £U L‘°0 C j(a„)% a „ J _ c t L"\ oc [(a n )*+a„J - 

t " a ’“ V ‘ Vr C - C " "«'«•* (*»)*- c(««)X «„ c„ 



M~'(<) ]- c t[i<F‘<>c-UF>)t$.F‘ ] 

C ‘[- ‘ loc +%(& t L't 2 L" K c ] 


(G. 46) 


fr c < F ’<f (”»* - £jJl‘(v.)\ (±J f"q c w £ _ 

Cj (*) F LU>C + [~ U < Fco c -H/i, tv r )*J + 

C( ( |) w ‘ [ 26 , r > * *, F"*fr F’S(v r ) *] (G 47) 


t 6 ' C 7 L l-UJ c (V,) + -A C„(v r )*i: < ^ i ( C Oc-l) 

‘ C? /|)/ 2 ^ + «C i-'VjS, fvj)% 0 c lJ\oc,l 

~J C 7 L ’s c to c (V r ) % C„ ( vA* 

7 ^ G " 

is given by Eq. G.79. 
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(G. 49) 


F _ C z ( lF l & c -F i W) 

1 ol - 

a 


F °i-- 2*1 Si 

£ 


(G. 


Hoi- Cq \qc & c L 

£ 

Hol , A..) 3 . 


(G. 


(G. 


b - _ •/ tea. ■£ -/ i (W)*co 4 C,F 9 _ ± (\r r fc L F -h Co c F" c L t & c 

I ii 2 6 ^ °° 2 2 * 

P - = .(^f\, c c ( L".uv r )* ull e c L > c 1 4c(v r )* Cotgt L‘ 7 c,^e c?/ 1 


«. / /o 

fk) e c c e 


(G 


fa 


= iiW c c^ij <■<?*-£ v r CiF ri -u c F"ci(i')9 c 


t 


/V - 


// <7 

- lT r > 0C C 6 L" 4 ± bJiVrOcL C 7 - Cu) c Vr^ o ^n L 

- i. Cu> c 9 ^7 L- i- VJ- 

i <7o 2 


(G 




>7„ = ^ LC xF n 8 + r C ^ c ^ 

^7 rr do u 


(G 


/V s -# c * 


i'(?c 5f /- 


ll 


(G 


50) 

51) 

52) 

(G. 53) 
? 

C ? 

.54) 
(G. 55) 

.56) 

.57) 

.58) 
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(G.59) 


/>/$ = £z &c F ' 
2 

ho 


(G.60) 


(G.61) 


(G.62) 


hi- [l ^>45 F Q c + a >s Clf:9 -J Cz FG (««, + ^ is-) 

+ t fa £ a F' h 1 

00 J (G.61) 

fa - [-c 7 L QcCijs- - Xcdts cj^ L! 7 Crj ~l b iS 6 c L.'^ Crj 

+ l ^ A -c« L"\. c (a : s , a - e )] 

/ > 2S'= ^iw c c, (a, 3 - a ~ ) 2 C, F' l 6 c (l j cu c [ (uj c -i) i~ 

+ <,< ° cY i<v j * L ( 1 - u c) >>43 v rC H 4 i(*-fUc) C v t~ (V,)* 

+ C < CD fac-‘) V, b~ + (ucoc) (v r )*C, + 

K-*K 3 + «: 7„ «m)]-4 ^i ! * + ££ Ut f vl . (v fjc l £ F »J 

(G.63) 

^26 r { Zc ° cC^-^c) ^ «, 3 CrjL^i _ 2,aj t (l-HOc)(v,)*a~^ Crjll 1 f J_ 

R R 

4 C 7 (~j ® c c [ ( 1-Uk) a 7j - ( <+Q)c) an + ( i-a) c ) vr,fa _ 

& 4aic ) ky w] + 1 c 7 L^QcCtoc £ir r -(w)*] + 

2 «i8 c£*L - U h >'cCeL'_ l -r,Z ti c„-£( 1 r,)*a- Ctt J (G . 64) 


(G.63) 


(G.64) 
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hi- [(%-i) F 11 C, fo\ Q _ 

U 0C ‘ c a, s- C V^«< t-.s-i B c c, f ' 3 

+ + Z(ir,)*0J r , t P-/JT 

tVf '/ F +^<Z, s z,(l\F* 

4 I «7j C, /=■%, e-i) 4. j.C ( mfc F* - £ 

V 2 ~ Z /3 A -/ 

2 C 2 1 ij^ L - F' 1 * J 


(G.65) 




' /~ f ^-Wc) a, 3 0 £ z/c 7 , a ~ ^ c 7 ^'7 • _ 

‘/(V-,)*c U '. c r c„L‘ f /i) i ~ „ 

@}- 2w ‘ q 's£< 3 c C ',{- 


(« 

/* 


+ C 'W\*K Ca-ZlfK)**- Cj l 

[*u> c (fa)*F"c y £ 6 C _ c u c a, 3 c v / 4 Wf t7 3 l c,F 13 Q 

4 ^ S (+*> * c, F^j ^ lwr 


(G.66) 


~ W ‘(H' C >±3.^ + lc( V -rU c t l L F « 

2 « 


(G.67) 
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v e c c 1 (« l )*i'L' 1 c 1 + 

+ i(V,f l u>c c 7 l! 7 - 2 ( lT r f Ol e a H C n ll + 

2(t?r) (a) c trj j l~ — Cue 3/3 ^JLj @c ^7 ^ ~h 
^ C ®c <-7 L! 9 4- (co c -i)CUc( v,)* hi &ct<1 JL 3 + 

tU)c(^i) O'^CfjL -f> d L (Vs) CijcBc Crj L + i- 

-kf X>c ~~ i(v^) Clti^n-h M c G> 3 Cn 1 (G.68) 


K 4 r f><1TrY 6o t jo z X Ho 

Kt - - ( fytfr Y<ro f j>y ^-Hoj j K 4 

- i pit Vr %0 i pjk -A.H0 ) / Ki 
K, - ( pi i ^y YrO -b pn Y Hb 'j j K4 

«■ ~ ( p s V~ r Y$ 0 + p6 -&HO ) / ' & 4 
H-tf ~ ~ (v r Yfropj + Y-HO ho)/** 

^7 ~ ~ (V* ^tfo hi i Xrtopiy) / *• 

k $ ~ ~ (^r ^drO pit ~t X//(7 y£zC ) I k 4 

" f^T Y(irO pi<) -t Yho Plf'j j 'k>4 


(G.69) 
(G.70) 
(G. 71) 

(G.72 

(G. 73) 

{G. 74) 

(G-75) 
(G. 76) 

(G. 77) 
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(G. 78) 


K <o ~ ~ (hs)*t -^ H0 (ho) ]/*< 

c /2 = c 6 ( L lo O c -l\o C L j 



~h XhqL'^C ^”^7 L ^OC 



= - A- 

f Cuj c v r Y Gq , 

(-CiCv F* ±c,F") 

K< 

L 



~t X HO ^COc t~ ^ 7 L Oc ’^'y ~ £v~r Cr/ L J ^ 


, Z 
6)p,o 


(G.79) 


(G. 80) 


(G. 81) 



APPENDIX H 


THE SOLUTION OF EQUATIONS 5.101 AND 5.102 


For convenience, Eqs. 5.101 and 5.102 can be rewritten as 

(<tz+ A- * n )i + (*< -tu/Sh* (H ' 1> 

<L% 

where 

d*- ) * (K~ \oc) + n ( h . 2 ) 

d-X- (ft-ieijj (9-9 c) (Xo-Aoc) + /is K*I (H.3) 

and 

!„=/.*« +/* 3 K 7k ,„.4) 


^lo~ yU ^vx + Jj? fcyr 


(H. 5) 


The homogeneous system is given by Eq. H.l with £ =0 =0. The solution 

o o 

of the homogeneous system denoted by subscript h can be written as 


| A = d, e 5 '% + e Ss *‘ 

(H.6) 

ZX*r,e s '%z> 2 (T 2 e s ^ 0 
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where D 1 and D 2 are arbitrary constants and 


(T,= $4 - (di+/*? Kir) 

- <*Z 4 fL Kx)£ 


(H. 7) 


^2 = s 2 - ( <*<+ tc<m) 

~ dt -4 Jj? K/iX 

and S , S are the roots of the characteristic equation given by 
12 

s 1 - 2 4s+r^d 2 z +yi«,i z ]=o 


(H.8) 


(H.9) 


Therefore 


(ft* X,)n(^o~^oc)+^o k 6K 


(H.10) 


The particular solution of Eq. 5.103 is given by 


„ - I o(dr Kqi) 


(H. 11) 


7 l p - - 'Tjo ( dl+jL? k<w)+ %e (dl+/* L fyz) 

A 

A-( di-h/Sfyz) (d<- K<ff jx) +(d 2 +y Kqr) (dz'/Skil) 

The complete solution is given by 

i * h+ip 

Wh+I? 

so that from Eqs. H.6, H.13, and 5.100, the solution can be written as 


(H.12) 


(H.13) 
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Xj = f tit) = ft, (•»««) * % (l f *Hi>) 


- D. e^°+ ft, +■ 




(H.14) 


and 



-ti )£<& 

e 


£' /l 


(H.15) 
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APPENDIX I 


COORDINATE TRANSFORMATIONS AND PROOF OF THE EQUIVALENCE OF ROTATIONS 
1.1 Coordinate Transformations 

In this section, the various coordinate transformations for the various 
coordinate systems mentioned in Subsection 8.3.2 will be given, together with 
the various approximations involved. 

The relations between the i, j, k unit vectors and the i , j , k, unit 
vectors is shown in Fig. 29a. 



The relations between the i , ^ and the i 2> j 2> k 2 coordinate 

systems (see Fig. 29b) is given by 



The last two coordinate transformations are orthonormal. Thus, the 

T T 

inverse transformations can be obtained by using [R, (3)] and [R 2 (0,$)1 , 
respectively. 

If one assumes that the undeformed blade coincides with the i 2 direc- 
tion, the position of the elastic axis in the deformed state can be specified 

by two displacements v _, w _ [in the y_ and z directions, respectively, and 

e2 e2 2 a 

for simplicity u = 0] . These displacements for the same loading conditions 
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(1.3) 


will be related to v , w (see Fig. 27) by 

e e 

= V* +■ V0 e i-$) 

W ze = — + W e C Os(0--»-$) 

Then for 

tan-' fei) ; 22s?- W S 

\ 3 **/ «x, ) 

a coordinate system I , J , K can be defined. 

>%.X -a^X 

A coordinate system I , j , K, can be attached to the cross section 

•~1 ~1 -‘'I 

of the of the deformed blade at the elastic axis such that I, is tangential 

~1 


to the deformed elastic axis, coincides with the chord and is normal 

teen -fc ? 2J ^ 

V 3 **-/ a** t v*w 


to both I and J . For 

r>j X X 


LM 


F 


i 

^Vei. 

■aWei. 

ix7 


f \ 

u 

A/ 

j. 1 


_ ‘3 ,r e>i 

ax t 

1 

0 

j 

i* 

l*j 


?Vi ex 

__ 3x t 

' 

0 

1 

V 


rj 


(1.4) 


Also, for 3 small: 3/9 x = 3/3x_ = 9/9x . 

/ n 

In the inverse transformation, it will be assumed that 


i "fef + (^r) 1 s 1 
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Then the inverse transform is again given by using [R^] . 

Finally, the pertinent relations for the transformation between the 

i, j, k, I_, J . K . and I , J_, K„ systems defined in Subsection 8.3.2 are 
-v -v ~ ^2 ^2 v3 -v3 

given below 






\y J 


1 


3VI 

■ax 


. ?v Q»» 

1 ax 


- ^ 0 


-aw 

i>T 



f \ 


l 

r* 

- 

I 


/v 

^ J 


(1.5) 


f > 

Is 


10 0 

| 

h 

i J, 

3 

0 cm(S+-$) $no(e-+$) 

< 

T. 

V 

1C, 

v. J 


0 tos(&+$) 


K* 


(1.6) 


In the inverse transform of Eq. 1.5, again it will be assumed that 

, , 2 . 


f'aw 

ax 


+ + 1=1 

\lX/ 


(1.7) 


It is important to note that unless this assumption is made on the 
various coordinate transformations involving slopes, the nonlinear problem 
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becomes almost untreatable. 


1.2 Proof of the Equivalence of Rotations 

In this subsection, it will be shown that with the assumptions that 
0 , <J> are small and that 



(1.7) 


The final orientation of the blade cross section in space will be independent 
of the order in which the deformations occur. It will also be shown that it 
is immaterial if the rotations 0 and $ of the cross section take place about 
the feathering axis (i axis) , i axis or the I„ axis. 

«1 - s /2 


First, let us assume the following hypothetical displacement patterns: 


(1) The undeformed point on the elastic axis is located 

in the i. , j, , k system. 

X. J- 'v.L 

(2) The v . w deformations occur in the i_, j_, k system. 

e2 e2 «,2 ^2 ^ 

(3) The final orientation of the blade cross section in space 

is given by I , J , K unit vectors . 

^,1 ~-l .vl 

Measuring the final orientation of the I, , J, , K, system by the 

/wl ~1 ~1 

i , j, k vectors, the following relation can be written: 


/ ^ 



1 

r_ i 

:) 

r 

ii 

^ -B- 

n 

Ka 



'C' J 


l~J 


( 1 . 8 ) 


The transformation matrix [R^] , obtained by going through the various coordi- 
nate transformations and using Eqs. 1.1,- 1.2, 1.4, and 1.3, is given by 
Eq . 1.9 on page 

In the second stage, the following pattern of transformations is 
assumed: 
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( 1 ) 


A point on the elastic axis is assumed to have the 
displacements v, w (Eqs. 8.4 and 8.5). The new 
orientation of the blade is given by the I_, J_, 

/v 2 ^2 

K system. 

(2) The cross section is rotated by an angle of (0 + $) 

about the I_ axis. The final orientation of the cross 

section in space is given by the I„, J„, K„ system. 

^3 ~3 ~3 

Then the following relations can be written 




r 0 


i - - 

fo 1 ^ 


( 1 . 9 ) 


= ft, h 

if 


L J 

*v 


k 3 ' 


J* 

1 




J 


Performing the coordinate transformation, it is found that fR,.] is 
given by Eq. I. 10, page 

Comparing the two matrices [R^ ] and [R,_] , it can be seen that all of 

the elements agree, except R^(l,l), R 4 (2,3), and r 4 ( 3 j 3) when compared with 

the same terms of R_. The discrepancies between the two matrices always 

5 

involve only second-order terms which according to Eq. 1.7 are negligible 
compared to 1. 

By varying the displacement pattern similar results are obtained. 

This is not surprising since these are all the result of the assumptions 
3 small and Eq. 1.7, which is equivalent to taking 

e “( If) ° (Dr) s 1 

It should be noted that this approximation is related only to calculating 
the orientation in space of the deformed blade. This is required for the evalu- 
ation of the aerodynamic loading terms. It has not been used in calculating the 
inertia loads which therefore are consistently accurate up to O(e^). 
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APPENDIX J 


COMPLETE EXPRESSIONS FOR THE ACCELERATIONS 



(J.l) 


261 



In this report, whenever a x is used its meaning will be Eq. J.l without 
tiie terms marked by 


“j - * 


r + 2 SLv — J2? v + x z - (c os & - 


— ( £ i'n O' + $ - (s'\r> 8- + § cos *) $ - Si ( cos & ~ $ sir> ©■) 

+ 2 SI f ~ ( cos & - <§ s\r> 4- / jpr / 5 i^8-+ |cos £ 

1 c V pX 0 ^ ' 


+ f siTi e- C05 0-) § - / sinfr 4- $Cosfr) 

-3X. V X ' \ J 




'iVl 

1X a 


( Cos8 ^ § jJ 


(J.2) 


O.^ = w - X x ( sin 8 + $c/s*'j *■ 

+ X x ( cos 8 - $ s\ n ) 8 +■ X I ^ cos fr - § s /^ ) $ 


(J.3) 
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APPENDIX K 


DERIVATION OF THE MOMENT LOADS PER UNIT SPAN ABOUT THE FEATHERING AXIS 


K.l Derivation of q., 

xl 


In order to derive the quantity q^, the quantities and q^ are re- 
quired. From the geometry of Fig. 28 


i 


= M. + 


>r w ( H ’- w W, r ( K^J 


(K.l) 


where moments in the nose-up direction, for the cross section are considered 
positive. 

From Eqs. 8.11, 8.12, 8.20, and 8.23 


l 


■ 


( V + S\r)& + $^ + 2110 -Sl\) C 0 S&J X z 7n 


X Z/ m 


, $ si-T) 2 9 - SL* ^Cos&si’r)^ - 
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— cosfr sVw& — 2&1> SiprJ^COS^- — ^ si'v\ x & 


+ $ Cos& siryd) -f 2. SI ( COS0- Si'wG’) & 

^ ^ )^0 /- 

+ 2.Ji<§& —y^os^si'ndj +|^— 


4> cog^&'J + si^^ose- - & + 2Sl$ Qj^s'&j 


+ < - / *»»X 


- $si->>9) — 'mx.^ | ($ co.s x & - §><jT ^i^y^cos& 


+• 0 f COS X fr - 5 



) 


(K.2) 

The terms in Eq. K.2, marked by arrows, are of 0 (e^) , according to 
considerations (1) through (6) Subsection 8.4.1, and therefore are negligible. 
Thus 

f * M a - x .Sl* ^ si‘h& ■*-$cos^>' +- 'mx J . £'G'^swr>fr -*-<j5cos0^ 

- WC0S& + 2ft0(svn9-+i£coi:£j J 

— + SL x s.\v>Qcoifr 4'Jl a $cos7®- + si-n^ cos5 

+ ^ SiVfl)! - Z-Il§ — ■ COS 2^ 

-ax^t -1 3x>at 

+£e ^e%^dcose-2il^s^9-|~c«s^)^sin^] . 


(K.3) 
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Due to the model selected to represent the blade cross section, terms 
due to the polar moment of inertia of the cross section itself are missing. 
The missing terms are : 

(a) The regular angular inertia term given by 

(K. 4) 

(b) The missing part of the propeller moment effect, derived 
separately in Section K.3 of this appendix, given by 

-mjkJtffsinfrcosQ- (K. 5) 


From Eqs . K.3 through K . 5 

~ pj a - w X x fl Z (si*©" +§COS©J V + 7T) X x S? [> ^Sir> fr 

+ § cos - W cos & +• 2 5 ( & + ■£ c 0 S&- )J 
- vr\ [k* + X*) § Sl~ - nrn ( Jz 0 l + frcosb- +$ cos 



- vnX*Sl Z f2 ( shi 9- c.os& +- s-ivi x & ) - 2 §coS 2-©l 

* L \2X 0 J 

w o 

- x/ )Vil t - B- 2 S\r>frC0$& 

- 2 (f£ sin *-lc cos ^* s ''"*J 

Next, will be derived. From the geometry of Figs. 27 and 28 

(K.7) 
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where a positive q^ is one which gives compression in the nose fiber of the 
blade. From Eqs. K.7 and 8.18 


z = ttj x x £l [t? ^ cos & - J - fx D +■ ^ si-n O' J 

- U ( cos© - - 2u| CO 5 $ - ) j 


a a 

-f nmx x Sl 


'IX.'St’ 


cos 


x & — Sm6"tos6' 4-^0: COXES' 






^bL sVn9'coc©- > ) +2$ sivt&CGS&l + f~ §> { £;■* frees O' 

^x 0 ) J L [ox/zy 2 - 


_ -2!*L COff^©" 
'ZX 0 'd'f^ 


- c«*fr) 
•JX^t 7 


■+- cos^fr - 22L si'wO'Cos©') + 2, 5 f — -— ~ si-»>8cas8 

2X» 2X a / \^x 0 3f 

+? ( if si ' B8cojft - t£ twl9 ) ] 


3X« 


3x e 


r smfrcosfr - cos*®") & 4- 

Ux/aY 

'dX^Y / 


* / 

cos z fr^ 

+ 2 8/ .£«>¥» 8 C.o£ 9" -h $ coi 3. 

( 

\ 'ZX&Y 1 ' 

si<n9 - cofi8’ — t 

■2xpY^ J 


* -1 r 

— 2-S1 + | — < 
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The quantities marked by the arrows 
fore negligible . Thus 


are 0(£ D ) 


or higher and are there 


% 


= - ThSl x 


x £(x.+o+e,)c«e- - (x„+e,) 


+ HOX, 


- U COS& - -zv-cos&'j 


+ 



2 Xco s* 9 - 

^K c 


■f- QL si'n&cos® ~ 


DX a 


+ 2 .$ s\v> 9 " C-OS frj 


-g 3 v- c o5 l 9 — vJ — 


-h 


f%( z s'iyi&cosS- 

I— \ 


~ '£aL COS, 2.^ 



-f $ Coi2.fr — 2$ sit\&cos9- + 2$ ( St , hfreo.S& 
x "•' x — ■ UxTY 


3X f 


•ax,, 


- 2J-. COS'-®) + $’ (2S si^cosS- - ^ cos*e-) 1 

9X a 9^f / V**° ^ X “ ^ 

4 [ f 2 2l*L sr*frcos& — 'L*±_ 6os 2 &) 9 

L V '3X e 3t / 

4 - $7 22L s\y>Vcos9 - co5 2 ^) + 29 (sn.*c«s» 

V^x,, rx „ / v 


-f. <|> 


»} 


(K.9) 


In order to evaluate q , 6q^ is required. In this analysis, 3 is a 

small quantity, again of 0 (£ ) and when evaluating 3q , using Eq. K.9, many 

D z 

terms can be neglected. The result can be written as 


<j 2 ji s - vnSl*X x £ (*, + X 0 )pcosPj — 2 W! X_Jl X VjiCos& 

+ ”-*X [ ( + It. Sl,,ft “ se ty 


- '2j£T C os^fl - — — A £ir>$COS& 7 

/ w/av* r -J 


+• 


0 ( e»*) 


(K. 10) 
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K.2 


The expression for can be obtained from Eqs. K.6, K.10, and 8.35. 
Derivation of Mp 

From Eqs. 8.17 and 8.18 





pf> - - ~ (x.**,)/* - 2^] 


- *.X ACU* 

1 / 2X/3T7 


_ 2lbL a ( cos& + — si'nQ-) 1 

MM 1 -* I [ VXo > \ 


(K. 11) 
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From Eqs . K.ll, 8.23, and 8.24 



v L * - yv\& ^ w v* + $ vco&J - 'vr>x x & *si cose- 


- e,)^3v 


(K. 12) 


From Eqs. 8.20 and 8.21, neglecting u3 = 0(E d ) 
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-ffr* + zf§)(w-^X„) COS& _ 8*^Sil>& + 


+ *(£*». 


-| Zcos&'jt^-px.) 


* ^Tj~ / . ~ 

+• 2-i & — COS 9-( \AJ -px o J 


(K. 13) 


From Eqs . 8.33, K.12, and K.13 


H p * uL 

* t\ 


2. /* £ \ 

•V»J1 WV +^f 'U-COS O' ) 


- wx x W Xos$ - onJI ( / X 0 + e,^V — L ^ ^WJ-^X 0 J 

+■ 'm(yi-/3X 0 )(i? - v-Jsi 2 ~t 'mx^Sl j - C0£d-(w~^xj 

+■ $ £\r)&^ -^Xo') — ^iri'n^VM - jZX^ 




2('vJ-/3X.)f — let - 2-2 si«fr) 

' ' A 9X.9T 3X/jy J 


9X„3Y 3X^y 




7 


(K. 14) 
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In Eqs. K. 11 through K.14, all of the terms marked by arrows are 
negligible 0(£p) quantities. 


K.3 Derivation of the Propeller Moment 




FIG. K.l GEOMETRY FOR PROPELLER MOMENT EVALUATION 

The geometry of the problem is shown in Fig. K.l; f| is the chordwise 
coordinate (see Fig. 28) . 

Assume that the mass of the blade is compressed in a plate-like struc- 
ture with p(x,y) denoting the mass per unit area. 

From the geometry, the centrifual force acting on an element of mass 
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is given by 


<Lf = j=(x ; ^)fyS^ [C (Lx J 


Its component in the y-direction is 


<*F + SL'fr^) si^dLUrt 




The propeller moment per unit length of the blade is given by 

c(.M PRaP = ~J j > ( x ) r [) *}* otrj^x 

t-€ 

l.-e 

_ __J mrj t cos(& *§) ( & +• 

t.-e 


Using 


and 
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J 


'm - / k l 


,*■ Jr? + ?n X z 

* E 


t.e 


COS + •Sl’vi ^ + =^COS& - $>sir) ^ sivff’ *" §^ 0 ^& J 

~ si-^^cos^- - § + § cos’-fr - <5^sn>^co^6- 

clM p Rop = “ j^Sl'Vj&COS^' 4-^Cos2.aJ 

- m ^ £si'*9'C0£&- +$cos 26 -JJ 


(K. 15) 
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(K. 16) 



Comparing Eqs. K.16 with the appropriate term in Eq. K.3, it can be 
seen that the last term of Eq. K.16 is included in Eq. 3 while the first 
term of Eq. K.16 is missing and must be added to Eq. K.3. 
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APPENDIX L 


THE ELASTIC COUPLING EFFECT 


The elastic coupling effect is due to the angle of pitch setting 6 

which causes the bending perpendicular to the hub plane to be coupled with 

the bending parallel to the hub plane. In order to be able to represent the 

stiffness of the blade by rotating flap and lag frequencies, which is the 

usual practice in rotary wing-stability analyses, it is necessary to replace 
2 2 

cos 0 by (1 - sin 0) in Eq. 8.72. The appropriate expressions in Eqs. 8.72 
can then be rewritten as 

(ei) =(n) 3 * [ (ei) i - (ei) ? ] si,V (l.i) 


(ei)j W* HEi) a * -[(«), s ™ 1 ' 9 'M Er ^ (I " 2, 


It is convenient to define 


E ci 2 [H-ter),] s '" 1 9 

and 

E et = [( E1 ), -( er ) ? ] s "' 9cos9 


When applying Galerkin’s method in order to eliminate the spatial 
variable in the flap and lag equations , it is necessary to evaluate the 
following expressions 


- 






i: *«• 




(L.5) 



zJ fa '') 1 <*«. 

h a ‘ 


£ A 




(L. 6) 
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Note that throughout this report it will always be assumed that the first 
elastic mode in flap is the same as in lag (i.e., y^ = * thus Eq. L.6 

is satisfied. 

In rotary wing work, it is useful to express E^, E c2 by rotating or 
nonrotating first flap and lag frequencies. 


From elementary consideration 

. cr * k 




PIN ’ 'FI 


fJ (Et V {*')* *' « W * Mu 

r t a* L '“ 


(L.7) 


(L.8) 


Using Eqs. L.7 and L.8, together with the relation y = Eqs. L.5 
and L.6 can be rewritten as 


E e , - n c , J - “ f ,n J " M u -<C. ] «••») 

In order to be consistent with the approach used in this report, 

Eqs. L. 9 and L.10 must be expressed in terms of the rotating flap and lag 
frequencies of the blade. The rotating flap and lag frequencies of the blade 
will be calculated using Eqs. A. 2 and A. 4. The use of these equations implies 
that the rotating flap and lag frequencies of the beam are defined as those 
which would occur at a pitch setting of 0 = 0. 


Equations A. 2, A. 4, L.7, and L.8 yield 






it 




(L. 11) 
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r_J_[ __ 

+ Xj, nL ‘ ” 


(L. 12) 


With the assumption that = and Eq. M.6, Eqs. L.ll and L.12 can 


be rewritten as: 


— - z — i B 

W r,„ = W r, - 1 


Ft 


(L. 13) 


- a - r B‘ 

^UN = W U +!-■=- 
Thus, Eqs. L.9, L.10, L.13, and L.14 finally yield 

— Z 1 


OJ 


~ 1 — ul 1 = Ui - U) -f- 1 
lin Pin t-i Ft 


(L.14) 


(L. 15) 


and 


E tl * M F , si'n*d[w L * - u) p * f (l. 16) 

e~ cx * r%* si^cosej^wj -cr^ f ij * n u , si* 0 cos&[w L * (l.i7> 


ri . .. . . • Q ^ □ a ^ , the elastic coupling effect 

Using the approximation sm0 = 0, cos0 = 1, 

can be rewritten as 

il, - n tl »‘(<o L r -H>1) = F u •'[aj- =c + 1) (t - 18 ’ 

E Cl - M pl - UJp, +1^ - IL-19) 

Finally, it is important to note that the treatment of the elastic 
coupling in this appendix implies the use of the so-called Southwell coeffi- 
cient for rotating beam problems (Ref. 40) . As pointed out by Brantwell in 
Ref. 20, this treatment can lead to considerable error, depdning upon the 
assumed mode shapes in the flap and lag degrees of freedom. 
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For the mode shape defined by Eq. 7.1, it can be shown that the present 

treatment of the elastic coupling is a good approximation for u) > 1.2. 

FI 
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APPENDIX M 


- DEFINITION OF- THE- ADDITIONAL -GENERALIZED MASS INTEGRALS AND 
INTEGRALS ENCOUNTERED IN THE AIRLOAD EVALUATION 


In order to evaluate the various integrals required for the applica- 
tion of Galerkin's method on the flap and lag equations, various integrals 
which represent generalized masses must be considered. These expressions 
are defined in Eqs. M.l through M.14 below. For the feathering equation, 
Galerkin's method is not used; still, spanwise integrals similar to those 
occurring in the flap and lag equations are encountered. These are defined 
in Eqs. M.14 through M.18 given below. 

1 4 

B' [f +? .) <L *° (M - 1) 

k ° *o 



i / [ I -m (?,+?,) Ax, 

j-b 0 V 




(M.2) 



(M.3) 


(M.4) 


(M.5) 


(M.6) 


(M. 7) 


279 




(M. 8) 


{M. 9) 


(M. 10) 


(M. 11) 


(M. 12) 


(M. 13) 


(M . 14 ) 


CM. 15) 


CM. 16) 


(M. 17) 
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B 




iH nf f -W"**]-* 


(M. 18) 


The integrals due to the integrations of the airloads associated with 
the feathering equation are given below: 



(M. 19) 


(M. 20) 


T 3 



(M. 21) 
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APPENDIX N 


QUANTITIES REQUIRED FOR THE CALCULATION OF THE STATIC EQUILIBRIUM 
POSITION AND DEFINITION OF THE FLUTTER DERIVATIVES 


N.l Quantities Required for the Calculation of the Static 
Equilibrium Position 

The quantities associated with the calculation of the static equilibrium 
condition are defined below. 


S « ^ci 
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Q = - E 


(N.2) 
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(N. 8) 
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S u - 1(1 ♦<3.*)-fr(f*/* r V - I -I txj 


I ( 4 ) V ^' 2 


2 K 


C, - - &' * l 




F’6- - F z >, 
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(N.9) 


(N. 10) 
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(N. 12) 


F s „ ' B*Kf, 


(N. 13) 


L « = B Y*. - » , i*» 
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t £N - + fv^.T*.-M u *. to)' 
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N. 2 Definition of the Flutter Derivatives 

In defining the flutter derivatives , it is important to note that some 
are due only to elastic or inertia forces, while others can be due to both 
aerodynamic and inertia forces. In order to be able to keep track of the 
relative importance of the aerodynamic and inertia effects, the flutter 
derivatives in which these effects combine, will be split into two parts; the 
aerodynamic part will be denoted by the superscript A and the inertia part by 
the superscript I. 

The various flutter derivatives required for treating the flap equation 
are given below. 
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F,. = 
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The various flutter derivatives required for treating the lag equation 
are given below. 
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L*, = ( +E cl 
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T 2 UJ m m 


L 


I 1 . ft 

L- Y + 1_y 


3 


= 2- CO 


©2_ 


Ll 


*? 

1 Cl I 


K 


.17) 
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Finally, the flutter derivatives associated with the feathering 
equation _are~given_below_: 


T v = - {^T - 2 B 'Vf r - -Pr 


(N.18) 


I E ■ I - Txf 


(N. 19) 


T t ■ S ' z s:(nr) in - 20) 


V *H b3r « T ‘ + *^)*C' L '^’- 2T 7v* 


T » T 1 + T a 

» vp 1 y 1 Y 


K E = i(i*w.*)-T t 


K. - * 


(N. 21) 


(N.22) 


(N.23) 


(N. 24) 


t*. = r s ( 67 * g,l M) 


(N. 25) 


T s, = r E [‘ 1,a fr ( Ptf) ■ ( "- 26 ’ 


tj, w^h.)-zkk jj 

+ f ati { k h*- ^ d) ar* f *< i %'[ 


(N. 27) 
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x, T i, + T i, 


(N. 28) 


T*. = ^[b 7 ^ 4- B ,l J 0 (p+f r )+M F1 2J*S, + B» ( s (»+§.)] 




__ __X -r 

T = i + T 

X, ‘-X, ■>, 


x = i ( _g >r - : *y) 

t *L - 5^ - 

T- A = ~ f (!) tx „{ '• 

i E l / ^ — 


(N. 30) 


(N. 31) 


(N. 32) 


(N . 33) 


[>(»♦*.) w 

■x [£f± If >o r(e* $,)+z%r ] $; + 1 (froS+tf F'°) $A 


Ti = T* +• T* 

\ *, X 


T \ ' U [ B " + *1} ’ ' * Ml ' k ' + S ' c f J 


(N. 35) 


(N. 36) 
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T x " = - J- X (if f L“ ($■ * §„) - L' >. +• L " h,2°] 

z 2 I. UR L l “/ J 


(N. 37) 


T. = Tv 1 + T, 




* -2 r^J ♦ fli T*(l -*-,)( 5 ?. 5 ) j 


+ fQ-z,)b[-L'IJ-F‘f% -rfye 

+ X (t‘pp + f j ? ;) + 1 5 > f z ]-( fyeL«+ g y) K-L't;) 1 } 


?« = - T ? 
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APPENDIX O 


ALGEBRAIC EXPRESSIONS FOR THE VARIOUS QUANTITIES USED 
IN THE NONLINEAR EQUATIONS OF MOTION 


t « 5 ? 




B 7 X, 


(O.l) 


( 0 . 2 ) 


= [-n u ( x,2 K ♦ Z) -M'X ft. ' » F , fo'\ ♦*;)] f E <°- 3 ’ 


f ,l ( \ , r, *) " £ j * xsy ♦ t ? m X,i 


'FI 


+ F 




1 


NS 


( 0 . 4 ) 


L ^KK*i,Kr.r) m jjr {- g ’ 2 *.r -2(/i,), li x I x i 


- b ’X,Y> + |(|f[|fX,L r - |L->. y x, *^) 

-(sM*,)* + it i-x.ji’t'fx,]* 1 -™} 


( 0 . 5 ) 
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T Hl ± ^ M L) Y-M L ,h 4 0 AX z 'f + {& r 7 X 4 *S 

- ^ X^JFFrU ^ ^ ^ x/ ^ ^/)7 ^ _s/z / 2 / >_ ^ y 

- Ml, - * Lt x z[*i($o+r)4 - 2{v* p>f>*<*4 -h&x-lizfo 

+ Mfi *-,%<( p+(if>) + B° ^;x ( ^ f %, & n + 8'Y A y % v- 

+ p& r uZ 1 Y + % %*>{**[- ^) Z (^) z L i5 ^ 2 F ,/ x*X z 

-t % L 1 (x^4* z < f )] + fe(4-x»)(x~»-os)(-^fZ z y )]J 

- t (|) 2 ^Vt * 2 '($)[* (6+&+Y)L*-lo L ' 3 ] ** *2. 

_ L &y> 1 /V /. i *z.} + £ ft) V A *,/3 p Fy __ £ A 0 y> £ z A/3/> + K F V X */ 

/? j A l ^ / t /? /? _/ 


' £ ( 1)7 ( l'*) u '* 1 Y + ( F<e - FlXo ) x ’Y fL " , '?L*i'f +*’■?) 

+ ^-i % cK(r$oV (^ */) r]- il i r F '($°->r)'f *■> - %[ fM F> 


+ <<■?. F g (i^'f)] + L 1 * + (l-xjjiy [?*($/*, 

•h *,£ + ) — £.(hf FXz'j — iz^lpGL 

* 5ft)7 |*V/ F '' %*'*’[ F e (e+&r) -^of 3 ]-^° F ‘'^l° 
- L [kf" (6+$,+y)+ *■ C M F '°l*'Xi - il x 2'(+ x %'f)[t i i’ e ^‘ l 

ft L a, J R. 

+ (l*+ x <) eF '°] ~ %(* z $o+r^) x * eFlo +(%-fy l >H 

4 i ] - L 6 Tti ?c z [ {lf> F +($"■*■ Xi) L jJ + 7 ~ 5 N^} 


(0.6) 
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APPENDIX P 

COEFFICIENTS OF THE CHARACTERISTIC EQUATION 


The coefficients of the characteristic equation obtained by expanding 
Eq. 10.6 are given below. These coefficients were obtained by using the 
FORMAC algebraic manipulative system. 


V - t 5?, l Y f v, - L y, \ ~ L * r * \~ L ", f y t ¥ F y T s; + 1 


(P.l) 


* U + L *f _ L y,'h» F *x\ ~ F r : 


t 5T, L * \ ~ L J, R x-. \ L y F J, L ", \ ~ L f T ”, 


L l, F ? \ ’ 4* \ L y \ L y f y t s. 


F f T i* ” F f \ f Y t *. 


(P.2) 


^ + F x, v L x 1 +co --' ~ L f T vXr L r \3» 


M ] -^3- ~ L~T~ Q - L~ T*+ F 

®I>1 ?J>2<Td 1 T ^2 Fl y 
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*" x, % 


L x* K e F S; " T x, L ? F f, ' T i. L f F “ 1 


T-* L. T"«* L» Fw* 
X, Y * x X, x x 


- q F* - T* U*F« 

X, Oj>3 %j *1 T X l 


T v, 4 F i, ~ 4,4 T v l L ” 4 ~ 4,4 4 "4 


444 " 4, 4 T v. " 4*44 ~ 4 "4 


■ L i, 4 \ ~ 4*4 4 ' 4* 4, ~ 4*, F ? T *. 

+ }„J„ * F r 4 4 ' *-*, 4 4 ■ 4 F f 4 


F *4. - 4 4 4, - 4 4 ~ 4 4 


(P.3) 


3 s - 1 r** F - L * T> \ - L« T ®l +• J. 5 J 

a Y X 2 d-p» y \ <J p ’ Y Xx «*' <Ti>t<ipz<> 


D3 


+ K xL " 4,4 + 4 }- ~ L f 


-4V" + 44 *3 


“ A w a 


■f x 


j pi 


Ft 


"444 


l« T* F M F v M F v - L v ^ % 

Vp % t y^\CPS 
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L ♦ K F«* 
*. * X t 


T L* F*« 1— F«* ~ L*s F* 

*• Y x 4 *, t _ x t x, <r P3 % 4 


' L *x, * e \ - \ L ?V T i, L ?\ " T *, L -t\ 


l OJk~ 


- L i. F *> ' L x, F ? T K - l £, f y \ - V; ' L x F y T j, 


' L ?\' L i F T^ 


L ", F, f T x t ~ L X,. ?ra + W >-' + * F ? T X, L X, 


+ F ? T J, L X, - Q U Tj. - f„ F v T.. - c3 t * F** T. 
" L, T x, ~ *> Ti, ' T x, ~ ^ \ 


° l L r T x 1 } w L ?V., \U- + td| -' Wi 

+ L*4 - L,T«m* - L. T^w* - L„T X *;, 


+ 1 1 CO 2 -*- K OJ 2 - ) +W 1 W 1 1 -L F*F 

Cyi 0 V3 n e FI ^X* fi FI l' V X 2 X, 
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L * r ifx7' MrV*.' + l»x 3 F x, + V F *r 


W u F'x, ^X.^E F X t T X l L 'f F r 2 


L X, K e \ T x, L r F ^ " T x, L 'f \ L x,5' M F x l 
L X, ^E F x t ~ T x, l ? Fx i " T *i L £ Fx z ~ ~^x, L r F x* 

L xSi " L X, F v> T x t ~ L x, F y> 4^"^,“ L r\ 


L X, F ^ \ “ L X, F y\ t L X r K E * 


n *• 


LI 


+ F y T %: L x, + F f T i, L x l + F v T *, L *, -^‘vy 
-hT*Tj ~IX\ ' °u l| >VU> T *. 


-FT 

X, 


(P.5) 


*. = ^^4., - L x, l< J > , +*u K el, ~ L >\< 


L f \ < ' L xX 3 ^' + W » K ' ^ + 3 >* K < ^ 
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- L %5jv, 


-L t T jk F * 1 


+IAV L *.‘^ ' X, ’ ? 


- *L L vp F x x l x.1m Xl - 


PZ 


- r _ T i_,F -T?tF x l % F y T *-* 

- t - 5 ,K* F *> *■ * X ' ' 


- U F, 


.T„ - 4 Vh * 


♦ F..T* L, *■ F $>\ L *» 


X,' ? '*1 


- Eof.Fs.T- 


(P.6) 


-* p t, - 1 F. T x ~ Wur T 
* U u ' f X, T ’ 


= * u x,^ F ' 


-* - u. K.DJ, * w u K e ujri L 'f T »i F, ‘' 

-u k f„ ♦ 


- T,, 


K * 'x, 


- L F T- + Vx f l x, X x 

X, 'P X x 1 


F-T, W ~ w “ *f r% 


(P.7) 
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APPENDIX R 


EVALUATION 'OF' THE EXACT DIVERGENCE BOUNDARY IN FLAP-PITCH 


The exact divergence boundary can be obtained from Eq. 9.22 after 


suppressing the lag degree of freedom by taking h^ = 0. 
Then from the first of Eqs. 9.22: 


*->// 


From the last of Eqs. 9.22 


S 3i ( J + ^33 ~ ~^ N + ^3 


(R.l) 


(R.2) 


Using Eqs. N.15 and R.2, one obtains, after some algebraic manipulation, 

a third-order polynomial for <& . Thus 

o 


( — K m f, ) $° 3 - Jr M "t) 


+ ^ j _ S'jjJ § o + j' Cj_ C,J = O 


(R.3) 


where 


J- a 'V/V) + £( zf Fl>i ° + 


(R.4) 


Equation R.3 can be easily solved using Subroutine DPRQD available 

o 

in the IBM Scientific Subroutine Package. Thus, 4> ,g., are known. The exact 

o 1 

divergence boundary is obtained from Eq. 10.4, and can be written as 


= 0 - £ u d SsiJ ) - S 1$ ( £ - 3 / 5V 


(R. 5) 
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